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EDITORIAL 


Los efectos del cambio climático global se están dejando sentir cada vez con mayor inten- 
sidad en el mundo. Las inundaciones cada vez más frecuentes en Asia, Europa y América 
alcanzan dimensiones casi apocalípticas y sus consecuencias tienen efectos a más largo 
plazo que antes. En gran medida, estas inundaciones se deben a la fuerza descomunal 
de los huracanes, ciclones, tormentas, tifones y demás que azotan a las diferentes 
regiones. En fechas recientes hemos sido testigos de grandes tormentas que se originan 
en el trópico y terminan de descargar su furia en regiones tan septentrionales como 
Nueva Escocia o Canadá. Ante estas graves situaciones, los seres humanos tenemos el 
enorme reto de proponer y poner en práctica nuevas medidas de mitigación entre las que 
la reforestación es una de las más importantes. Otra medida importante que está siendo 
atendida con gran ímpetu en algunas regiones del mundo es la denominada construcción 
sustentadle. Este tipo de construcciones aspira a proporcionar una nueva forma de cons- 
truir que sea más amigable con el ambiente y que tenga un menor impacto en el uso de 
recursos. Las tres maneras en que se puede atender este punto son, en primer lugar, el 
desarrollo de nuevos métodos de diseño, buscando el uso más eficiente de los materiales 
estructurales, por ejemplo, propiciando que su respuesta estructural sea bajo cargas 
axiales, más que a flexión, que requiere por naturaleza de secciones más robustas: otra 
manera sería a través de nuevas formulaciones de los materiales de construcción, 
usando menos elementos contaminantes y sustituyéndolos por materiales que emitan 
menos gases de efecto invernadero y que generen menos residuos tóxicos o dañinos 
para el ambiente: finalmente, la tercera forma sería a través de un mayor uso de mate- 
riales “verdes”, como la madera y el bambú, principalmente. El uso de estos materiales 
apoyaría la captura y almacenamiento de carbono, reduciendo la presencia de ese 
elemento en la atmósfera. 

En este contexto, en este número incluimos varios artículos de investigaciones tecnoló- 
gicas que prometen aportar mejoras a los usos y métodos de diseño y fabricación de 
elementos estructurales con base en materiales lignocelulósicos. El primer manuscrito se 
refiere a las técnicas de decisión multicriterio usadas en la selección de componentes 
estructurales a partir de la tecnología de la madera, para construcción de viviendas 
sociales en Venezuela. Haciendo uso de la Metodología del Diseño Ambientalmente Inte- 
grado (dAi) y de las Técnicas de Decisión Multicriterio (mcdm), los autores de este trabajo 
llegaron a consolidar el desarrollo proyectual y conceptual de un nuevo componente 
constructivo con calidad estructural, bajo criterios técnicos de innovación tecnológica con 
el uso de materiales alternativos y nuevos procesos de fabricación. Como resultado, se 
logró seleccionar una opción que cumplía con los requerimientos de diseño en referencia 
a costo, calidad, resistencia y los principios del Ecodiseño. 

El siguiente artículo trata de la influencia de las características y técnicas constructivas, 
en las propiedades mecánicas de vigas laminadas encoladas de madera de Pino Caribe. 
Dicho trabajo es una investigación centrada en valorar las potencialidades de esta 
madera en la elaboración de vigas laminadas encoladas con calidad estructural. Este 
estudio permitió determinar la posible influencia de las características de crecimiento y de 
las uniones longitudinales “dentadas” fingerjoint en las propiedades mecánicas. 
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Otro manuscrito estudia el contenido fenólico y acción antioxidante de extractos de 
acículas de cuatro especies de pino de Durango, México. Se practicaron extracciones 
sucesivas, con acetona y metano! acuosos, de acículas de pino desengrasadas previa- 
mente con hexano. Se realizó una purificación cromatográfica a los extraíbles acetónicos 
generando cinco fracciones. Los ensayos con los extractos mostraron la supresión de la 
oxidación de las diferentes especies. Así se pudo determinar cuáles de los extractos 
tenían la mayor actividad inhibitoria. Los extractos acetónicos de dos especies presentan 
potencial como fuente de polifenoles bioactivos. 

Otro trabajo trata el tema de la variación estacional de compuestos fenólicos foliares en 
Quercus síderoxyla en diferentes tipos de suelo, como parámetro de deterioro del género. 
Se encontró para todos los suelos una disminución del contenido fenólico foliar en prima- 
vera y un aumento de éste en verano. Se determinó una correlación alta y positiva entre 
el contenido fenólico y la conductividad eléctrica en suelos de textura arcillo arenosa y 
también una correlación alta y positiva entre la retención de humedad del suelo y el conte- 
nido fenólico foliar, independiente de la textura. 

El siguiente artículo se refiere a la anatomía de la madera de las especies arbóreas de 
Ipomoea. Aquí se describe la anatomía de la madera de ocho especies de la serie Arbo- 
rescentes del género Ipomoea, grupo con porte arbóreo o de enredadera leñosa y flores 
de color blanco o amarillo. Los resultados coinciden con lo publicado para el grupo, pero 
se presentan diferencias estadísticas en la longitud y el diámetro medio de las fibras, en 
la altura y el número de células de los radios unisonados y en la altura y número de series 
de los radios multiseriados. 

Finalmente se incluye un trabajo sobre la modelación espacial de área basal y volumen 
de madera en bosques manejados de pinos en el estado de Hidalgo, México. Se mode- 
laron estas variables en rodales coetáneos de Pinus patula y P. teocote utilizando varia- 
bles de porcentaje de cobertura arbórea, índice de área foliar en sitio y alométrico, varia- 
bles cartográficas y de refiectancia obtenidas de una imagen multiespectral del satélite 
SPOTS. Se integró un Sistema de Información Geográfica (sig) constituido por varias 
capas de información. Con datos de campo y del sig se ajustaron modelos de regresión 
lineal múltiple. 

Esperamos que este interesante y valioso conjunto de resultados de investigación aporte 
datos importantes para que el uso y manejo de los recursos forestales de Latinoamérica 
se haga de manera más eficiente en el futuro. 


Raymundo Dávalos-Sotelo 

Editor 



EDITORIAL 


The effects of global climate change are being feit with increasing intensity in the worid. 
The increasingly frequent floods in Asia, Europe and America have reached almost 
apocalyptic dimensions and its consequences have longer term effects than before. To a 
large extent, these floods are due to the enormous strength of hurricanes, cyclones, 
storms, typhoons and others that plague the different regions. Recently we have 
witnessed great storms that origínate in the tropics and finish off discharging their fury as 
far north as Nova Scotia in Ganada. Given these serious situations, human beings have 
the enormous challenge of proposing and implementing new mitigation measures from 
which reforestation is one of the most important. Another important measure that is being 
treated with great emphasis in some regions of the worid is called sustainable construc- 
tion. This type of construction aims to provide a new way of building that is more environ- 
mentally friendiy and has a minor impact on the use of resources. The three ways you can 
address this point are, first, by developing new design methods, seeking the most efficient 
use of structural materials, for example, encouraging their structural response under axial 
loading rather than in bending, which requires by nature, more robust sections; another 
way would be through new formulations of building materials, using less polluting 
elements and replacing them with materials that emit less greenhouse gases and gené- 
rate less Waste, toxic or harmful to environment and finally, the third way would be through 
greater use of “green” materials such as wood and bambeo, mainly. The use of these 
materials helps the carbón capture and storage, reducing the presence of this element in 
the atmosphere. 

In this context, this issue ineludes several papers of technological research that promise 
to make improvements to the uses and methods of design and manufacture of structural 
elements based on lignocellulosic materials. The first paper deais with multicriteria deci- 
sión making techniques in the selection of structural components from wood technology, 
for construction of social housing in Venezuela. Using the methodology of Integrated Envi- 
ronmental Design (ied) and the Muticriteria Technical Decisión method (mcdm), it was 
possible to achieve the successful consolidation of design and conceptual development of 
a new constructive component of structural quality with technical innovation criteria with 
the use of optional materials and new manufacturing processes. First, alternativas to the 
structural elements were generated with the creativa brainstorming technique. Secondiy, 
a selection technique was carried out according to the environmental, socio-economic and 
technological change, which carne to define one alternativa. Using specialized software, 
it was chosen an option that met the design requirements in terms of cost, quality, strength 
and Ecodesign principies. 

The next paper deait with the influence of the growth characteristics and fabrication tech- 
niques on the mechanical properties of pinus caribaea timbar glued-laminated beams. 
This paper is a research, focused on assessing the potentials of this wood for the produc- 
tion of structural quality glued-laminated beams. The study allows determining the 
possible influences on the mechanical properties of glued laminated beams, especially the 
finger joints. 
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A third paper presents the results of a study on the phenolic contení and antioxidant acti- 
vity of extracts from needles of four pine species from Durango, Mex. Successive extrac- 
tions were performed with acetone and methanol aqueous pine needles previously degre- 
ased with hexane. Chromatographic purification was performed on acetone extractadle 
generating five fractions. The tests with the extracts showed the removal of the oxidation 
of different species. This could determine which of the extracts had the highest inhibitory 
activity. The acetone extracts of two species have potential as a source of bioactive poly- 
phenols. 

Another paper deais with the seasonal variation of foliar phenolics in Quercus sideroxyla 
in different soil types, as a parameter of deterioration of the genre. It was found for all soils 
decreasing foliar phenolic contení in spring and increasing in summer. A high and positiva 
correlation between phenolic contení and electrical conductivity was determinad in sandy 
clay soil textura and aiso a high and positiva correlation between the retention of soil mois- 
ture and leaf phenolic contení was found, independent of soil textura. 

The following article refers to the wood anatomy of tree species of Ipomoea. It describes 
the wood anatomy of eight species of the series Arborescent of the genus Ipomoea. This 
group has arboreal habit or woody vine and flowers white or yellow. The results coincide 
with those published for the group, but there are statistical differences in the length and 
average diameter of the fiber, the height and number of uniseriate ray cells and the heíght 
and number of sets of multiseriate rays . 

The last paper ineludes work on spatial modeling of basal area and volume of wood in 
managed forests of pines in the State of Hidalgo, México. These variables were modeled 
with similarly aged stands of Pinus patula and P. teocote using variables of percentage of 
tree cover, leaf area Índex and allometric site, mapping and reflectance variables obtained 
from a multispectral satellite image spots. The authors integrated a Geographic Informa- 
tion System (gis) consisting of several layers of information. With fieid data and gis linear 
regression modeis were adjusted. 

We hope that this interesting and valuable set of research results provide important data 
for the use and management of forest resources in Latin America is made more efficiently 
in the future. 


Raymundo Dávalos-Sotelo 

Editor 



Madera y Bosques 16(3), 2010:7-22 
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lARTlCULO DE INVESTIGACION! 


Las técnicas de decisión multicriterio 
en la selección de componentes estructurales, 
a partir de la tecnología de la madera, 
para construcción de viviendas sociales 

en Venezuela 

Multicriteria decisión making techniques in the selection 
of structural components from wood technology, 
for construction of social housing in Venezuela 

Wilver Contreras-Miranda'', Vicente CloqueII-Ballester2 
y Mary Owen de Contreras^ 

RESUMEN 


Las metodologías de diseño requieren aplicar técnicas que permitan, en la fase de generación 
de alternativas, generar métodos que logren seleccionar la mejor de éstas, y con ello, la satisfacción 
de las principales necesidades de los distintos requerimientos del proyecto. En el sector forestal, el 
caso de la industria mecánica forestal de la madera laminada encolada no escapa a estas exigencias. 
Haciendo uso de la Metodología del Diseño Ambientalmente Integrado (dAl) y de las Técnicas de 
Decisión Multicriterio (MCDM), para la etapa mencionada, se llegó a consolidar con éxito el desarrollo 
proyectual y conceptual de un nuevo componente constructivo con calidad estructural, bajo criterios 
técnicos de innovación tecnológica con el uso de materiales alternativos y nuevos procesos de fabri- 
cación. Primero, se generaron un total de 53 alternativas de elementos estructurales con la técnica 
creativa tormenta de ideas. Segundo, se realizó una selección técnica según los factores ambientales, 
socioeconómicos y tecnológicos, con la cual se llegaron a definir siete alternativas. Haciendo uso del 
software Expert Choice 2000 para MCDM, se logró seleccionar de las siete, la número 39. Esa 
cumplía con los requerimientos de los factores señalados en cuanto a costo, calidad, resistencia y los 
principios del ecodiseño, tales como: menor consumo de energía, uso de materias primas alternativas 
como las gramíneas, reciclaje y reutilización de materias primas, entre otros. 

PALABRAS CLAVE: 

Ecoeficiencia, Ecología Industrial, desarrollo, innovación, materiales alternativos. 


ABSTRACT 

The design methodologies require appiying techniques in the phase of generating options, selecting 
methods to achieve the best of them, and thereby meeting the needs of the various requirements of the 


1 Universidad de Los Andes, Facultad de Ciencias Forestales y Ambientales, Centro de Estudios Forestales y 
Ambientales de Posgrado, Laboratorio Nacional de Productos Forestales. Mérida, Venezuela. wilver@ula.ve 

2 Universidad Politécnica de Valencia (UPV), Departamento de Proyectos de Ingeniería e Innovación, 
Valencia, España, cloquell@dpi.upv.es 

3 Universidad de Los Andes, Facultad de Arquitectura y Diseño, Escuela de Diseño Industrial, Mérida, Vene- 
zuela. marowen3@hotmail.com 
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Las técnicas de decisión multicriterio, en la selección de componentes estructuraes 


project. In the forest industry, the case of 
engineering glued laminated timber is not 
immune to these demands. Using the 
methodology of Integrated Environmental 
Design (lED) and the Muticrltena Technical 
Decisión (MCTD) method, referred to the 
stage, it was possible to achieve the 
successful consolidation of design and 
conceptual development of a new construc- 
tivo component of structural quality with 
technical innovation criteria based on the 
use of optional materials and new manufac- 
turing processes. First, a total of 53 alterna- 
tivas to the structural elements were gene- 
rated with the creativa brainstorming tech- 
nique. Secondiy, a selection technique was 
carried out according to the environmental, 
socio-economic and technological changa, 
which carne to define alternative 7. Using 
the Expert Choice 2000, it was chosen 
among the last 7, option number 39. This 
one met the requirements of the above 
factors in terms of cost, quality, resistance 
and Ecodesign principies, such as reduced 
energy consumption, use of alternative raw 
materials such as grass, recyding and 
reuse of raw materials, among others. 

KEY WORDS: 

Eco-efficiency, Industrial Ecology, develop- 
ment, innovation, alternative materials. 


INTRODUCCIÓN 

Las Técnicas de Decisión Multicriterio 
(mcdm) se han ido transformando en herra- 
mientas importantes en la actualidad, 
especialmente en la alta gerencia de las 
organizaciones, asesores y consultores 
técnicos para empresas privadas y guber- 
namentales. Su fin ha ido trascendiendo 
en la medida en que los requerimientos 
van siendo cada vez más complejos, ya 
que su naturaleza permite, con base en 
criterios diversos, poder establecer 
patrones de comparación y, de manera 
metodológica, científica y tecnológica, 
llegar a una decisión consensuada para la 
mejor alternativa. Sus campos de aplica- 
ción son muy variados, desde la selección 
de la localización de una industria en un 
territorio determinado, hasta llegar al 
complejo campo del diseño industrial para 
la selección de un mejor concepto de 
diseño (Ramon-Solans y Ferriz, 2006). 


En ese contexto, el presente trabajo 
tiene como objetivo principal la aplicación 
de las Técnicas de Decisión Multicriterio 
(mcdm) como una herramienta metodoló- 
gica de gran importancia en la toma 
consensuada de decisiones importantes 
para una organización o proyecto, que en 
este caso particular se aplica a la selec- 
ción de un componente estructural, tipo 
viga, a partir de la tecnología de madera 
laminada encolada; todo ello enmarcado 
dentro del contexto del Diseño Ambiental- 
mente Integrado (dAi), que según 
Contreras y Cloquell (2006), es análogo 
al ecodiseño y, en este caso particular, al 
diseño industrial de un producto forestal 
de alto valor agregado como lo es el 
diseño de un producto forestal de 
madera laminada encolada con calidad 
estructural (pfml), con fines sociales, 
como la construcción de edificaciones 
prefabricadas para familias de menores 
recursos. 

En la actualidad no se puede negar 
que uno de los grandes retos de la 
sociedad moderna es procurar la equidad 
en la forma de vivir de todos sus ciuda- 
danos, especialmente los de menores 
recursos económicos, a los cuales se les 
debe procurar vía los gobiernos o admi- 
nistración central, regional y local, así 
como de las empresas privadas que 
conforman la industria de la construcción 
e inmobiliaria de un país, quienes deben 
aportar soluciones al grave problema del 
déficit habitacional. Según Nájera (2009), 
América Latina es un recinto geopolítico 
donde resalta el perfil urbano de sus 
ciudades con los contrastes entre las 
edificaciones y urbanizaciones opulentas 
y los barrios o chabolas que circunscriben 
las montañas y sitios inhóspitos para vivir. 

De ahí que la aplicación de las 
Técnicas de Decisión Multicriterio en el 
diseño industrial y arquitectónico permite 
abrir nuevos horizontes técnicos, más 
certeros en el cribado de alternativas 
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(Romero, 1993), donde un proyectista (dise- 
ñador industrial, arquitecto o ingeniero), 
debido a la falta de herramientas o de la 
selección de alternativas, y si su decisión 
final es errónea, puede llegar con ello a incu- 
rrir posteriormente en la pérdida de recursos 
financieros hacia la empresa contratante 
(Aragonés et al., 2004). 

Asimismo, se puede reflexionar a 
partir de lo dicho por Barba-Romero y 
Pomerol (1997), acerca de que la difusión 
de estas herramientas técnicas, o el uso 
de programas especializados de aplica- 
ción, los cuales facilitan los procesos de 
selección, como es el caso del Expert 
Cholee 2000 u otra versión más actuali- 
zada, es seguro que propiciarán un 
mayor intervalo de calidad, seguridad y 
satisfacción entre las partes involucradas 
en la selección de la mejor alternativa por 
desarrollar, esto en cualquier ámbito de 
acción del hombre moderno. 

MATERIALES Y METODOLOGÍA 

Contexto de trabajo de la investigación 

El contexto geopolítico y tecnológico de 
aplicación de la presente investigación 
fue Venezuela, país que en la actualidad 
tiene un déficit habitacional estimado en 2 
millones de viviendas (Fundación Vivien- 
da Popular, 2009). Siendo las más afec- 
tadas las familias pertenecientes a los 
estratos socioeconómicos C, D, E y F, 
que de acuerdo a lo expuesto por Venam- 
cham (2007), las clases socioeconómicas 
en Venezuela tienen ingresos estimados 
en dólares a cambio de la moneda 
nacional por 2.150 bolívares (Bs) por 
cada dólar americano, es decir: clases A- 
B ingresos (4.693,37 dólares); clase C 
(1.039,00 dólares); clase D (493,00 
dólares); clase E (277,80 dólares); clase 
F menores ingresos a la clase E. Para el 
año 2010 esta realidad cambia de forma 
negativa con la devaluación anunciada 


por el gobierno nacional y el cambio de la 
moneda a Bolívar fuerte (BsF) de 2.600 
BsF por cada dólar americano. Según 
Camacho (2007), eso estratifica a las 
clases sociales Ay B de la sociedad vene- 
zolana a la categoría alta. 

El uso de sistemas constructivos son 
variados con tecnologías que van desde 
el tradicional de hormigón o concreto 
armado de vigas y columnas; estructuras 
de acero con mamposterías de bloques 
de cemento y arcilla; hasta, quizás el de 
mayor proyección en el último decenio, 
sistema túnel de concreto armado. Este 
último, por su rapidez de fabricación en la 
edificación de unifamiliares, bifamiliares y 
multifamiliares, representa la alternativa 
más competitiva respecto a los sistemas 
más tradicionales y comúnmente em- 
pleados en Venezuela para la construc- 
ción de viviendas sociales. 

Recientemente el gobierno nacional 
ha venido proyectando, desde el año 
2000, el uso masificado de las planta- 
ciones de pino caribe {Pinus caribaea var. 
hondurensis (Sénéclauze) Barret y Golfari) 
de la Orinoquía. Este recurso estimado en 
más de 600 mil hectáreas, representa la 
masa arbórea más grande del mundo 
concentrado de una sola especie forestal 
(CVG Proforca, 2008). De igual forma, 
existen grandes propuestas de uso de las 
principales gramíneas como el bambú 
(Bambú sa vulgaris Schrad. ex J.C. 
WendI.), la guadua (Guadua angustifolia 
Kunth, 1822) y la caña brava (Gynerium 
sagittatum (Aubl.) P. Beauv.), en la manu- 
factura de productos forestales con calidad 
estructural para edificaciones con fines 
sociales (Cloquell etal., 2004). De ahí que 
un pequeño grupo del gremio de arqui- 
tectos, diseñadores industriales e inge- 
nieros, de manera individual o institucional, 
a través de centros universitarios de pres- 
tigio como la Universidad de Los Andes, la 
Universidad Central de Venezuela, o la 
empresa estatal CVG Proforca, por medio 
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de la gerencia del proyecto Vivienda con 
Madera, estén abocando sus esfuerzos 
técnicos para aportar la mejor solución de 
prototipos de edificaciones con madera 
sólida, mezclada con otros productos 
forestales y materiales tradicionales para 
la construcción. El plan de CVG Proforca, 
en convenio a 25 años con la filial petrolera 
Petróleos de Venezuela (pdvsa), es a 
mediano plazo construir con madera un 
promedio de 50 mil viviendas por año 
(Contreras et al., 2007). 

Por consiguiente, se plantea en la 
presente investigación, dentro de la 
visión del Diseño Ambientalmente Inte- 
grado (dAi), una nueva proyección en el 
contexto de los productos forestales con 
calidad estructural como son las vigas y 
viguetas para pisos, entrepisos y techos 
del tipo forjado, a partir de una mezcla 
de partículas de caña brava y adhesivo 
fenol formaldehído. La misma vendría a 
mejorar el diseño del producto industrial 
anteriormente desarrollado por Con- 
treras et al. (1999), referido a los 
tableros aglomerados de partículas de 
caña brava y adhesivo urea formal- 
dehído (uf) con resinosidades R10% y 
R13%. 

Además, la propuesta alternativa 
viene a solucionar, entre otros proble- 
mas, las prestaciones de resistencias 
físicas y mecánicas de los tableros, su 
proyección de uso y muy especialmente, 
la disminución de la posible escasez en 
un futuro cercano de especies autóc- 
tonas de bosque natural, por las grandes 
presiones debidas a la protección del 
medio ambiente: razón por la cual es 
vigente la propuesta de utilización de la 
gramínea caña brava como materia 
prima para la fabricación de tableros 
aglomerados de partículas y demás 
productos forestales para ser usados en 
cerramientos y sistemas estructurales, 
disminuyendo con ello la presión sobre 
el bosque natural. 


Otra razón importante de la presente 
investigación es que toda la suma de crite- 
rios de fabricación de los productos fores- 
tales de madera laminada encolada, invo- 
lucra por igual al material lignocelulósico de 
madera sólida de pino caribe, mezclada 
con otros materiales, variedad de especies 
de madera y tipos de adhesivo que sirvan 
para el mismo fin estructural. 

Contexto metodológico 

Una vez aplicados los principios metodo- 
lógicos del Diseño Ambientalmente Inte- 
grado (dAi) en la iteración de desarrollo de 
alternativas, con la técnica creativa de 
Brainstorming o Tormenta de Ideas 
(Cloquell, 2003; Contreras y Cloquell, 
2006), se pueden desarrollar un total de 
53 alternativas. En una primera etapa se 
seleccionaron 33 expertos, que a criterio 
de los autores tenían el mayor prestigio e 
inherencia técnica y científica en la temá- 
tica tratada en el ámbito iberoamericano, 
dadas sus publicaciones de libros y artí- 
culos en revistas de alto impacto. Final- 
mente respondieron 15 expertos a los 
cuestionarios aplicados, lo que permitió 
interactuar e incluir sus respuestas en el 
procedimiento metodológico y la conse- 
cución para alcanzar los objetivos plan- 
teados en el diseño. Se establecieron los 
siguientes pasos teniendo bien claro el 
problema a resolver: 

1 . Claridad en los objetivos del proceso 
de diseño. 

2. Se designó a un responsable de regis- 
trar todas las ideas promulgadas. 

3. A fin de que los comentarios de otras 
personas actuaran como estímulos de 
sus propias ideas, en una especie de 
reacción en cadena, se designó a un 
responsable que asegurara cumplir las 
siguientes reglas: suspensión de 
juicios o críticas improcedentes; 
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aceptar todas las ideas y registrarlas; 
animar a los expertos a construir sobre 
las ideas de los demás; animar a que 
se expresen los expertos con sus ideas 
“locas" o “fuera de lugar” (Contreras y 
Cloquell, 2006). 

Siguiendo la secuencia metodológica del 
dAi, se llega a la Iteración del cribado de 
ideas. Finalmente los autores del presente 
trabajo, con la generación de la matriz de 
requerimientos de diseño del producto y los 
análisis de ciclo de vida según el método 
Acv-Coclowen de diagnóstico propuesto 
por Contreras y Cloquell (2006), proce- 
dieron a la búsqueda de la mejor alternativa 
de diseño de un elemento estructural para 
la construcción de vigas y viguetas estruc- 
turales para viviendas sociales. Final- 
mente, haciendo el cribado de alternativas 
por descarte técnico, pertinencia y niveles 
de innovación antes planteados en la 
matriz de requerimientos de diseño del 
producto y ACV-Coclowen, los autores 
llegaron a la selección de siete alternativas 
(figuras 1 y 2). A éstas se les definieron 
cada uno de los criterios de valoración a 
partir del uso de la Técnica de Decisión 
Multicriterio (mcdm), método de Analytic 
Hierarchy Process (ahp) propuesto por 
Saaty (1977) (1980). 

La aplicación del ahp fue por medio 
del uso del software Expert Choice 2000, 
las figuras 3, 4, 5, 6, 7 y 8, exponen de 
forma detallada los resultados obtenidos. 
La aplicación del Expert Choice 2000 al 
proceso de selección de ideas, partió del 
análisis consensual sobre cada una de 
las comparaciones en cada alternativa. 

La designación de pesos de impor- 
tancia o preferencia de alternativas se 
realizó con la opinión consensual de dos 
expertos, donde se obtenía una sola 
respuesta, la cual era finalmente introdu- 
cida en el software para el desarrollo de 
las matrices. Los niveles de inconsis- 
tencia se encuentran dentro o cercanos a 


los límites recomendados por Saaty 
(1997, 1980), llegándose a obtener las 
siguientes inconsistencias: 

a) Para los pesos de los criterios de 
importancia entre alternativas (0,05). 

b) Criterio de influencia del factor 
ambiental (0,12). 

c) Criterio de influencia del factor 
socioeconómico (0,08). 

d) Criterio de influencia del factor 
tecnológico (0,09). 

RESULTADOS Y DISCUSIÓN 

El desarrollo metodológico de la investiga- 
ción permitió obtener resultados finales 
para cada una de las consideraciones 
recomendadas por García (2004) para la 
participación consensual de los expertos, y 
dado que la mejor propuesta de diseño de 
PFML para el contexto geopolítico y tecno- 
lógico de Venezuela fue la Alternativa 39 
(Figura 2), seguida por la Alternativa 15 y 
Alternativa 42 (Tabla 1) (Figura 3). 

Se determinó que en la selección de 
la mejor alternativa de las últimas siete 
seleccionadas, existió una uniformidad de 
criterios de evaluación para cada alterna- 
tiva, prevaleciendo en el método ahp, 
según fuere la alternativa, y en orden de 
importancia el factor tecnológico (64,4%), 
socioeconómico (27,1%) y ambiental 
(8,5%) (figuras 4, 5, 6). 

La tabla 1 expone el orden de impor- 
tancia de la selección de las mejores 
alternativas con la proyección de mayor a 
menor en la valoración de los pesos 
finales obtenidos de la aplicación del soft- 
ware Expert Choice 2000. La figura 1 
representa las Alternativas 39, 15 y 3, que 
tienen, en la opinión de los expertos, un 
menor nivel de sensibilidad del factor 
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Tabla 1. Comparación de los resultados de la selección de alternativas en el proceso 
de “cribado de ideas”, por el método AHP de Saaty (1997, 1980). 


Orden de importancia 

Selección de Alternativas por el Método AHP, 
usando el software Expert Cholee 2000. 

1 

Alternativa 39 (0,193) 

2 

Alternativa 15 (0,167) 

3 

Alternativa 42 (0,155) 

4 

Alternativa 6 (0,155) 

5 

Alternativa 31(0,139) 

6 

Alternativa 21 (0,102) 

7 

Alternativa 3 (0,090) 



OHIU/l 
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PARTlCUL&S GRAMÍNEAS O MADERA ♦ 
MADERA SÓLIDA » ADHESIVO 
ESTRUCTURAL 


partículas gramíneas o madera + 

MADERA SÓLIDA + ADHESIVO 

ESTRUCTURAL ^ 


partículas GRAMÍNEAS O MADERA + 
TIRAS DE GRAMINEAS + ADHESIVO 
ESTRUCTURAL + LAMINA METALICA 


PARTÍCULAS SiAWiNEAS 0 MADERA < 
ADHESIVO ESTRUCTURAL •* REFUERZOS 
DE FIBRA DE VIDRIO, CARBONO. 
ARAMIDAO ACERO 


-flatoso 


8cm 


Figura 1. Representación parcial de las alternativas de diseño seleccionadas para un 
producto forestal industrial laminado encolado con calidad estructural. Se aprecian seis 
de los conceptos realizados con su numeración en concordancia con lo expuesto en la 
figura 3. Éstas forman parte del total de 53 alternativas generadas en la 
“Tormenta de ideas”. (Fuente: Elaboración propia) 
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ambiental. Ha sido mayor la sensibilidad, 
respecto al factor socioeconómico de las 
Alternativas 6, 42, 31 y 21, dado que 
éstas pueden presentar mayores niveles 
de riesgo por impactos negativos en la 
obtención de sus materias primas, 
procesos de manufactura y seguridad 
industrial de los trabajadores. 

Se observa que el programa 
empleado permite proyectar una visión 
clara respecto del uso y diseño de las 
relaciones matriciales, escalas de valora- 
ción, niveles de complejidad y representa- 
ción gráfica. Dado que el nivel de conoci- 
mientos de los expertos es alto en los 
aspectos científicos y tecnológicos del 
problema planteado, el uso del método 
AHP ha permitido que disminuya el conte- 
nido de subjetividad en la decisión de 
cada valoración de aspectos, respecto a 
las acciones y atributos de una propuesta 
de diseño, llegándose a obtener resul- 
tados excelentes y confiables. 

La formulación de la metodología del 
Diseño Ambientalmente Integrado (dAi) y 


de la aplicación de las Técnicas de Deci- 
sión Multicriterio, dan un verdadero 
soporte y consistencia técnica al producto 
final propuesto, definido como pfml 39, 
permitiendo tener de forma segura una 
mayor aproximación a una certera conclu- 
sión técnica, con sus fortalezas y debili- 
dades, así como la proyección y su facti- 
bilidad real de fabricación. 

De esta forma, el industrial que 
pretenda manufacturarlo puede llegar a 
emprender de forma definitiva el proyecto 
industrial de manufacturar el producto con 
sus respectivas previsiones de recursos 
financieros, de infraestructura, de tiempo, 
etc., o bien si tiene que recurrir a desarro- 
llar otro proyecto de diseño de producto 
forestal a partir de las seis alternativas 
restantes que fueron generadas en el 
proceso de la iteración de cribado de 
alternativas. 

Se debe inferir técnicamente la rela- 
ción entre la trascendencia que debe 
alcanzar el desarrollo del nuevo producto 
industrial pfml 39, es decir, un producto 



Figura 2. Representación conceptual de mejores alternativas de diseño propuestas a 
partir de la aplicación de la técnica AHP de Saaty (1997, 1980) en el software Expert 
Choice 2000, según lo expuesto en la figura 3. (Fuente: Elaboración propia) 
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Figura 3. Definición del Expert Choice 2000 de la estructura del árbol jerarquizado y los diferentes valores obtenidos 
para la selección de la mejor alternativa de diseño de un PFML, definido según los expertos, de los criterios 
y aspectos más importantes de cada criterio, con sus respectivos pesos de importancia. 





Figura 4. Definición del Expert Choice 2000 de los diferentes valores obtenidos y nivel de importancia 
de cada alternativa de diseño de un PFML respecto al factor ambiental. 
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Figura 5. Definición del Expert Choice 2000 de los diferentes valores obtenidos y nivel de importancia de cada alternativa 

de diseño de un PFML respecto al factor ambiental con su respectiva inconsistencia. 





Figura 6. Definición del Expert Choice 2000 de los diferentes valores obtenidos y nivel de importancia 
de cada alternativa de diseño de un PFML respecto al factor socioeconómico. 



Figura 7. Definición del Expert Choice 2000 de los diferentes valores obtenidos y nivel de importancia 
de cada alternativa de diseño de un PFML respecto al factor tecnológico. 



Figura 8. Definición del Expert Choice 2000 de los diferentes valores obtenidos y nivel de importancia 
de cada alternativa de diseño de un PFML respecto al factor tecnológico y su respectiva inconsistencia. 
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ecoinnovador, englobado dentro de la 
visión de la Integración Ambiental Total 
que procura el establecimiento del 
desarrollo sostenible global, no sólo en la 
sociedad, sino que la industria forestal 
venezolana y sus proyectos, procesos y 
productos, tengan una relación equili- 
brada y armoniosa con el entorno. 

De ahí la necesidad de retomar la 
temática ambiental como eje fundamental 
para el desarrollo del dAi. Porque todo 
hace suponer que el desarrollo del proto- 
tipo del elemento estructural laminado a 
partir de la mezcla de partículas de caña 
brava y adhesivo fenol formaldehído, 
cuya tecnología puede ser análoga con la 
de tableros aglomerados de partículas de 
madera maciza, permite vincularla, por 
igual, a la ingente necesidad de con- 
servar de manera sostenible el medio 
ambiente, lo cual es un hecho que actual- 
mente preocupa a personas de todo el 
mundo. Por ello, la ciencia y la tecnología 
de los productos forestales debe estar un 
paso adelante de las propuestas tradicio- 
nales ofertadas en el mercado. 

El movimiento mundial que se ha 
venido consolidando también ha llegado 
al sector forestal y todas las industrias 
que le conforman, tal es el caso de la 
madera laminada encolada. De ahí, que 
en la industria de puertas, como en otras, 
ya existen claras definiciones de incluir el 
factor ambiental a la producción indus- 
trial. La Norma une 56200;2003 ex (2003), 
expone la visión de futuro, lo que hace 
ganar a los consumidores sensibles al 
deterioro ambiental que demandan cada 
vez más información acerca de la inci- 
dencia o impacto que sobre el medio 
ambiente causan productos y servicios 
que consumen habitualmente, esto con el 
fin de seleccionar aquellos menos perjudi- 
ciales para el entorno. Ese es el fin 
primario de la Ecología Industrial, que 
como ciencia pretende hacer una 
analogía entre los ecosistemas naturales 


que se encuentran en equilibrio dinámico 
respecto de un sistema de producción 
industrial de un producto (Fiksel, 1997). 
Además, esta ciencia contempla entre 
sus pilares fundamentales la Ecoefi- 
ciencia, y ésta a su vez el Ecodiseño, el 
cual es el sinónimo español del Diseño 
para el Medio Ambiente (díe). 

Este contexto se centra en el 
proyecto, procesos y productos, mientras 
que el nuevo paradigma del Diseño 
Ambientalmente Integrado (dAi), como 
una herramienta de diseño, se contextúa- 
liza desde la concepción de la Integración 
Ambiental Total (iAj), haciéndolo más 
holístico e integrador con todos los 
factores que interactúan, como son: 
procesos industriales, economía, socie- 
dad, política, cultura y medio ambiente, 
en conjunción con el resto de factores 
que se deben considerar para poder 
alcanzar el ansiado desarrollo sostenible 
global (Capuz y Gómez, 2002; Gómez- 
Orea, 2002; Cloquell, 2003). 

Por ello, la participación de las 
Técnicas de Decisión Multicriterio en el 
diseño industrial de multiplicidad de 
nuevos productos forestales industriales 
que se vienen desarrollando y manufactu- 
rando en los últimos años en la industria 
forestal mundial, especialmente la de los 
países desarrollados industrialmente con 
sus novedosos proyectos de l+D, no 
escapan a este planteamiento filosófico. 

El Diseño Ambientalmente Integrado 
(dAi), haciendo uso de las Técnicas de 
Decisión Multicriterio, hace resaltar por 
igual, entre muchas de sus considera- 
ciones, que estas propuestas deben 
surgir de proyectos de l+D, donde ha sido 
positiva la interrelación entre la industria 
forestal y centros de investigación, que 
han laborado en conjunto para alcanzar 
objetivos tecnológicos y, con esto, el éxito 
económico y comercial. 
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CONCLUSIÓN 

Se aplican de manera exitosa las Técnicas 
de Decisión Multicriterio según el método 
Analytic Hierarchy Process (ahp) a la meto- 
dología de Diseño Ambiental mente Inte- 
grado (dAi), para el desarrollo de un 
producto forestal de madera laminada 
encolada con calidad estructural identifi- 
cado como PFML 39, realizado a partir de 
partículas de gramíneas como la caña 
brava, el bambú o la guadua, y proyectado 
por igual para ser manufacturado con las 
partículas de madera maciza de pino 
caribe, de las plantaciones forestales 
ubicadas al sur de los estado Monagas y 
Anzoátegui de Venezuela. 

Las Técnicas de Decisión Multicriterio 
aplicadas por medio del software Expert 
Cholee 2000, son una herramienta prác- 
tica, de mucha objetividad, consensual y 
con proyección de una buena certeza en la 
selección de alternativas de un problema 
determinado, siempre y cuando exista un 
alto nivel en los conocimientos científicos, 
técnicos y humanísticos de los expertos 
consultados. La presente aplicación es 
original dentro del sector de la industria 
forestal, sobre todo en Venezuela. 
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lARTlCULO DE INVESTIGACION! 


Influencia de las características 
y técnicas constructivas en las propiedades 
mecánicas de vigas laminadas encoladas de 
madera de Pinus caribaea var. hondurensis 

Influence of the growth characteristics and fabrication 
techniques on the mechanical properties of Pinus caribaea 
var. hondurensis timber glued-laminated beams 

Pablo Ninin'i, Wilver Contreras-Miranda2 y Styies W. Valero^ 

RESUMEN 


El presente trabajo es la culminación de un amplio espectro de investigaciones realizadas en el 
Laboratorio Nacional de Productos Forestales, Mérida, Venezuela, centradas en valorar las potencia- 
lidades de la madera de Pinus caribaea var. hondurensis en la elaboración de vigas laminadas enco- 
ladas con calidad estructural. Este estudio permitió determinar, a partir de la inferencia de ios análisis 
estadísticos, la posible influencia de las características de crecimiento y de las uniones longitudinales 
“dentadas" fingerjoint, en las propiedades mecánicas: esfuerzo en el límite de proporcionalidad (ELP), 
módulo de ruptura (MOR) y módulo de elasticidad (MOE) de dichas vigas laminadas. Se logró deter- 
minar, en orden de importancia, qué defectos son relevantes para estas propiedades mecánicas, 
llegando a conclusiones y recomendaciones muy concretas que llaman a la reflexión a los fabricantes 
venezolanos de vigas laminadas encoladas, para la mejora continua de sus procesos, productos y 
servicios, a fin de poder seguir manteniendo la oferta de productos estructurales seguros y de alta 
calidad en los niveles de resistencia mecánica. 

PALABRAS CLAVE; 

Análisis estadísticos, esfuerzo en el límite de proporcionalidad, esfuerzos de diseño, módulo de elas- 
ticidad, módulo de ruptura. 


ABSTRACT 

This paper is the culmination of a wide range of researches carried out at the Laboratorio 
Nacional de Productos Forestales, Mérida, Venezuela, focused on assessing the potentials of the 
Pinus caribaea var. hondurensis timber for the production of structural quality glued-laminated beams. 
Henee, the study allows determining, from the statistical analysis inference, the growth characteristics 
and fabrication techniques, the possible ínfluences on the mechanical properties of glued laminated 
beams, especially the finger joints, with respect to the stress at Proportional Limit (ELP), Modulus of 
Rupture (MOR) and Modulus of elasticity (MOE). From the results, the influence of those defeets can 
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Influencia de las características y técnicas constructivas... 


be determined in order of importance, 
which allows to get very concrete conclu- 
sions and recommendations to make the 
glued-laminated beam Venezuelan 
manufacturers reflect about the conti- 
nuóos improvement of their processes, 
Products and Services with the aim of 
keeping the supply of high quality and 
reliable structural products. 

KEY WORDS; 

Statistical analysis, stress at proportional 
limit, design parameters, modulus of 
elasticity, modulus of rupture. 

INTRODUCCIÓN 

La calidad de la madera como material 
orgánico anisotrópico, ya sea en su 
condición de madera sólida o como un 
producto forestal de valor agregado y alto 
valor agregado, caso de la madera lami- 
nada encolada, depende de toda una 
suma de factores intrínsecos y extrín- 
secos propios de la anatomía del mate- 
rial, niveles de calidad de los procesos de 
transformación para la obtención y manu- 
factura del producto forestal final, condi- 
ciones de uso y los esfuerzos a los que 
estará sometido el mismo (Contreras et 
al., 2005; Cloquell et al., 2007). 

Desde el punto de vista de sus 
defectos anatómicos y de las propiedades 
mecánicas, la calidad de la madera es 
evaluada en función de los defectos que 
ésta contenga, principalmente las grietas, 
rajaduras, desviación de la fibra, tamaño 
de nudos, acebolladuras, ataques de 
agentes xilófagos y otros más específicos 
como el grosor de la pared celular, la 
presencia o no de parénquima, presencia 
de extractivos o contenido de humedad 
(León y Espinosa de Pernía, 2002; 
Valero, 2001). La influencia de estos 
defectos, en el caso de la madera sólida 
o madera laminada encolada con fines 
estructurales, se considera mediante el 
uso de normas de clasificación que varían 
mucho, podiendo ser visuales (los más 
usados) o mecánicos destructivos y no 


destructivos (Smuiski, 1997). De ahí que 
el desarrollo de nuevos productos fores- 
tales con aplicaciones estructurales se 
orienta a la eliminación de la influencia de 
los defectos en la resistencia mecánica 
de la madera. La madera laminada, 
microlaminada o el tablero contrachapado 
limitan la influencia de sucesivos defectos 
(Ross, 1985; Ross et al., 1996). 

Por consiguiente, la proyección de 
uso de la madera del Pinus cañbaea var. 
hondurensis de las plantaciones fores- 
tales del oriente de Venezuela en la 
manufactura industrial, ya sea en 
productos forestales de madera sólida o 
productos forestales derivados, obliga de 
manera inmediata a la profundización de 
las investigaciones científicas y tecnoló- 
gicas para la mejora continua del diseño 
industrial de estos productos. 

OBJETIVOS 

Estudiar la influencia de los defectos o 
características de crecimiento de la madera 
y de la técnica de uniones longitudinales 
“dentadas” o de finger joint, en las propie- 
dades mecánicas de las vigas laminadas 
encoladas, fabricadas con madera de Pinus 
cañbaea var. hondurensis y adhesivo isocia- 
nato dimetil (mdi). A partir de este estudio se 
podrá llegar a inferir técnicamente la carac- 
terización, por orden de importancia, de los 
posibles factores negativos que más perju- 
dican en alcanzar los altos estándares 
exigidos según las normas internacionales, 
a los productos laminados con calidad 
estructural manufacturados en Venezuela. 

METODOLOGÍA 

La investigación se realizó en las instala- 
ciones del Laboratorio Nacional de 
Productos Forestales Sección de Ensayos 
de la Universidad de Los Andes (Mérida, 
Venezuela). La misma se desarrolló a partir 
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de los resultados obtenidos por Contreras et 
al. (2007), en la determinación de ios 
esfuerzos de diseño de las vigas laminadas 
encoladas de Pinus caribaea var. hondu- 
rensis y adhesivo mdi. La tabla 1 expone los 
resultados obtenidos del estudio previo y la 
descripción de las características dimensio- 
nales de las vigas. El presente estudio se 
centra en el análisis detallado de 18 vigas, 
de las cuales nueve fueron manufacturadas 
a partir del uso de la técnica de unión finger 
Joint identificadas como fj, y las otras nueve 
fueron fabricadas con láminas (tablas/ 
tablones) completas, identificadas como ec. 
La industria fabricante es la única de este 
ramo, hasta el presente, en Venezuela, 
encontrándose localizada en la zona indus- 
trial de Matanzas, de la Ciudad de Puerto 
Ordaz, Estado Bolívar. 

Una vez recibidas las vigas en las 
instalaciones del laboratorio, se procedió 
a su identificación con números que van 
de 1 hasta el 18 (Tabla 1). Posteriormente 
cada viga fue analizada morfológica- 
mente según los aspectos anatómicos 
(color, tamaño y localización de nudos en 
cada viga, dirección de la fibra, presencia 
de defectos por podrición); aspectos 
técnicos según lo definido por Freas y 
Seibo (1954), de la manufactura de las 
vigas (calidad física de las líneas de cola, 
aserrado-labrado mecanizado para la 
conformación de las láminas o tablas/ 
tablones de madera-uniones dentadas). 
Las dimensiones nominales de las 
láminas que conforman cada una de las 
vigas a partir de tablas son de 2,5 centí- 
metros y de los tablones 5 centímetros. 
Las dimensiones nominales de las vigas 
las exponen detalladamente Contreras et 
al. (2007). La determinación de los 
esfuerzos de diseño se realizó a partir de 
los ensayos de flexión estática de 
acuerdo con la Norma astm d- 198 - 84 . La 
información se reportó de manera deta- 
llada en planillas descriptivas y gráficas. 
El análisis de estadística descriptiva, 
análisis de varianza y análisis multiva- 
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riado de componentes principales se 
desarrolló a partir de los valores 
expuestos en la tabla 1 que fueron intro- 
ducidos, entre otros valores, en la base 
de datos desarrollada con el software 
Minitab versión 12. 

RESULTADOS Y DISCUSIÓN 


Estadística descriptiva de las propie- 
dades mecánicas de las vigas lami- 
nadas de Pinus caribaea var. hondu- 
rensis 

En las tablas 2 y 3 se observan las esta- 
dísticas descriptivas de las propiedades 
mecánicas de Deformación (mm) en el 
Esfuerzo al Límite Proporcional, Esfuerzo 
al Límite Proporcional (elp) (MPa), Módulo 
de Ruptura (mor - MPa) y Módulo de Elas- 
ticidad (MOE-MPa), tanto para vigas 
compuestas por elementos unidos por 
medio de la técnica de finger joint, como 
por aquellas compuestas por elementos 
completos. 

A partir de los resultados de la esta- 
dística descriptiva se puede afirmar que 
las vigas elaboradas con finger joint 
presentan menores valores promedio en 
las propiedades mecánicas Deformación 
(mm) en el Esfuerzo al límite Proporcional 
(MPa) y Módulo de Ruptura (MPa) que los 
promedios obtenidos en vigas elaboradas 
con láminas completas. 

Se destaca que el Módulo de Elasti- 
cidad (MPa) es ligeramente mayor en su 
valor promedio para vigas elaboradas con 
finger joint que para vigas elaboradas con 
láminas completas. Considerando la 
desviación estándar como medida de 
dispersión de los valores de las propie- 
dades mecánicas en estudio, se destaca 
que los valores son mayores para vigas 
elaboradas con láminas completas que 
para vigas elaboradas con finger joint. 



Tabla 1. Determinación de los esfuerzos de diseño a partir de la Norma ASTM D-1 98-84, de las vigas laminadas 

de Pinus caribaea var. hondurensis encoladas con adhesivo de MDI. 
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Tabla 2. Estadística descriptiva de las propiedades mecánicas de vigas laminadas 
compuestas por elementos unidos por fingerjoint. 


Propiedad mecánica 

n 

Promedio 

Desviación 

Estándar 

Deformación al Limite Proporcional (mm) 

9 

27,3 

5,6 

Esfuerzo al limite proporcional (ELP) (MPa) 

9 

17,61 

2,37 

Módulo de ruptura (MOR) (MPa) 

9 

23,73 

3,53 

Módulo de elasticidad (MOE) (MPa) 

9 

3112,4 

235,7 


Tabla 3. Estadística descriptiva de las propiedades mecánicas de vigas laminadas 

compuestas por elementos completos. 


Propiedad Mecánica 

n 

Promedio 

Desviación 

Estándar 

Deformación en el Esfuerzo al Límite Proporcional (mm) 

9 

27,7 

8,3 

Esfuerzo a! límite proporcional (ELP) (MPa) 

9 

19,21 

3,24 

Módulo de ruptura (MOR) (MPa) 

9 

30,42 

7,50 

Módulo de elasticidad (MOE) (MPa) 

9 

3376,6 

609,9 


Tabla 4. Análisis multivariados de componentes principales, variación explicada 
para la propiedad mecánica Deformación (mm) en el Esfuerzo al Límite Proporcional, 
considerando variables de manufactura y anatómicas en vigas laminadas (n=18 vigas). 


Aspectos de crecimiento y manufactura 

Valor propio 

Proporción 

Acumulada 

Número de Nudos 

2,5623 

0,512 

0,512 

Número de Uniones Fingerjoint 

1,7401 

0,348 

0,860 

Tipo de Viga 

0,4102 

0,082 

0,943 

Leño Juvenil 

0,2670 

0,053 

0,996 

Deformación (mm) en el EsfUerzo al Límite Proporcional 

0,0204 

0,004 

1,000 


Tabla 5. Análisis multivariado de componentes principales, variación explicada para la 
propiedad mecánica Esfuerzo al Límite Proporcional (MPa), considerando variables de 
manufactura y anatómicas en vigas laminadas (n=18 vigas). 


Aspectos de crecimiento y manufactura 

Valor propio 

Proporción 

Acumulada 

Número de Nudos 

2,4121 

0,482 

0,482 

Número de Uniones Finger Joint 

1 ,6775 

0,335 

0,818 

Tipo de Viga 

0,4963 

0,099 

0,917 

% Leño Juvenil 

0,3926 

0,079 

0,996 

Esfuerzo al Límite proporcional (MPa) 

0,0216 

0,004 

1,000 
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Figura 1. Relación de variables de manufactura y anatómicas, con respecto 
a la propiedad mecánica Deformación (mm) en el Esfuerzo al Límite Proporcional, 
mediante el resultado gráfico de análisis multivariado de componentes principales. 


Entre las vigas elaboradas con finger joint 
y las vigas elaboradas con láminas 
completas no existen diferencias estadís- 
ticas significativas a un 95% de confiabi- 
lidad, para las propiedades mecánicas en 
estudio. Esta afirmación se fundamenta 
en los resultados de los análisis de 
varianza de una vía, diferencia de medias 
de Tukey y Fisher, considerando como 
variable dependiente el tipo de viga y 
como variables dependientes la Deforma- 
ción (mm) en el Esfuerzo al Límite 
Proporcional (MPa), Módulo de Ruptura 
(MPa) y Módulo de Elasticidad (MPa). No 
existen diferencias estadísticas significa- 
tivas en las propiedades mecánicas 
causadas por el tipo de manufactura de 
vigas en estudio. Los mayores valores de 
desviación estándar, de las propiedades 
mecánicas correspondientes a vigas 
elaboradas con láminas completas, 
sugieren que el control de calidad, para 
valores admisibles de carga de refe- 


rencia, generará mayores dificultades que 
para vigas elaboradas con finger joint. 

Análisis muitivariado de componentes 
principaies de ias propiedades mecá- 
nicas de las vigas laminadas de Pinus 
caríbaea var. hondurensis 

El análisis de componentes principales 
(acp) es una técnica utilizada para reducir 
la dimensionalidad de un conjunto de 
datos. Intuitivamente la técnica sirve para 
determinar el número de factores subya- 
centes explicativos tras un conjunto de 
datos que expliquen la variabilidad de 
dichos datos. Se emplea sobre todo en 
análisis exploratorio de datos y para 
construir modelos predictivos. 

La variación de la propiedad mecánica 
Deformación (mm) en el Esfuerzo al Límite 
Proporcional en vigas laminadas es expli- 
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Figura 2. Relación de variables de manufactura y anatómicas, con respecto 
a la propiedad mecánica Esfuerzo al límite proporcional (MPa), mediante el resultado 
gráfico de análisis multivariado de componentes principales. 


cada en 51 ,2% por la variable número de 
nudos, 34,8% por el número de uniones, 
8,2% por el tipo de viga y apenas 5,3% por 
el leño juvenil (Tabla 4). En la figura 1 se 
observa la relación de variables de manu- 
factura y anatómicas, con respecto a la 
Deformación (mm) en el Esfuerzo al Límite 
Proporcional, donde se destaca que la 
propiedad mecánica en cuestión resultó 
con un comportamiento indirectamente 
proporcional con la variable anatómica 
número de nudos, siendo esta variable la 
que explica la mayor variación. Es decir, en 
la medida que se incrementa el número de 
nudos en una viga laminada los valores de 
Deformación (mm) en el Esfuerzo al Límite 
Proporcional se hacen menores. La posi- 
ción de los vectores de las variables de 
manufactura Tipo de Viga y Número de 
uniones, con respecto a la posición del 
vector de la variable mecánica en cuestión, 
expresa una mínima incidencia en su 
magnitud. 


En el tabla 5 se observa la explica- 
ción de la variación propiedad mecánica 
Esfuerzo al Límite Proporcional (MPa) en 
vigas laminadas. La variable Número de 
Nudos explica 48,2% de la variación de la 
propiedad mecánica en cuestión, seguida 
por 33,5% de la variable Número de 
Uniones, siendo apenas de 9,9% y 7,9% 
de Tipo de Viga y %Leño Juvenil, respec- 
tivamente. En la figura 2 se destaca un 
patrón similar al observado para la 
variable de propiedad mecánica Deforma- 
ción (mm) en el Esfuerzo al Límite 
Proporcional. 

En atención a las propiedades 
mecánicas. Deformación (mm) y Esfuerzo 
al Límite Proporcional (MPa), una se 
corresponde con la otra, vale decir, es la 
deformación o deflexión (mm) de la viga 
como respuesta a una carga (MPa) límite 
que permite mantener sus propiedades 
elásticas. Si se aplica una carga (MPa) 
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Tabla 6. Análisis multivariado de componentes principales, variación explicada para la 
propiedad mecánica Módulo de Ruptura (MPa), considerando variables de manufactura 

y anatómicas en vigas laminadas (n=18 vigas). 


Aspectos de crecimiento y manufactura 

Valor propio 

Proporción 

Acumulada 

Número de Nudos 

2.2775 

0,456 

0,456 

Número de Uniones Finger Joint 

1,3886 

0.278 

0,733 

Tipo de Viga 

1,0546 

0,211 

0,944 

% Leño Juvenil 

0,2572 

0,051 

0,996 

Esfuerzo al Limite proporcional (MPa) 

0,0220 

0,004 

1,000 


Tabla 7. Análisis multivariado de componentes principales, variación explicada para la 
propiedad mecánica Módulo de Elasticidad (MPa), considerando variables de 
manufactura y anatómicas en vigas laminadas (n=18 vigas). 


Aspectos anatómicos y manufactura 

Valor propio 

Proporción 

Acumulada 

Número de Nudos 

2,4028 

0,481 

0,481 

Número de Uniones Finger Joint 

1 ,6001 

0,320 

0,801 

Tipo de Viga 

0,5915 

0,118 

0,919 

% Leño Juvenil 

0,3841 

0,077 

0,996 

Esfuerzo al Límite proporcional (MPa) 

0,0215 

0,004 

1,000 



Figura 3. Relación de variables de manufactura y anatómicas, con respecto a la 
propiedad mecánica Módulo de Ruptura (MPa), mediante el resultado gráfico 
de análisis multivariado de componentes principales. 
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mayor, la viga presentará una deforma- 
ción permanente. Para ambas propie- 
dades mecánicas la variable de mayor 
incidencia se corresponde con el Número 
de Nudos, representando este defecto 
discontinuidad o interrupción en los 
tejidos de madera que conforman una 
viga y que soportan una carga específica. 
Una viga con pocos nudos tendrá una 
mayor magnitud de deflexión frente a una 
carga específica, puesto que se trata de 
un elemento más homogéneo y menos 
rígido frente a la aplicación de cargas. 

Las variables de manufactura y 
anatómicas se relacionaron con la 
variable Módulo de Ruptura (mor) (MPa) 
mediante el análisis de componentes 
principales, del cual se destaca que 
45,6% de la variación se debe al Número 
de Uniones, 27,8% al Tipo de Viga, 
21,1% para Número de Nudos y apenas 
5,1% por causa de! %Leño Juvenil. En la 
figura 3 se observa la relación de varia- 
bles de manufactura y anatómicas, con 
respecto al Módulo de Ruptura (mor), 
donde se destaca que la propiedad 
mecánica resultó con un comporta- 
miento indirectamente proporcional con 
la variable de manufactura Número de 
Uniones, siendo esta variable la que 
explica 45,6% de la variación del Módulo 
de Ruptura (mor) (MPa). De la posición 
de vectores se deduce que en la medida 
que se incrementa el Número de 
Uniones, los valores de Módulo de 
Ruptura (mor) (MPa) se hacen menores. 
La posición de los vectores de las varia- 
bles anatómicas Número de Nudos y el 
%Leño Juvenil, con respecto a la posi- 
ción del vector de la variable mecánica 
en cuestión, expresan una mínima inci- 
dencia en su magnitud. 

El Módulo de Ruptura (MPa) expresa 
la carga necesaria para producir la 
ruptura de la viga por falla en flexión 
(deformación elástica, deformación 
permanente o plástica, ruptura). En la 


medida que el número de uniones fínger 
joint se incrementan en los elementos 
constitutivos de las vigas laminadas, la 
carga necesaria para generar su ruptura 
o falla por flexión en el Módulo de Ruptura 
(mor), es menor. Es decir, que en la 
medida que se minimicen el número de 
uniones finger joint entre la luz de carga, 
se incrementará la resistencia de la viga. 
La unión de finger joint es menos elástica 
que la madera, y es más rígida que la 
madera sin unión, es por ello que el incre- 
mento de número de uniones determina 
que el valor de la deformación elástica 
sea muy cercano al valor de ruptura. Una 
viga sin uniones tiende a presentar 
valores de deformación elástica notable- 
mente menores que su correspondiente 
valor de ruptura. 

La variación propiedad mecánica, 
Módulo de Elasticidad (MPa) en vigas 
laminadas es explicada en 48,1% por la 
variable Leño (madera) Juvenil, 32,0% 
por el Número de Nudos, 11,8% por el 
número de Uniones y apenas 7,7% por 
el leño juvenil (Tabla 7). En la figura 4 se 
observa la relación de variables de manu- 
factura y anatómicas, con respecto al 
Módulo de Elasticidad (MPa), donde se 
destaca que la propiedad mecánica en 
cuestión resultó con un comportamiento 
indirectamente proporcional con la 
variable anatómica %Leño Juvenil, 
siendo esta variable la que explica la 
mayor variación. 

Es decir, en la medida que se incre- 
menta el %Leño Juvenil en una viga lami- 
nada los valores de Módulo de Elasti- 
cidad (MPa) se hacen menores. La posi- 
ción de los vectores de las variables de 
manufactura Tipo de Viga y Número de 
uniones, con respecto a la posición del 
vector de la variable mecánica en cues- 
tión, expresan una mínima incidencia en 
su magnitud. 
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Figura 4. Relación de variables de manufactura y anatómicas, con respecto a la 
propiedad mecánica Módulo de Elasticidad (MOE) (MPa), mediante el resultado gráfico 
de análisis multivariado de componentes principales. 


El Módulo de Elasticidad (MPa) indica la 
capacidad que tiene un material para recu- 
perar su forma y tamaño original después de 
remover los esfuerzos, es decir, es una 
medida de resistencia del material. La 
variable anatómica que repercute de 
manera más contundente en la capacidad 
de una viga laminada de soportar cargas sin 
deformarse luego de ser retirada la carga se 
corresponde con el Leño Juvenil presente. 
La madera juvenil del Pinus caribaea var. 
hondurensis de la Orinoquía Venezolana se 
caracteriza por presentar menor densidad, 
menor número de anillos de crecimiento 
{Ninin, 2007), y menores longitudes de 
traqueidas que la madera adulta (Trejo, 
2006). Este conjunto de variables afecta la 
resistencia mecánica de la madera juvenil. 

CONCLUSIONES 

El presente estudio permite concluir que 
las propiedades mecánicas de las vigas 


laminadas encoladas de Pinus caribaea 
var. hondurensis de la Orinoquía Venezo- 
lana, se ven afectadas según los defectos 
evaluados, en et siguiente orden de 
importancia: 

1 . Número de Nudos afecta la capacidad de 
deformación o deflexión (mm) de la viga. 

2. Número de Uniones afecta los valores 
de Módulo de Ruptura (mor) (MPa). 

3. Porcentaje de Leño Juvenil afecta los 
valores del Módulo de Elasticidad (MPa) 

De acuerdo con los objetivos planteados, la 
determinación de la caracterización por 
niveles de importancia de los defectos en las 
vigas laminadas encoladas de Pinus cari- 
baea var. hondurensis, permiten reco- 
mendar a los fabricantes de este tipo de 
vigas que se debe tener presente la mejora 
en los procesos de fabricación, especial- 
mente en la técnica de finger joint, la cual 
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Figura 5. Apreciación visual y técnica de la cantidad y tamaño de los nudos en las 
láminas de madera que conforman una de las vigas laminadas realizadas con la técnica 
de fingerjoint Se puede apreciar el mal uso de esta técnica al quedar, en este caso de 
estudio, las uniones muy cercanas entre sí, lo cual ocasiona la falla y una menor 
resistencia del elemento laminado al momento del ensayo. 

(Foto: Wilver Contreras Miranda) 



Figura 6. Apreciación visual de una unión fingerjoint ubicada en la zona de tracción, de 
una de las vigas estudiadas, contribuyendo a ocasionar la falla del elemento laminado 
al momento del ensayo. (Foto: Wilver Contreras Miranda) 
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viene siendo mal empleada según los crite- 
rios de uso formulados, entre otros, por 
Freas y Seibo (1954), American Association 
of State Highway OfRcials (1973), American 
Institute of Timber Construction (1984), 
American Institute of Timber Construction 
(1985), Contreras y Owen de Contreras 
(1999) o Contreras et al. (2000). 

La mejora del uso de la técnica permi- 
tirá aumentar los valores en los esfuerzos de 
diseño y una mejor proyección de seguridad 
en el uso de las piezas estructurales en las 
edificaciones. Respecto a los porcentajes en 
la presencia de madera de leño juvenil, inhe- 
rente a la calidad anatómica de la madera 
de Pinus caribaea var. hondurensis, se reco- 
mienda no incluir piezas con importantes 
porcentajes de este defecto y de tal forma 
obtener una calidad de uso estructural de 
forma consistente. Para evitar de manera 
contundente la participación del leño juvenil, 
se recomienda para la obtención de piezas 
aserradas que conforman las vigas, la apli- 
cación del reaserrado con la conicidad, vale 
decir, que siga la conicidad de las trozas, de 
tal forma que al encajonar el leño juvenil en 
una cuña central se evitará que éste se 
mezcle con madera aserrada madura en 
una misma pieza (Ninin, 2007). 

En referencia a los nudos, es inherente 
la cantidad de éstos en el fuste del árbol de 
Pinus caribaea var. hondurensis de la Orino- 
quía Venezolana, se tiene una gran 
presencia y cercanía de los verticilos de 
nudos. Esto afecta notablemente la calidad 
estructural de esa madera. 

Se recomienda desarrollar nuevas 
investigaciones que permitan evidenciar y 
ponderar el efecto del porcentaje de leño 
juvenil en piezas que conforman las vigas 
laminadas encoladas de Pinus caribaea 
var. hondurensis, siendo para ello impres- 
cindible controlar o congelar el efecto del 
resto de defectos sobre las resistencias 
mecánicas, tales como nudos y número 
de uniones de finger joint. 
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Contenido fenólico y acción antioxidante 
de extractos de acículas de Pinus cooperi, 
P. durangensis, P. engelmanii y P. teocote 

Phenolic contení and antioxidant action in needle extracte 
from Pinus cooperi, P. durangensis, P. engeimanii 

and P teocote 


María de los Ángeles Sáenz-Esqueda'', Martha Rosales-Castro2, 
Nuria Elizabeth Rocha-Guzmán'i, José Alberto Gallegos-Infante'' 
y Rubén Francisco González-Laredo'' 

RESUMEN 


Durango es uno de los estados con mayor riqueza forestal, siendo el pino el género dominante 
en sus bosques, de ahí el interés por evaluar química y biológicamente extractos de acículas de Pinus 
cooperi (COO), P durangensis (DUR), P. engelmanii (ENG) y P teocote (TEO). Se practicaron extrac- 
ciones sucesivas con acetona y metanol acuosos de acículas de pino desengrasadas previamente 
con hexano, en las cuales DUR y ENG mostraron respectivamente los mejores rendimientos. Los 
valores más altos del contenido de fenoles totales fueron los extraibles acetónicos de TEO y los meta- 
nólicos de ENG. Se realizó una purificación cromatográfica a los extraibles acetónicos de TEO, gene- 
rando cinco fracciones; F3-TEO con 489,6 EC presentó el mayor contenido fenólico. Los extractos 
acetónicos y metanólicos de TEO y la fracción purificada F3-TEO presentaron la mayor concentración 
de taninos condensados. En el ensayo de DPPH, los extractos más relevantes según su concentra- 
ción efectiva media (CE'so) fueron los acetónicos de TEO (223,5 ppm) y las fracciones F3-TEO (145,0 
ppm) y F4-TEO (150,8 ppm). De manera similar, en el ensayo de ABTS, los extractos de TEO en 
acetona (27,7%) y en metanol (18,2%), asi como las fracciones F3-TEO (60,1%) y F4-TEO (55,8%) 
mostraron la supresión más elevada de la oxidación. En el ensayo de la oxidación de LDL, la mayor 
actividad inhibitoria fue mostrada por extractos acetónicos, particularmente por ENG (99,4%) a una 
concentración de 90 ppm. Los extractos acetónicos de TEO y de ENG presentan potencial como 
fuente de polifenoles bioactivos. 

PALABRAS CLAVE; 

ABTS, acículas. antioxidante, DPPH, LDL, polifenoles. 


ABSTRACT 

Coniferous trees are a vast resource in Durango, México forests and a potential nontimber 
source of these Chemicals, thus. our study for analyzing phenolic contení and biological activity in 
needle extracts from Pinus cooperi (COO), P. durangensis (DUR), P engelmanii (ENG) and P. teocote 
(TEO). Pine needles were dewaxed with hexane and extractad with aqueous acetona and methanol. 
DUR and ENG extracts showed the higher yields, while acetone TEO and methanol ENG extracts 
showed the highest catechin equivalent (CE) contents. A partial chromatographic purification on 
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acetone TEO extracta produced five 
fractions, the third one (F3-TEO) had the 
highest phenolic content (489,6 CE). In 
condensed tannins contenta, acetone 
and methanol TEO extracta, aa well frac- 
tion F3-TEO were the beat (24,3, 15.7 
and 44,4 mg/g extract, reapectiveiy). In 
DPPH activity, the more effective CE 50 
were obtained from TEO acetone extract 
(223,5 ppm) and fractiona F3-TEO 
(145,0 ppm) and F4-TEO (150,8 ppm). In 
ABTS aaaay, TEO acetone (27,7%) and 
methanolic (18,2%) extracta, and frac- 
tiona F3-TEO (60,1%) and F4-TEO 
(55,8%) were the beat. In LDL aaaay the 
higheat inhibitory activity waa ahown by 
acetone extracta, particularly ENG 
(99,4%) at 90 ppm. Therefore acetone 
extracta from TEO and ENG may have 
utilization aa a potential aource of bioac- 
tive phenolica. 

KEYWORDS; 

ABTS, Antioxidanta, DPPH, LDL, 
needlea, polyphenola. 


INTRODUCCIÓN 

El particular hábitat de la sierra mexicana 
ha favorecido la mayor variedad de espe- 
cies de Pinus en el mundo. El bosque de 
pino se caracteriza por ser una comu- 
nidad siempre verde asociada con enci- 
nares y otras especies, y es el de mayor 
importancia económica para la industria 
forestal de México, sustentando activi- 
dades forestales como aserrío, resina- 
ción, producción de pulpa celulósica, 
postaría y recolección de semilla. Sin em- 
bargo, el follaje conifero está por recibir 
aún un adecuado aprovechamiento como 
recurso natural no maderable. 

La vegetación forestal está domi- 
nada por pinos con alturas promedio de 
15 m a 30 m. Estos pinares tienen un 
estrato inferior relativamente pobre en 
arbustos, pero con abundantes gramí- 
neas, esta condición los expone a 
frecuentes incendios y a tala inmoderada. 
Las especies más comunes en México 
son pino (P. hartwegií), ocote blanco (P 
montezumae), pino lacio (P pseudos- 


trobus), pino (P rudis), pino escobetón (P. 
michoacana), pino chino (P teocote), 
ocote trompillo (P oocarpa), pino ayaca- 
huite (P ayacahuite), pino (P pringleí), R 
duranguensis, P. chihuahuana, P. engel- 
manii, P. lawsonii, R oaxacana, entre 
otros (Conafor, 2000-2006). 

Durango, México es uno de los 
principales productores de madera de 
pino, contando con aproximadamente 
veinte especies, de las cuales diez son 
las más aprovechadas por su abun- 
dancia y calidad de madera, algunas de 
ellas son Pinus durangensis, P. arizo- 
nica, P. engelmannii y P cooperi (García 
y González, 2003). Otras especies 
representativas son P. leiophylla, P. 
ayacahuite y P. chihuahuana. Del apro- 
vechamiento de la madera se generan 
altos volúmenes de corteza que repre- 
sentan de 10% a 15% del peso total del 
árbol (Kofujita et al., 1999). Adicional- 
mente se generan volúmenes no cuanti- 
ficados de follaje de alta riqueza fitoquí- 
mica. Los extraíbles o metabolitos 
secundarios presentes en la madera, 
corteza y follaje de las pináceas, son 
múltiples compuestos orgánicos de 
naturaleza diversa que incluyen grasas, 
ceras, ácidos resínicos, terpenoides, 
esteroides y polifenoles (flavonoides y 
taninos). Todos ellos son producidos por 
el árbol como protección o defensa ante 
ataques microbianos o de insectos 
(Perry, 1991). A los metabolitos secun- 
darios se les consideraba sustancias de 
desecho o errores de las rutas metabó- 
licas primarias de las plantas. Sin 
embargo, su gran abundancia, su 
variedad estructural y el descubrimiento 
de rutas biosintéticas y de enzimas 
específicas para su síntesis han 
cambiado este enfoque. Un caso parti- 
cular son los polifenoles o taninos vege- 
tales que poseen al menos un anillo 
aromático con grupos hidroxilo en su 
estructura química. Los polifenoles son 
micronutrientes abundantes en nuestra 
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dieta y existen evidencias de su rol en la 
prevención de enfermedades degenera- 
tivas como el cáncer y enfermedades 
cardiovasculares y neurodegenerativas 
(Markham, 1989). Los polifenoles tienen 
diversa función, tamaño y arreglo 
estructural y se pueden clasificar en 
ácidos fenólicos, flavonoides, estilbenos 
y lignanos. Los flavonoides comparten 
una estructura común de dos anillos 
aromáticos unidos por una cadena 
propánica (Bors et al., 1990). En función 
de sus características estructurales se 
pueden clasificar en chalconas, 
flavanos, flavanoles, flavonas, antocia- 
nidinas e isoflavonoides. Entre otras, se 
les atribuyen actividades antifúngica y 
bactericida, e importante capacidad 
para fijar metales como hierro y cobre, 
además de conferir color que promueve 
fenómenos de polinización (Fórmica y 
Regelson, 1995). Algunas patentes 
norteamericanas registran el uso de 
extractos de hojas de pino como ingre- 
dientes para combatir la alopecia 
(Uchida et al., 1996), como ingrediente 
para preparar tofu (Im, 2000), bebidas 
tonificantes (Kang, 1987) y nutracéutico 
cardiprotector (Ji, 2001). 

OBJETIVO 

A través de la evaluación del contenido 
fenólico y de la acción bioquímica de los 
extractos de acículas de Pinus cooperi 
(coo), P. durangensis (dur), P engelmanii 
(eng) y P. teocote (teo), determinar el 
potencial antioxidante del follaje de estas 
especies, que fundamente su eventual 
aprovechamiento. 

METODOLOGÍA 

Material de estudio. Las acículas se 
colectaron en el ejido Chavarría Viejo, del 
municipio de Pueblo Nuevo, Durango. 
Las muestras fueron identificadas y depo- 
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sitadas en el Herbario del ciidir-ipn 
Unidad Durango. Las acículas se secaron 
a la sombra, se cortaron en segmentos 
pequeños y se molieron. 

Extracciones sucesivas. El material 
previamente desengrasado con hexano 
se extrajo con acetona acuosa (70%) y 
se separó por filtración; el extracto se 
concentró en rotaevaporador y se liofi- 
lizó. El material remanente fue extraído 
con metanol acuoso (50%), filtrado y el 
extracto, rotaevaporado y liofilizado. 

Evaluación de fenoles totales. Se 

prepararon muestras de extractos crudos 
y fracciones purificadas a 100 ppm en 
agua libre de CO 2 . A una muestra (125 
|aL) se le adicionó agua (500 \xL) y reac- 
tivo Folin Ciocalteu 2N (125 ¡.iL), se agitó 
y dejó reposar 6 minutos, luego se agregó 
Na 2 C 03 7% (1,25 mL) y agua libre de 
CO 2 (1 mL). Los tubos se agitaron nueva- 
mente y dejaron reposar (90 minutos), 
protegidos de la luz (Waterman y Mole, 
1994). Después se midió la absorbancia a 
760 nm en un espectrofotómetro Cary50. 
Se construyó una curva de calibración de 
catequina como estándar (40 pg/mL -120 
pg/mL). Los resultados se reportan como 
miligramos equivalentes de catequina por 
gramo de extracto seco. 

Evaluación de taninos condensados. 

Los extractos crudos y fracciones purifi- 
cadas se prepararon a 1000 ppm. A 50 pL 
de muestra se le agregaron 3 mL de vaini- 
llina (4% en metanol), 1,5 mL de HCI 
(concentrado), se agitó y dejó reposar 15 
minutos en la oscuridad. Después se 
midió la absorbancia a 500 nm contra un 
blanco de metanol. Al mismo tiempo se 
construyó una curva de calibración con 
catequina en un rango de 100 ppm a 600 
ppm (Heimier etal., 2005). Los resultados 
se reportan como miligramos equiva- 
lentes de catequina por 100 mg de 
extracto. 
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Evaluación de capacidad antioxidante 

Método del dpph (2,2-difenil-1-picrílhidra- 
cilo). Los extractos acetónicos se 
evaluaron a 200, 250, 300 y 350 ppm, los 
de metanol a 500, 750, 1000 y 1200 ppm 
y las fracciones purificadas a 1 00 y 2000 
ppm. Se colocaron 2900 fiL del reactivo 
de DPPH y 100 fxL de muestra, se agitó y 
midió la absorbancia a 515 nm a tiempo 
cero y a 30 minutos (Brand-Williams et 
al., 1995). Se calculó la concentración 
efectiva media (CE 50 ) de cada muestra 
como indicador de la concentración 
requerida para estabilizar 50% del radical 
presente. Los resultados se reportan 
como CE 50 en ppm. 

Método del abts (ácido 2,2'-azinobis-3- 
etilbenzotiazolin-6-sulfónico). Las mues- 
tras se probaron a 100 ppm en etanol 
absoluto. El radical se diluyó a una absor- 
bancia de (0,70±0,1) nm a 754 nm. Se 
utilizaron 40 ¡.iL de muestra y 960 ¡.iL de 
radical y se leyó a 754 nm; la absorbancia 
se monitoreó a 1 y 6 minutos después de 
adicionar el radical (Re et ai, 1999). Los 
resultados obtenidos se reportan como 
capacidad atrapadora de radicales libres 
(%rsa). 

Inhibición de la oxidación de Lipoprote- 
ínas de Baja Densidad (ldl). Las ldl se 
precipitaron a partir de plasma de volun- 
tarios sanos. Las muestras se probaron a 
30, 60 y 90 ppm. En tubos falcon de 15 
mL, se colocaron 400 ¡.iL de buffer de 
fosfatos, 100 |.iL de sulfato de cobre, 100 
\xL de muestra y 400 ¡.iL de ldl. Se 
agitaron e incubaron a 37 °C por 180 
minutos. Luego se agregó 1 mL de ácido 
tricloroacético (5%) y 1 mL de ácido 
tiobarbitúrico (0,37%), se agitó por 20 
minutos en baño María a punto de ebulli- 
ción, después se enfrió en hielo por 10 
minutos y se agregaron 2 mL de butanol, 
se mezcló bien y centrifugó a 3500 rpm 
(10 min). Se extrajo el sobrenadante 
(fracción soluble en butanol) y se leyó a 


532 nm contra un blanco de butanol (Re 
et al., 1999). Se utilizó catequina como 
referencia. Los resultados se reportan 
como porcentaje de inhibición de oxida- 
ción de LDL. 

Purificación en columna. Una parte del 
extracto crudo de acetona se purificó en 
una columna con Fractogel hw- 40 . Se 
eluyó primero con metanol, separando 
dos fracciones de 100 mL (F1, F2), poste- 
riormente se eluyó con acetona 70%, y se 
obtuvieron tres fracciones más de 150 mL 
(F3, F4 y F5). 

Identificación por hplc. Se empleó un 
sistema hplc (Waters 2489) con columna 
C 18 (250 mm x 4,6 mm), usando flujo de 
1 mL/min de fase móvil A, agua acidifi- 
cada (pH 3,3) con ácido fosfórico, y fase 
móvil B, acetonitrilo (Sakakibara et al., 
2003). Se inició 5 minutos con 100% de A, 
luego 30 minutos con 25% de B, 10 
minutos a 50% y 20 minutos con 100% de 
A. Se emplearon los estándares, cate- 
quina, catequina galato, galocatequina, 
galocatequina galato, ácido caféico, ácido 
gálico y ácido elágico. 

Análisis estadístico. Los datos obtenidos se 
valoraron mediante análisis de varianza de 
una sola vía (a= 0,05). Se realizó una 
prueba de comparación de medias por el 
método de Tukey (a= 0,05). En todos los 
casos los experimentos fueron realizados al 
menos por duplicado. 

RESULTADOS 

Los rendimientos obtenidos en las extrac- 
ciones sucesivas de aciculas de pino con 
acetona 70% y metanol 50% se muestran 
en la tabla 1 . Salvo por el extracto coo, 
los extractos acetónicos presentaron un 
rendimiento similar cercano a 10%, mien- 
tras que ENG sorprende con un rendi- 
miento aún mayor en el extracto meta- 
nólico. 
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Tabla 1 . Rendimiento obtenido en las extracciones sucesivas con acetona (70%) 
y metanol (50%) de acículas de Pinus teocote, P. engelmanii, 

P cooperi y P durangensis. 


ESPECIE 

Acetona (70 %) 

% RENDIMIENTO 
Metanol (50 %) 

Acumulado 

TEO (P. teocote) 

9,41 ± 0,81a 

1,12 ±0.10^ 

10,53 

ENG (P engelmanii) 

9,22 ± 1,06a 

9,63 ± 0,85b 

18,85 

COO (P cooperi) 

4,90 ± 0,05b 

2,52 ± 0 , 33 c 

7,42 

DUR (P. durangensis) 

10,08 ± 0,56a 

6,08 ± 0,50d 

16,16 


Superíndices diferentes en cada columna significan diferencia estadística. Cada valor representa el 
promedio de dos experimentos realizados por duplicado (Tukey, p< 0,05) 


La evaluación de fenoles totales se 
muestra en la tabla 2. Los extractos 
acetónicos oscilaron entre 169,6 y 399,2 
equivalentes de catequina (ec), mientras 
que los extractos metanólicos resultaron 
considerablemente más bajos (p<0,05) 
de 49,0 a 140,6 ec. El contenido de 
fenoles totales en las fracciones purifi- 
cadas a partir del extracto acetónico de 
TEO varió entre 163,5 y 489,6 ec (Tabla 3), 
en un rango similar al que presentan los 
extractos acetónicos. 

El contenido de taninos condensados 
en los extractos acetónicos se presentó de 
18,4 a 24,3 mg/g y en los metanólicos de 
13,2 a 15,7) de extracto seco (Tabla 2), 
siendo diferentes unos de otros (p<0,05). 
Los taninos condensados en los extractos 
acetónicos fueron mayores, mostrando la 
misma tendencia que su contenido de 
fenoles totales. De manera similar, en las 
fracciones purificadas del extracto acetó- 
nico de TEO se encontró un contenido de 
taninos condensados disperso de 6,8 a 44,4 
mg/g de extracto seco, siendo F3-teo y fs- 
TEO las fracciones con la mayor concentra- 
ción (Tabla 3). 

En la evaluación de capacidad antio- 
xidante por el método del dpph, la CE50 
más eficiente que se obtuvo de los 


extractos crudos correspondió al extracto 
acetónico teo con un valor de 223,5 ppm 
(Tabla 2). Entre los extractos metanólicos, 
esta misma especie fue la única que 
mostró actividad antioxidante suficiente 
por encima del umbral de 50%, aunque a 
concentraciones muy altas (784,1 ppm), 
ya que para las condiciones del análisis, 
los otros tres extractos mostraron una 
actividad máxima del orden de 30%, que 
es insuficiente para estimar su CE50. 

Entre las fracciones purificadas, 
destacan la F3-teo y la f4-teo con CE50 de 
145,0 ppm y 150.8 ppm, respectivamente, 
como las más eficientes de las muestras 
analizadas (Tabla 3). Sin embargo, estos 
valores están aún alejados del estándar 
catequina con una CE50 de 90,5 ppm. Se 
encontraron diferencias significativas 
(p<0,05) entre los extractos acetónicos 
excepto para coo y eng. De igual manera 
se analizaron los resultados de las frac- 
ciones purificadas, mismas que 
mostraron diferencias significativas 
(p<0,05) entre los extractos excepto entre 
F3-TEO y F4-TEO. 

Los resultados de la actividad antioxi- 
dante de los extractos acetónicos ante el 
radical libre abts, expresados como 
porcentaje de actividad antioxidante, se 
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Tabla 2. Cuantificación de fenoles totales, taninos condensados y actividad antioxidante 
de los extractos crudos de aciculas de Pinus teocote (teo), R engelmanii (eng), 

P. cooperi (coo) y P durangensis (dur) con acetona (70%)* y metanol (50%)°, 

comparados con catequina (cat). 


ESPECIE 

FENOLES TOTALES 
Equivalentes 
de catequina (EC) 

TANINOS CONDENSADOS 
mg/g ext seco 

DPPH 

ppm 

ABTS 

% 

TEO* 

399,20 ± 2,42a 

24,31 ± 1,84a 

223,54 ± 2,04a 

27,70 ± 1,20^ 

ENG* 

312,39 ± 4,66b 

23,10 ± 1,10a 

278,66 ± 4,70b 

23,36 ± 2,75b 

COO* 

169,59 ± 9,71c 

20,22 ± 1,03b 

295,60 ± 3,43b 

27,09 ± 0,24® 

DUR* 

217,81 ± 2,52d 

18,39 ± 1,36b 

403,13 ± 1,45c 

18,77 ± 1,35c 

TEO° 

94,99 ± 7,48e 

15,71 ± 0,92c 

784,12 ± 7,60d 

18,18 ± 2,86c 

ENG° 

140,61 ± 9,45f 

13,72 ± 0,32d 


10,27 + 2,25d 

COO“ 

48,99 ± 3,429 

13,17 ± 0,62d 


10,42±2,42d 

DUR° 

79,91 ± 6,91e 

13,16 ± 0,45d 


13,38 ± 1,97d 

CAT 



90,47 ± 2,20® 

66,83 ± 6,57® 


Superíndices diferentes en cada columna significan diferencia estadística. Cada valor representa el 
promedio de dos experimentos realizados por duplicado (Tukey, p< 0,05) 


Tabla 3. Cuantificación de fenoles totales, taninos condensados y actividad antioxidante 
de extracto acetónico de aciculas de Pinus teocote (teo), de sus fracciones purificadas 

(F1, F2, F3, F4, F5) y catequina (cat). 


ESPECIE 

FENOLES TOTALES 
Equivalentes 
de catequina (EC) 

TANINOS CONDENSADOS 
mg/g ext seco 

DPPH 

CEso 

ppm 

ABTS 

% 

TEO 

399,20 ± 2,42a 

24,31 ± 1,85a 

223,54 ± 2,04a 

27,70 ± 1,20® 

F1-TEO 

163,51 ± 2,48b 

6,85 ± 0,54b 

1930,05 ± 12,12b 

33,72 ± 1,68® 

F2-TEO 

379,30 ± 8,68c 

35,98 ± 5,94c 

297,47 ± 2,40c 

32,51 ± 3,52® 

F3-TEO 

489,65 ± 3,23d 

44,37 ±4, 61d 

145,03 ± 1,21d 

64,33 + 8,00b 

F4-TEO 

476,14 ± 0,99d 

28,50 ± 2,57a 

150,78 ± 0,35d 

58,31 ± 4,44b 

F5-TEO 

438,07 ± 9,92® 

42,16 ± 3,44cd 

268,29 ± 5,86® 

36,03 ± 3,53® 

CAT 



90,47 ± 2,20f 

66,83 ± 6,57b 


Superíndices diferentes en cada columna significan diferencia estadística. Cada valor representa el 
promedio de dos experimentos realizados por duplicado (Tukey, p< 0,05) 


encontraron en el rango de 18,8% a 
27,7%, mientras que los metanólicos de 
10,3% a 18,2% (Tabla 2) muy alejados del 
estándar catequina (66,83%). Sin 
embargo las fracciones purificadas de 
TEO oscilaron entre 32,51% a 64,33%, 
donde F3 y F4 presentaron una actividad 
comparable a dicho estándar (Tabla 3). 


Los porcentajes de inhibición de la 
oxidación de las ldl mostraron tendencia 
similar a las determinaciones anteriores, 
ya que los extractos acetónicos 
mostraron los mejores resultados (Tabla 
4). Sin embargo, los extractos de eng 
presentaron el mayor porcentaje de inhi- 
bición (99,4%), superando incluso al 
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estándar catequina (93,0%)(p<0,05). Por 
el contrario, los extractos metanólicos 
(Tabla 5) de teo y eng, aunque pueden 
considerarse buenos con una inhibición 
de casi 40%, dado su carácter residual, 
comparados con los extractos acetónicos 
son apenas la mitad. 

En el análisis cromatográfico por 
HPLC de las 13 muestras analizadas, 11 
resultaron positivas cualitativamente para 
el compuesto catequina-galato. Sin 
embargo, sólo en cinco de ellas se les 
encontró de manera significativa, como 
fue el caso de los extractos crudos acetó- 
nicos de TEO, ENG y DUR, así como en las 
primeras fracciones purificadas F1 y F2, 


del extracto acetónico de acículas de teo, 
aunque la F3 fue la que mostró mayor 
respuesta antioxidante (Figura 1). 

DISCUSIÓN 

Como era previsible el rendimiento de los 
extractos acetónicos es mayor que los 
metanólicos, excepto por los extractos de 
ENG que muestran contenidos similares 
en ambos casos (p<0,05). El rendimiento 
de las extracciones acetónicas de las 
cuatro especies es similar excepto en los 
extractos de coo que corresponden apro- 
ximadamente a la mitad del resto de los 
pinos (p<0,05). En los extractos metanó- 


Tabla 4. Inhibición de la oxidación de las ldl por extractos con acetona (70%) de 
acículas de Pinus teocote (teo), P. engelmanii (eng), P cooperi (coo), 

P durangensis (dur), comparados con catequina (cat). 


ESPECIE 

30 ppm 

60 ppm 

90 ppm 

TEO 

22,15 ± 5,01^ 

35,46 ± 1,85^ 

50,98 ± 3,31^ 

ENG 

84,75 ± 0,71b 

91,30 ± 0,65b 

99,37 ± 0,11b 

COO 

16,19 ± 0,03^ 

60,65 ± 0,12c 

73,62 ± 0,70c 

DUR 

23,23 ± 0,68a 

30,87 ± 2,59d 

71,72 ± 1,71c 

CAT 

84,53 ± 0,71b 

91,27 ± 0,65b 

99,37 ± 0,11b 

Superíndices diferentes en cada columna significan diferencia estadística, 
promedio de dos experimentos realizados por duplicado (Tukey, p< 0,05) 

Cada valor representa el 


Tabla 5. Porcentaje de inhibición de la oxidación de ldl, por extractos secuenciados con 


metanol (50%) de acículas de Pinus teocote (teo), P engelmanii (eng), 
P cooperi (coo), P durangensis (dur), comparados con catequina (cat). 

ESPECIE 30 ppm 

60 ppm 

90 ppm 

TEO 

10,34 ± 1,86^ 

18,12 ± 0,53® 

39,29 ± 0,87® 

ENG 

10,55 ± 0,25b 

14,94 ± 1,24b 

38,63 ± 0,88b 

COO 

9,04 ± 1 ,70® 

17,37 ± 0,30a 

19,18 ± 0,17c 

DUR 

4,49 ± 0,77a 

8,26 ± 0,98c 

21,45 ± 0,83c 

CAT 

78,97 ± 0,23c 

89,36 ±0,1 2d 

96,96 ± 0,67d 


Superíndices diferentes en cada columna significan diferencia estadística. Cada valor representa el 
promedio de dos experimentos realizados por duplicado (Tukey, p< 0,05) 
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Figura 1. Perfil cromatográfico (hplc) de extractos acetónicos de acículas 
de Pinus teocote y Pinus engelmanii, así como fracción purificada F3 de P teocote, 

conteniendo catequina galato. 
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licos, ENG presentó el máximo valor 
acumulado cercano a 19%. Cabe señalar 
que la extracción de compuestos fenó- 
licos de material vegetal está influenciada 
por la naturaleza química de éste, por el 
método de extracción, el tamaño de partí- 
cula, el tiempo y las condiciones de alma- 
cenamiento (Prior y Cao, 1999). 

Se puede observar la mayor selecti- 
vidad por polifenoles del solvente acetona 
acuosa al comparar los valores obtenidos 
con metanol acuoso, aunque se debe 
recordar que son extracciones sucesivas. 
El metanol como segundo solvente sola- 
mente extrae los componentes fenólicos 
residuales que no fueron extraídos por la 
acetona (Heimier et al., 2007). Esta dife- 
rencia en el solvente de extracción deter- 
mina el rendimiento y calidad del conte- 
nido fenólico de ambos extractos (Stratil 
et al., 2006; Ninfali etal., 2005). 

La determinación de taninos 
condensados es un indicador del conte- 
nido de proantocianidinas presentes en 
los extractos crudos tanto acetónicos 
como metanólicos. Este es un grupo 
muy importante de flavonoides que 
brindan en su mayoría pigmentación a 
las plantas que los contienen. Al analizar 
individualmente los resultados de los 
extractos crudos de cada especie de 
pino, se encontró que teo, tanto en su 
extracto acetónico como en el metanó- 
lico, resultó con el mayor contenido de 
taninos condensados (p<0,05). En las 
fracciones purificadas a partir de teo se 
observó una mayor concentración de 
taninos en las fracciones F3-teo y f5-teo, 
lo que sugiere al menos dos fracciones 
polifenólicas oligoméricas de diferente 
grado de polimerización. Se encontraron 
diferencias significativas (p<0,05) entre 
los extractos acetónicos de teo y eng 
respecto a coo y dur, y los extractos 
metanólicos de teo y eng, así como en 
las fracciones f2-teo y f4-teo. 
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De acuerdo con el contenido de 
taninos condensados presente en 
extractos crudos y fracciones purificadas, 
se observa una relación entre el conte- 
nido de taninos condensados y la capa- 
cidad antioxidante de los mismos repre- 
sentada por su CE 50 del ensayo dpph. 
Esto indica que la actividad antioxidante 
de los extractos está en función directa 
con su contenido de taninos, y éstos 
juegan un papel muy importante en el 
metabolismo del género y sobre todo de 
la especie Pinus teocote (Katsube et al., 
2004). 

Al analizar la posible correlación entre 
el contenido fenólico y la actividad antioxi- 
dante de los extractos y considerando que 
los métodos utilizados para el análisis de 
ambos se fundamentan en el mismo meca- 
nismo de propensión a donar hidrógenos, 
se puede tener una lectura predictiva de 
que a mayor contenido fenólico en el 
extracto, mayor será su capacidad antioxi- 
dante (Katsube et al., 2004). Sin embargo, 
esta correlación no se presenta en todos los 
casos. La correlación lineal obtenida entre 
el blanqueamiento del dpph y el contenido 
fenólico de las fracciones purificadas fue de 
r 2 = 0,6836. 

El ensayo del radical libre abts se 
recomienda para medir la actividad antio- 
xidante de gran diversidad de sustancias, 
tanto en medios acuosos como lipidióos, 
sin embargo, en cada caso se encuentran 
diferentes radicales, por ejemplo, en 
medio lipídico se tienen radicales pero- 
xilos (Rice-Evans et al., 1996). Esta 
técnica se utiliza como referente para 
confirmar los resultados obtenidos de la 
prueba de dpph (Okonogi et al., 2007), ya 
que ambas se fundamentan en el mismo 
mecanismo de acción; la donación de 
electrones. Por ello se analizó una corre- 
lación lineal entre las dos técnicas, dando 
como resultado una r 2 de 0,684, lo cual 
nos indica que sí hay cierta relación entre 
las mismas. 
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En la inhibición de la oxidación de la 
LDL a 30 ppm de muestra, de acuerdo con 
el solvente de extracción, no se obser- 
varon diferencias significativas para las 
muestras acetónicas, pero sí para las 
metanólicas. A 60 ppm se encontraron 
diferencias significativas entre los cuatro 
extractos acetónicos y entre los metanó- 
licos sólo entre dur y eng. Finalmente, a 
90 ppm se encontraron diferencias signifi- 
cativas entre los extraíbles acetónicos y 
metanólicos de teo y eng (p<0,05). De 
acuerdo con lo anterior, se puede inferir 
que el extracto acetónico de eng es un 
buen inhibidor de la oxidación de ldl, sin 
embargo, se debe enfatizar su papel 
como inhibidor y no como agente 
previsor. Existen evidencias de que la 
formación de la placa aterosclerótica se 
da en la pared arterial, aislada de la 
acción de los antioxidantes plasmáticos, 
pero aún no hay sustento documentado 
de los beneficios aportados por el 
consumo de antioxidantes en la preven- 
ción o reducción del riesgo de enferme- 
dades cardiovasculares. Por otra parte, 
se ha demostrado que algunos antioxi- 
dantes con capacidad estabilizadora de 
radicales libres como fenoles simples y 
polifenoles pueden presentar efectos 
prooxidantes bajo determinadas condi- 
ciones (Pérez-Guerra, 2007). Sin 
embargo, el mecanismo por el cual se 
lleva a cabo la oxidación de ldl, es dife- 
rente al de las determinaciones anteriores 
y su actividad no debe subestimarse 
(Katsube et al., 2004). El mecanismo que 
se sigue es la quelación de metales 
óxido-reductores, lo que muestra que los 
antioxidantes actúan de manera diferente 
frente a radicales libres diferentes, 
presentando en cada caso un mecanismo 
de acción específico (Pérez-Guerra, 
2007). Así, resulta que los extractos de 
aciculas de teo mostraron la mejor acti- 
vidad antioxidante por atrapamiento de 
radicales libres, mientras que las de eng, 
la mejor actividad antioxidante por acción 
contra metales. 


La presencia dominante del compuesto 
catequina galato en extractos y fracciones 
purificadas activas sugiere una acción 
quelante y antioxidante por atrapamiento de 
radicales libres. Lo anterior confirma el 
carácter monomérico de las fracciones 
iniciales (F1 y F2). Sin embargo, la mayor 
actividad, contenido fenólico y de taninos 
condensados se presentaron en las frac- 
ciones subsecuentes (F3, F4, F5), resaltadas 
por el hecho que se trata de fracciones fenó- 
licas oligoméricas y poliméricas bioactivas 
que podrían ser evaluadas en ensayos de 
mutagenicidad (Kwak et al., 2006). 

CONCLUSIONES 

Las aciculas de pinos representan una 
fuente potencial de polifenoles aún por 
explotar. Los extraíbles polares presentes 
pueden representar hasta 19% en peso, 
como en el caso de R engelmannii. Estos 
extractos constituyen un significativo 
contenido fenólico rico en taninos 
condensados que se correlaciona con la 
actividad antioxidante mostrada contra 
los radicales libres dpph y abts, como 
resulta en el caso particular de los 
extractos de R teocote y sus fracciones 
purificadas. Por otra parte, se muestra 
también acción contra metales, inhi- 
biendo la oxidación de sistemas lipídicos 
como es el caso de los extractos de R. 
engelmanii. 

Catequina-galato se encontró de 
manera significativa en los extractos acetó- 
nicos de Rinus durangensis, R engelmanii y 
R. teocote al igual que en las fracciones puri- 
ficadas iniciales (F1 y F2) de este último. En 
estas fracciones fenólicas monoméricas, 
catequina-galato puede ser el componente 
fenólico más significativo que explique su 
acción estabilizadora sobre radicales libres 
y metales óxido-reductores. Además, esta 
capacidad quedó mas evidenciada en las 
fracciones fenólicas oligoméricas (F4 y F5). 
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La actividad antioxidante que 
presentó la especie R teocote fue la más 
relevante de las muestras evaluadas y 
está relacionada con el contenido de toni- 
nos condensados de la misma. Comple- 
mentariamente, la actividad mostrada por 
P. engelmannii en modelos lipidíeos 
resulta sobresaliente y puede relacio- 
narse con su alto contenido fenólico. 
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lARTlCULO DE INVESTIGACION! 


Variación estacional de compuestos fenólicos 
foliares en Quercus sideroxyla 
en diferentes tipos de suelo 

Seasonal varíation in the phenolic content of oak leaves 
(Quercus sideroxyla) in different soil texturas 

Jorge Armando Arámbula-Salazar'', B.l. Ibarra-Salinas2, 
Rubén Francisco González-Laredo'’-^, o.D. Muñoz-Galindo'' 

y H. Hernández-Vela2 

RESUMEN 


El propósito de este trabajo fue determinar la variación estacional y el efecto de la textura del 
suelo en el contenido fenólico (medido por el método de Folin-Ciocalteu) en hojas de encino colorado 
(Quercus sideroxyla). como parámetro de deterioro del género. Se muestreó el follaje de seis árboles, 
en ocho sitios diferentes, donde se determinaron los componentes edáficos en tres tipos de textura: 
franco arcillosa (12,5%), limo arcillosa (75%) y arcillo arenosa (12,5%), mismos que permanecieron 
inalterados durante el estudio. Se encontró para todos los suelos una disminución del contenido fenó- 
lico foliar en primavera y un aumento de éste en verano. Se determinó una correlación alta y positiva 
entre el contenido fenólico y la conductividad eléctrica en suelos de textura arcillo arenosa y también 
una correlación alta y positiva entre la retención de humedad del suelo y el contenido fenólico foliar, 
la cual fue independiente de la textura de los suelos. 

PALABRAS CLAVE: 

Encino, fenoles, foliar, Quercus sideroxyla, suelo. 


ABSTRACT 


The seasonal varíation and the soil texture effect were evaluated in the phenolic contente (Folin- 
Ciocalteu) of Quercus sideroxyla (red oak) leaves as a parameter of the genus detehoration. Leaves 
of six trees from eight sites were sampied and their edaphic characteristics determinad as soil texture: 
clayed (12,5%), lime-clayed (75%) and clay-sandy (12,5%). Soils and sites were unaltered during the 
four season collections. It was found in spring a notorious decrease in phenolic content of Quercus 
leaves followed by a substantial increase during summer. A high positiva correlation was found 
between phenolic content of leaves and electric conductivity of clay-sandy soils, and aiso between soil 
moisture retention and phenolic content, independently of the soil texture. 

KEYWORDS: 

Oak, phenolic, leaves, Quercus sideroxyla, soil. 
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INTRODUCCIÓN 

Los encinos, por su abundancia y distri- 
bución constituyen, después de los pinos, 
el segundo grupo de importancia entre las 
especies maderables de clima templado- 
frío de México (Conafor, 2007). Los 
encinos son árboles que pertenecen a 
uno de los géneros más importantes 
mundialmente: Quercus. Este género es 
el más grande de la familia Fagaceae 
(Rodríguez y Romero, 2007). Se conocen 
dos centros de diversidad para Quercus 
(encinos o robles), el primero al sureste 
de Asia con alrededor de 125 especies y 
el segundo en México, con una lista preli- 
minar de 161 especies para este género 
(Valencia-A., 2004). El género Quercus 
se encuentra distribuido en toda la Repú- 
blica Mexicana, excepto en los estados 
de Yucatán, Quintana Roo y Campeche, 
normalmente en las regiones templadas y 
frías o en lugares montañosos (Rze- 
dowski, 2006). Los bosques de encino 
han sido mal explotados maderablemente 
en nuestro país, el estado de Durango no 
es la excepción, aunque aproximada- 
mente 44% de su superficie tiene hábitat 
adecuado para el desarrollo del género 
(Bacon-Sobbe, 1997; Bonfil, 1998). El 
aprovechamiento del encino ha sido redu- 
cido debido al desconocimiento tanto 
taxonómico como tecnológico de las 
especies existentes (Zamora y Torres, 
2002), lo que ha propiciado la destrucción 
de grandes áreas forestales, así como el 
cambio de uso de extensas superficies 
que ha traído como consecuencia el dete- 
rioro ambiental (Martínez-Ruíz, 2005). 
Las plantaciones comerciales aún no 
tienen eco en Durango; los ejidos y comu- 
nidades se atienen a la reproducción de 
la naturaleza, pero son las plantaciones 
las que podrían redundaren un desarrollo 
sustentadle del bosque (Meléndez- 
Moreno, 2008). 

De ahí que sea importante conocer 
el tipo de suelo y las condiciones que 


favorecen la producción forestal. Las 
plantas carecen de un sistema inmuno- 
lógico, pero sintetizan compuestos que 
han sido clasificados como metabolitos 
secundarios (Hadacek, 2002). Estos 
compuestos han sido considerados 
como productos de desintoxicación, 
sobreexpresiones del metabolismo 
primario, productos de degradación y 
productos de almacenamiento (Hada- 
cek, 2002). Los metabolitos secundarios 
son clasificados de acuerdo con su 
origen biosintético en tres grandes 
grupos: polifenoles, terpenos y com- 
puestos nitrogenados (Taiz y Zeiger, 
1991). Los polifenoles son los metabo- 
litos secundarios más ampliamente 
distribuidos dentro del reino vegetal, son 
sustancias que poseen al menos un 
anillo aromático con un radical hidroxilo 
sustituyante en su estructura química 
(Markham, 1982). Además juegan un 
papel importante en la defensa contra el 
herbivorismo y suelen incrementarse en 
situaciones de estrés (Salminen et al., 
2004; Isaza-M, 2007). Estos com- 
puestos han sido objeto de estudio, 
evaluando su contenido en cortezas de 
pino y encino (Rosales-Castro y Gon- 
zález-Laredo, 2003), y recientemente 
valorando los compuestos fenólicos y 
propiedades nutracéuticas de hojas de 
encinos de Durango (Rocha-Guzman et 
al., 2009). Existe un interés creciente a 
nivel mundial en estudiar los com- 
puestos polifenólicos, ya que son 
capaces de inhibir procesos oxidativos 
asociados con alteraciones de la salud 
(Cala-Molina et al., 2007). Se han deter- 
minado los compuestos fenólicos totales 
en diversos tejidos vegetales, por 
ejemplo, las hojas secas de Coccinia 
grandis y Cymbopogon citratus con- 
tienen 1,6 y 1,0 mg/g de compuestos 
fenólicos totales, respectivamente; los 
frutos secos de Psidium guajava y 
Cucúrbita pepo contienen 0,60 y 0,20 
mg/g de compuestos fenólicos totales, 
respectivamente (Mongkolsilp et al., 
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2004). Los frutos secos de Prunos 
pérsica, Phaseolus vulgaris y Ficus 
carica contienen 0,50; 0,35 y 0,59 mg/g 
de compuestos fenólicos totales, 
respectivamente; y los tallos secos de 
Allium cepa y Allium porrum contienen 
0,36 y 0,27 mg/g de compuestos fenó- 
licos totales, respectivamente (Marinova 
et al., 2005). 

Los compuestos polifenólicos son 
sintetizados a partir de esqueletos de 
carbono, al igual que los rebrotes jóvenes 
que también necesitan y demandan 
reservas de carbono para su formación y 
crecimiento (Masashi et al., 2007). Esta 
competencia incide en la formación de 
follaje, de ahí la importancia de la medi- 
ción del contenido fenólico en los 
encinos, a mayor contenido fenólico 
menor producción de rebrotes. 

La pregunta de investigación que se 
planteó fue si había un mejor suelo para 
la plantación de los Quercus sideroxyla y 
la forma en que el factor suelo afectaba el 
contenido fenólico foliar, ya que el tipo de 
suelo se ha relacionado siempre con los 
nutrientes disponibles que asimila la 
planta y, por lo tanto, diferentes tipos de 
suelo podrían incidir en el contenido fenó- 
lico foliar de los encinos. 

La importancia de este estudio 
radica en la probable relación inversa 
entre la formación de rebrotes en 
Quercus sideroxyla y el contenido fenó- 
lico foliar. 

OBJETIVO 

Por lo anterior el objetivo planteado es 
determinar la variación estacional en el 
contenido fenólico foliar del encino rojo 
{Quercos sideroxyla), como un parámetro 
del deterioro del género, para que pueda 
servir de base en la implementación de 
una plantación comercial. 
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METODOLOGÍA 

El estudio se llevó a cabo en el Predio 
“Las Bayas”, localizado en el municipio de 
Pueblo Nuevo, Durango, a 100 km al sur 
de la ciudad de Durango, a una altitud de 
2500 msnm a 2900 msnm, por el camino 
Durango-La Flor. El predio cuenta con 
5 068 ha, de las cuales 3 200 ha se 
encuentran en producción, y está dividido 
en 36 rodales, de los que se seleccio- 
naron 8 subrodales ubicados en 6 rodales 
diferentes como área de estudio con 
centro 23° 28’ Norte, 104° 51’ Oeste. 

Muestreo foliar. Se realizó muestreo alea- 
torio estratificado, donde los estratos fueron 
las unidades de manejo donde se presen- 
taba la especie de interés, Quercus side- 
roxyla. El muestreo se realizó en sitios altos, 
por encima de los 2700 msnm, durante las 
cuatro estaciones del año; se inició con el 
invierno (5 y 6 de febrero, 2008), luego 
primavera (23 y 24 de mayo), después 
verano (20 y 21 de agosto) y finalmente 
otoño (4 y 5 de diciembre). El tamaño de los 
sitios de muestreo fue de un décimo de 
hectárea (1000 m2) en forma circular. 
Siempre se tomó como punto de referencia 
el árbol de mayor tamaño. Sólo se utilizaron 
árboles de apariencia sana y de la misma 
edad, con una altura promedio de 7 m a 9 m. 
Para calcular la edad de los árboles se esti- 
maron los parámetros, diámetro a la altura 
del pecho y altura total del árbol (Reed etal., 
1992). El muestreo de hojas se realizó en la 
parte media del árbol, evitando colectar 
hojas exteriores o en activo crecimiento. Las 
hojas se colocaron en bolsas de polietileno, 
se trasladaron en hielo y se mantuvieron en 
refrigeración en el laboratorio. Para la prepa- 
ración del extracto se analizaron hojas de 
cada individuo de manera separada, éstas 
se secaron bajo la sombra, a temperatura 
ambiente durante siete días y se molieron 
hasta tener un polvo fino homogéneo. 

Extracción de fenoles. Se tomó polvo de 
hoja seca (1 g) y se agregaron 10 mL de 
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metanol (50%). Se dejó reposar (1 h), se 
filtró y a 10 pL de extracto se le agregó 
1 mL de metanol (50%). 

Contenido fenólico foliar. El análisis de 
compuestos fenólicos totales se realizó 
por el método de Folin-Ciocalteu modifi- 
cado por Dewanto et al. (2002), utilizando 
ácido tánico como estándar. La cantidad 
de compuestos fenólicos totales se 
expresa como equivalentes de ácido 
tánico (eat, mg de ácido tánico/g de 
muestra), a través de una curva de cali- 
bración de ácido tánico. La curva de cali- 
bración tuvo rango de (20-120) mg/mL 
(R2 = 0,9996). El límite de detección del 
método fue de 20 mg/mL de ácido tánico 
que corresponde a 0,004 unidades de 
absorbencia. 

Este método espectrofotométrico 
consiste en agregar 0,1 mL de muestra a 
un matraz de 10 mL que contenga 6 mL 
de agua desionizada, luego se añade 
0,5 mL del reactivo comercial de Folin- 
Ciocalteu. Después de 60 s se agrega 1 ,5 
mL de Na 2 C 03 (20%), se afora a 10 mL 
con agua desionizada, se mezcla y 
después de 2 h se mide la absorbancia a 
760 nm. Para cada determinación se 
realizaron tres repeticiones. Para este 
estudio se utilizó un espectrofotómetro 
marca génesis. 

Muestreo del suelo. La toma de muestra 
se realizó a 25 cm de profundidad. Cada 
muestra correspondió a la ubicación de 
un árbol, de acuerdo con el método 
descrito por Rodríguez-Fuentes (2002). 

pH del suelo. Se utilizó el método as-02 
señalado en la nom-02i-semarnat, mane- 
jando un potenciómetro manual oakton 
CE. Para cada determinación se reali- 
zaron tres repeticiones. 

Conductividad eléctrica del suelo. Se 
utilizó una bomba de vacío marca felisa, 
con papel filtro de poro fino Whatman 42, 


de acuerdo con el método descrito por 
Rodríguez-Fuentes (2002). Para cada 
determinación se realizaron tres repeti- 
ciones. 

Retención de humedad del suelo. Se 
utilizó un horno marca arsa, según el 
método gravimétrico descrito por Rodrí- 
guez-Fuentes (2002). Para cada determi- 
nación se realizaron tres repeticiones. 

Textura del suelo. Se utilizó el método de 
Boyoucos (Rodríguez-Fuentes, 2002). 

Los datos registrados se analizaron 
estadísticamente utilizando el proc ttest 
y el PR 0 C_C 0 RR procedimiento del 
SAS/STAT® (SAS Institute Inc., 2004). 

Formulación de hipótesis. Se estable- 
cieron como base de la investigación 
cuatro hipótesis que se enuncian de la 
siguiente manera: 

Hipótesis 1. El contenido fenólico 
foliar varía en función de las caracterís- 
ticas del suelo, estacionalmente. 

Hipótesis 2. Existe una relación 

significativa entre el contenido fenólico 
foliar y la conductividad eléctrica del 
suelo. 

Hipótesis 3. Existe una relación 

significativa entre el contenido fenólico 
foliar y la retención de humedad del 

suelo. 

Hipótesis 4. Existe una relación 

significativa entre el contenido fenólico 
foliar y la acidez del suelo (pH). 

Operacionalización 

Para la operacionalización de las varia- 
bles en la Hipótesis 1 se estableció el tipo 
de suelo y se determinó el contenido 
fenólico foliar promedio (mg/g) en mués- 
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Tabla 1. Determinación de variables por tipo de suelo. 



VARIABLE INDEPENDIENTE 

VARIABLE DEPENDIENTE 

Tipo de suelo 

Estación 

Contenido fenólico foliar 

Determinado 

experimentalmente 

Invierno 

Primavera 

Verano 

Otoño 

Determinado 

experimentaimente 

Tabla 2. Operacionalización de Hipótesis 2, 3 y 4. 


VARIABLES INDEPENDIENTES 

VARIABLE DEPENDIENTE 

TIPO DE SUELO 

Conductividad eléctrica 
Retención de Humedad 
pH 

Contenido fenólico foliar 


tras tomadas en las diferentes estaciones 
del año (Tabla 1). 

La operacionalización de las hipótesis 
subsecuentes se realizó determinando 
conductividad eléctrica, retención de 
humedad, pH del suelo y estimando el 
contenido fenólico foliar de las muestras 
tomadas en las diferentes estaciones del 
año, en función del tipo de suelo (Tabla 2). 

Para la verificación de las hipótesis 
2, 3 y 4 se determinó la relación utilizando 
el coeficiente de correlación de Pearson 
(rx) entre las variables involucradas en 
cada una de las hipótesis, relacionando: 

1. El contenido fenólico foliar y la 
conductividad eléctrica. 

2. El contenido fenólico foliar y la 
retención de humedad. 

3. El contenido fenólico foliar y el pH 
del suelo. 


RESULTADOS 

Resultados por variable 

Tipo de suelos. A partir de la toma de 
muestras se determinó la presencia de 
tres tipos de textura (Tabla 3). 


Tabla 3. Tipo de suelo en muestras 
analizadas. 


Tipo de suelo 

% 

Franco arcilloso 

12,5 

Limo arcilloso 

75 

Arcillo arenoso 

12,5 


Contenido fenólico foliar. Los resultados 
obtenidos estacionalmente por tipo de 
suelo se señalan en la tabla 4 y figura 1 . 

El contenido fenólico foliar de 
Quercus sideroxyla varía de forma dife- 
renciada, en función del tipo de suelo a lo 
largo del año, observándose un compor- 
tamiento semejante en el contenido foliar 
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Variación estacionai de compuestos fenólicos foiiares en Quercus sideroxila 


Tabla 4. Contenido fenólico foliar* 


Textura 

Invierno 

Primavera 

Verano 

Otoño 

Promedios 


mg/g de hoja seca 

mg/g de hoja seca 

mg/g de hoja seca 

mg/g de hoja seca 

FRANCO ARCILLOSA 

0,271 

0,236 

0,274 

0,256 

0,259 

LIMO ARCILLOSA 

0,264 

0,233 

0,258 

0,249 

0,251 

ARCILLO ARENOSA 

0,227 

0,226 

0,281 

0,242 

0,244 


‘medias aritméticas 


CONTENIDO 

FENÓÜCO 

FOÜAR 



INVIERNO PRIMAVERA VERANO OTOÑO PROMEDIO 
FRANCO ARCILLOSA -<*«LiMO ARCILLOSA ARCILLO ARENOSA 


Figura 1. Contenido fenólico foliar* estacionalmente por tipo de suelo. 

‘medias aritméticas 


de los encinos en suelos franco arcilloso 
y limo arcilloso, manteniéndose siempre 
el limo arcilloso por debajo del franco arci- 
lloso con valores semejantes durante la 
primavera. Esto no ocurre en tos encinos 
de suelos arcillo arenosos que presentan 
un bajo contenido fenólico foliar en 
invierno y primavera, lo que garantiza la 


formación de rebrote, aunque durante el 
verano es el que presenta el valor más 
alto de contenido fenólico foliar que dismi- 
nuye radicalmente en el otoño para volver 
a colocarse por debajo de los otros dos 
tipos de suelo, situación que no se refleja 
en los valores promedio. 
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Conductividad eléctrica. Los resultados 
obtenidos estacionalmente por tipo de suelo 
se presentan en la tabla 5. Independiente- 
mente de la textura de los suelos, las 
conductividades más bajas se encuentran 
en el invierno, debido al menor contenido de 
humedad y a un menor metabolismo de la 
planta, mientras que lo contrario ocurre en 
el verano cuando hay mayor humedad y 
temperaturas más altas que incrementan el 
metabolismo vegetal. 

Humedad del suelo. Los resultados obte- 
nidos estacionalmente por tipo de suelo 
se señalan en la tabla 6. Para cualquier 
tipo de suelo la retención de humedad es 
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menor de 20% en la primavera y mayor 
de 40% en el verano, en las estaciones 
de invierno y verano no se presentó una 
diferencia significativa en la retención de 
humedad entre los diferentes suelos. 

Acidez del suelo. Los resultados obte- 
nidos estacionalmente por tipo de suelo 
se señalan en la tabla 7. Aquí el pH se 
presenta extremadamente ácido en 
invierno, aumentando hasta moderada- 
mente ácido en el otoño, para los suelos 
limo arcilloso y arcillo arenoso. Aunque el 
pH del suelo franco arcilloso incrementa 
casi en la misma proporción, no logró 
llegar al rango de moderadamente ácido. 


Tabla 5. Conductividad del suelo* 


Textura 

Invierno 

Primavera 

Verano 

Otoño 

Promedios 

FRANCO ARCILLOSA 

31,9 pS/cm 

36,5 pS/cm 

38,3 pS/cm 

37,4 pS/cm 

36,0 pS/cm 

LIMO ARCILLOSA 

25,6 pS/cm 

40,7 pS/cm 

43,7 pS/cm 

30,8 pS/cm 

35,2 pS/cm 

ARCILLO ARENOSA 

20,1 pS/cm 

36,7 pS/cm 

40,4 pS/cm 

27,6 pS/cm 

31 ,2 pS/cm 

‘medias aritméticas 







Tabla 6. 

Humedad del suelo* 



Textura 

Invierno 

Primavera 

Verano 

Otoño 

Promedios 

FRANCO ARCILLOSA 

28,82% 

12,25% 

41,21% 

28,16% 

27,61 % 

LIMO ARCILLOSA 

28,90% 

13,23% 

41 ,47% 

24,52% 

27,03 % 

ARCILLO ARENOSA 

30,59% 

19,44% 

40,38% 

32,18% 

30,64 % 

‘medias aritméticas 







Tabla 7. pH del suelo* 



Textura 

Invierno 

Primavera 

Verano 

Otoño 

Promedios 

FRANCO ARCILLOSA 

4,0 

4,1 

4,3 

5,2 

4,4 


EA 

EA 

EA 

FA 

EA 

LIMO ARCILLOSA 

4,3 

4,6 

4,7 

5,6 

4,8 


EA 

FA 

FA 

MO 

FA 

ARCILLO ARENOSA 

4,0 

4,3 

4,5 

5,5 

4,5 


EA 

EA 

FA 

MO 

FA 


‘medias aritméticas 

EA = extremadamente ácido, FA = fuertemente ácido, MO = moderadamente ácido 
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Variación estacional de compuestos fenólicos foliares en Quercus sideroxila 


Evaluación de las hipótesis 

Los resultados de la evaluación de la 
Hipótesis 1, considerando el contenido 
fenólico foliar promedio por tipo de suelo, 
comparado estacionalmente (invierno- 
primavera; primavera-verano y verano- 
otoño) se presenta en la tabla 8. 


Se encontró que existen diferencias 
estadísticamente significativas en el 
contenido fenólico foliar, en todos los 
suelos, en las diferentes estaciones del 
año, con excepción del suelo arcillo 
arenoso en invierno-primavera, en donde 
no hay diferencias estadísticamente signi- 
ficativas en el contenido fenólico foliar. 


Tabla 8. Contenido fenólico foliar estacional por tipo de suelo. 


Tipo de suelo 


Hipótesis 

Ho 


p valué 

^ O 

O C/D 

Inv-Prim 

Hq,i 

Xi = X2 

0,0001 

a < 0,05 

ü O 
Z -I 
< ^ 

Prim-Ver 

Hq,2 

X 2 = X 3 

0,0001 

a < 0,05 

01 O 
LL 01 
< 

Ver-Oto 

Hq ,3 

X 3 =X 4 

0,0001 

a < 0,05 

O 

(/) 

o o 

Inv-Prim 

Hq ,4 

X 5 = X 6 

0,0001 

a < 0,05 

S d 

ü 

Prim-Ver 

Ho ,5 

X e = X 7 

0,0001 

a < 0,05 

c¿. 

< 

Ver-Oto 

Ho,6 

X 7 = Xs 

0,0039 

a < 0,05 

° o 
3 o 

Inv-Prim 

Hq ,7 

X 9 “ ^10 

0,5178 

a > 0,05 

-1 Z 

ü 

(Y CL 

Prim-Ver 

Hq,8 

X 10 = X 11 

0,0001 

a < 0,05 

< < 

Ver-Oto 

Hq ,9 

X 11 = X 12 

0,0001 

a < 0,05 


Tabla 9. 

Correlación de Variable del Contenido Fenólico Foliar con Conductividad 
Eléctrica, Retención de Humedad y pH. 

FRANCO ARCILLOSO 


VARIABLE 

CONDUCTIVIDAD 

RETENCIÓN 

pH 




DE HUMEDAD 


CONTENIDO 

rx 

-0,149 

0,842** 

-0,076 

FENÓLICO 

A 

0,486 

< 0,0001 

0,720 

FOLIAR 

SD 

4,72 

10,83 

0,48 


N 

24 

24 2 

4 

LIMO ARCILLOSO 


VARIABLE 

CONDUCTIVIDAD 

RETENCIÓN 

pH 




DE HUMEDAD 


CONTENIDO 

rx 

-0,023 

0,474** 

-0,038 

FENÓLICO 

A 

0,783 

< 0,0001 

0,650 

FOLIAR 

SD 

8,04 

10,48 

0,49 


N 

144 

144 1 

44 

ARCILLO ARENOSO 


VARIABLE 

CONDUCTIVIDAD 

RETENCIÓN 

pH 




DE HUMEDAD 


CONTENIDO 

rx 

0,610** 

0.734** 

0,220 

FENÓLICO 

A 

0,001 

< 0,0001 

0,301 

FOLIAR 

SD 

8,25 

8,13 

0,57 


N 

24 

24 

24 
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Hipótesis 2, 3 y 4. A partir de ios resultados 
reportados en la tabla 9, se puede afirmar que 
hay una correlación entre el contenido fenó- 
lico foliar y la retención de humedad en los 
tres tipos de suelo (muy alta para los suelos 
franco arcilloso y arcillo arenoso). En tanto, 
que no se presenta correlación con pH en el 
rango en que se encontraron estos suelos 
(pH 4,0 a 5,6): en el caso de la relación de 
contenido fenólico y la conductividad, sólo es 
significativa para el suelo arcillo arenoso. 

DISCUSIÓN Y CONCLUSIONES 

En este artículo se ha analizado estacional- 
mente la producción fenólica foliar, como 
un parámetro de deterioro del género 
Quercos, observándose de forma general 
que, independientemente de la textura de 
los suelos, el contenido fenólico foliar fue 
más bajo en primavera, ganando en la 
competencia de las reservas energéticas la 
producción de hojas nuevas sobre la gene- 
ración de compuestos fenólicos foliares 
(Masashi et ai., 2007). 

Como ha sido observado por otros 
investigadores, el contenido fenólico foliar 
aumenta durante el verano independien- 
temente de la textura del suelo, lo que se 
atribuye a que la mayoría de estos 
compuestos cumplen funciones de 
defensa contra depredadores, patógenos 
e irradiación solar (Salminen et ai., 2004), 
característicos en esta estación; mientras 
que su disminución durante el otoño se 
atribuye a la moderación de esos 
factores. 

El mejor suelo será entonces el que, 
por un lado, promueva el menor contenido 
fenólico foliar durante la primavera para 
permitir un mayor número de rebrotes y que, 
por otro lado, produzca mayor cantidad de 
compuestos polifenólicos para protección 
del árbol en el verano. Bajo estas circunstan- 
cias se puede decir que el mejor suelo es el 
de textura arcillo arenosa, por ser el que 
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genera en el encino la menor producción de 
compuestos polifenólicos en primavera y la 
mayor en verano en relación con los otros 
tipos de suelo, además de que presenta el 
mayor contenido de humedad en invierno, 
primavera y otoño, siendo más o menos 
equivalente en el verano como resultado de 
precipitaciones abundantes y en donde la 
correlación de conductividad significativa 
pareciera ser una variable asociada al 
contenido fenólico foliar. 
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lARTlCULO DE INVESTIGACION! 


Anatomía de la madera de las especies 
arbóreas de Ipomoea (Convolvulaceae) 

Wood anatomy of the arboreal species 
of Ipomoea (Convolvulaceae) 

Jacqueline Ceja-Romero'' y Carmen de la Paz Pérez-Olvera'' 

RESUMEN 


Se describe la anatomía de la madera de ocho especies de la serie Arborescentes del género 
Ipomoea, grupo con porte arbóreo o de enredadera leñosa y flores de color blanco o amarillo. Las 
especies presentan cambia sucesivos que forman anillos concéntricos de xilema, floema y tejido 
conjuntivo. Todos los taxa tienen vasos solitarios a frecuentemente agrupados, redondeados a elíp- 
ticos y con diámetro tangencial medio; los elementos de vaso son cortos, con la placa de perforación 
simple y las punteaduras alternas y poligonales. Presentan tílides. Las fibrotraqueidas van de muy 
cortas a medianas en su mayoría con un diámetro medio y paredes muy delgadas. El parénquima 
axial es de tipo vasicéntrico y ocasionalmente en bandas angostas. El parénquima radial está formado 
por radios uniseriados y multiseriados, heterogéneos, con granos de almidón y drusas. Los resultados 
coinciden con lo publicado para el grupo, pero se presentan diferencias estadísticas en la longitud y 
el diámetro medio de las fibras, en la altura y el número de células de los radios uniseriados y en la 
altura y número de series de los radios multiseriados. 

PALABRAS CLAVE: 

Cámbium sucesivo, Ipomoea, serie Arborescentes, variante cambial, xilema secundario. 


ABSTRACT 


Wood anatomy of eight species of the series Arborescentes of Ipomoea is described. This group 
has arboreal or viny habit and white or yellow flowers. The species have a type of cambial variant 
which consists of concentric rings of xyiem, phloem and conjunctive tissue. Taxa have solitary pores 
or sometimos they are grouped, rounded to elliptical and with médium tangential diameter; vessel 
elemente are short with simple perforated platee, pits are altérnate, polygonal. Tyloses are present. 
Fibre-tracheids have a médium length and diameter, and thin walls. Axial parenchyma is vasicentric 
although occasionally narrow bands are present. Radial parenchyma is formed by uniseriate to multi- 
seriate raye, heíerogeneous. Abundant starch grains and druses are present. Resulte are consistent 
with those found in the literature for the group, but there are statistical differences in length and 
diameter of fibres, height and cell number of uniseriate raye and height and series number of multise- 
riate raye. 

KEY WORDS: 

Successive cambium, Ipomoea, series Arborescentes, cambial variant, secondary xyiem. 


1 


Universidad Autónoma Metropolitana-Iztapalapa, División de C.B.S., Departamento de Biología. Área de 
Botánica Estructural y Sistemática Vegetal. AP. 55-535, 09340, México, D.F. jcr@xanum.uam.mx. 


62 


Anatomía de la madera de las especies arbóreas de Ipomoea (Convolvulaceae) 


INTRODUCCIÓN 

La familia Convolvulaceae está integrada 
por 40 a 50 géneros con aproximada- 
mente 1600 especies de distribución 
cosmopolita, pero se localizan preferente- 
mente en las regiones tropicales y subtro- 
picales del mundo. En México se tienen 
registrados hasta el momento 16 géneros 
y 217 a 260 especies (Austin y Pedraza, 
1983; McDonald, 1991). Dentro de la 
familia, el taxón mejor representado es 
Ipomoea, con cerca de 600 especies, de 
las cuales 150 crecen en México (Austin y 
Pedraza, 1983; Miller et al., 1999). 

Ipomoea incluye plantas principal- 
mente herbáceas, aunque también hay 
sufrútices y leñosas; postradas, trepa- 
doras o erectas; con hojas enteras, 
dentadas o lobadas a digitado compues- 
tas. Las flores son solitarias, dispuestas 
en inflorescencias racemosas o cimosas, 
con la corola generalmente infundibuli- 
forme y los frutos capsulares (Matuda, 
1964; Austin y Pedraza, 1983; Standiey, 
1982). Debido a la diversidad morfológica 
entre las especies Austin (1979, 1980) 
propuso, con base en las características 
del gineceo y de las semillas, tres subgé- 
neros: Ipomoea, Quamoclit y Erios- 
permum, cada uno dividido en secciones 
y series, sin embargo, la clasificación del 
género sigue siendo objeto de discusión, 
señalándose en diversos trabajos (Austin 
y Huáman, 1996; Manos et al., 2001) la 
necesidad de revisarla no solo a nivel 
infragenérico sino también de géneros 
relacionados. 

Para las Ipomoea americanas, 
Austin (1979) y Austin y Huáman (1996) 
abordan un tratamiento infragenérico 
completo, incluyendo en el subgénero 
Eriospermum a taxa perennes de hábito 
variado, sépalos coriáceos, semillas 
pilosas y gineceos con dos lóculos. 
Dentro de este subgénero hay un grupo 
de especies, que se distingue del resto 


por su hábito arbóreo o de enredadera 
leñosa y por sus flores blancas, a veces 
algo rosadas o amarillas (figuras 1 y 2). 
Taxonómicamente la variabilidad en las 
características morfológicas de este 
conjunto de taxa, hace que existan dife- 
rentes opiniones en cuanto a su delimita- 
ción y número. Matuda (1964) los ubica 
dentro del “grupo Arbóreas", Austin (1979, 
1980) y Austin y Huáman (1996) los reco- 
nocen bajo la serie Arborescentes de la 
sección Eriospermum del subgénero del 
mismo nombre, McPherson (1981) los 
denomina como “el grupo de Ipomoea 
arborescens" y reconoce 10 especies, de 
las cuales nueve se distribuyen en 
México y una es endémica de Perú. 
McDonald (1991), en su estudio de las 
Convolvulaceae de México, coincide en el 
número de especies propuesto por 
McPherson, aunque difiere en cuanto a la 
delimitación de las mismas. Carranza et 
al. (1998) añaden al grupo una nueva 
especie del centro del país, quedando, 
hasta el momento, constituido por 11 
especies: Ipomoea arborescens (Humb. 
& Bonpl. ex Wild.) G. Don, /. chilopsidis 
Standiey, /. intrapilosa Rose, I. muru- 
coides Roem. & Schult., /. pauciflora M. 
Martens & Galeotti, I. populina House, 
/. praecana House, /. pulcherrima van 
Ooststroom, /, teotitlanica McPherson, 
/. rzedowskii S. Carranza, S. Zamudio & 
G. Murguía e /. wolcottiana Rose, En 
México a varias de ellas (/. arborescens, 
I. intrapilosa, I. murucoides, I. pauciflora e 
I. wolcottiana) se les da el nombre común 
de “cazahuate” y dado que son un 
elemento predominante de algunas 
selvas bajas, llegan a formar verdaderos 
“cazahuatales”. 

El estudio de las especies que inte- 
gran a la serie Arborescentes, ha sido 
abordado por disciplinas como la palino- 
logía, la taxonomía y la histoquímica 
(Pedraza, 1983; Pérez-Amador et al., 
1983, 1992). Desde el punto de vista 
anatómico se han estudiado las flores y 
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las semillas de algunas especies 
(Murguía-Sánchez, 1986,1995; Murguía- 
Sánchez et al., 1995) y en lo que se 
refiere a la anatomía del tallo se pueden 
citar los trabajos de Luján (1974), que 
describe el origen y desarrollo del creci- 
miento secundario en dos especies del 
centro de México (/. arborescens e /. 
pauciflora)] el de Cariquist y Hanson 
(1991), que estudian la anatomía de este 
órgano dentro de la familia, incluyendo 
tres de las especies del grupo de las 
arborescentes (/. pauciflora, I. muru- 
coides e I. arborescens) y el de McDonald 
(1992), que analiza las implicaciones 
evolutivas que la estructura interna de 
este grupo de especies tiene con 
respecto al resto de la familia, descri- 
biendo a detalle la anatomía /. muru- 
coides, como representante de las espe- 
cies con crecimiento secundario anómalo. 

OBJETIVO 

Realizar el estudio anatómico compara- 
tivo de la madera de: Ipomoea arbores- 
cens, I. chilopsidis, I. intrapilosa, I. muru- 
coides, I. pauciflora, I. populina, I. rze- 
dowskii e I. wolcottiana, grupo de taxa 
que integran a la serie Arborescentes de 
la sección Eriospermum del subgénero 
del mismo nombre, según la clasificación 
propuesta por McDonald (1991) para las 
Convolvuláceas mexicanas. 

MATERIALES Y METODOLOGÍA 

Se hicieron diversos viajes de recolección 
a lugares previamente seleccionados 
(Tabla 1). En cada localidad se esco- 
gieron dos individuos representativos por 
especie, de los que se obtuvieron mues- 
tras de la zona cercana a la corteza, a 
una altura de 1,30 m del fuste en los 
árboles y de la rama más gruesa en los 
arbustos, las cuales se fijaron en una 
solución de formalina-ácido acético- 
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etanol (faa, Sass, 1958). En el caso de 
los taxa con hábito arbóreo, también se 
recolectaron ramas con un diámetro de 
alrededor de 10 cm, con el fin de tener 
tablillas para la xiloteca. De todos los indi- 
viduos se hicieron los respectivos ejem- 
plares de respaldo, quedando estos 
depositados en el Herbario Metropolitano 
“Ramón Riba y Nava Esparza” (uamiz). 

Con el material fijado de árboles y 
arbustos, se hicieron cubos de 2 cm x 2 
cm; en el caso de las lianas, debido a su 
diámetro menor, los cubos fueron de 1 cm 
X 1 cm. De las muestras se obtuvieron 
cortes (transversales, radiales y tangen- 
ciales) de (20-30) pm de grosor, utilizando 
para ello un micrótomo de deslizamiento 
American Optical. Las secciones fueron 
deshidratadas con alcohol etílico a 
concentraciones graduales, teñidas con 
safranina-verde rápido (Sass, 1958) y 
montadas con resina sintética. El material 
disociado se obtuvo a partir de un solo 
derivado cambial utilizando la solución de 
Jeffrey (Johansen, 1940), los elementos 
se tiñeron con Safranina O y se montaron 
en gelatina glicerinada. 

La determinación de los caracteres 
macroscópicos se hizo de acuerdo con 
las tablas de Tortorelli (1956). El color se 
determinó con las tablas de Munsell 
Color Company (1990). Las caracterís- 
ticas microscópicas se determinaron a 
partir de las preparaciones permanentes 
de los cortes y del material disociado, 
tomando para cada carácter un total de 
25 mediciones. La denominación de 
cada uno de los caracteres se hizo con 
base en la media según la clasificación 
de Chattaway (1932) y del iawa 
Committee (1937, 1939 y 1989). Para 
cada una de las características cuantita- 
tivas (Tabla 2) se realizó un análisis de 
varianza y una prueba de Tukey, utili- 
zando el programa ncss 2001 (Number 
Cruncher Statistical System for 
Windows, Kaysville Utah). 
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Anatomía de la madera de las especies arbóreas de Ipomoea (Convolvulaceae) 


Tabla 1 . Material recolectado y procesado de Ipomoea serie Arborescentes. 


Especie 

Núm. de Colecta Registro de 

herbario 

Estado 
del tallo (cm) 

Hábito 

Diámetro 

1. arborescens 

JCR 1232 

UAMIZ 67052 

Morelos 

árbol 

25 


JCR 1234 

UAMIZ 67050 

Morelos 

árbol 

24 


JCR 1242 

UAMIZ 67048 

Morelos 

árbol 

25 


JCR 1422 

UAMIZ 67037 

Colima 

árbol 

20 

1. chilopsidis 

JCR 1702 

Sin registro 

Chihuahua 

arbusto 

6 


JCR 1706 

UAMIZ 65357 

Chihuahua 

arbusto 

8 

1. intrapilosa 

JCR 1318 

UAMIZ 67054 

Jalisco 

árbol 

25 


JCR 1319 

UAMIZ 67046 

Jalisco 

árbol 

25 

1. murucoides 

JCR 1237 

UAMIZ 67044 

Morelos 

árbol 

20 


JCR 1257 

UAMIZ 67053 

Oaxaca 

árbol 

25 


JCR 1320 

UAMIZ 67045 

Michoacán 

árbol 

25 

1. pauciflora 

JCR 1240 

UAMIZ 67056 

Morelos 

árbol 

25 


JCR 1316 

UAMIZ 67041 

Puebla 

árbol 

20 


JCR 1259a 

UAMIZ 67035 

Oaxaca 

árbol 

22 

1. populina 

JCR 1059 

UAMIZ 57056 

Guerrero 

enredadera 

1.5 


JCR 400bis 

UAMIZ 67038 

Morelos 

enredadera 

1.5 

1. rzedowskii 

JCR 1246 

UAMIZ 67047 

Querétaro 

arbusto 

10 

1. wolcottiana 

JCR 1565 

UAMIZ 67065 

Oaxaca 

árbol 

20 


JCR 1266 

UAMIZ 67033 

Oaxaca 

árbol 

20 


RESULTADOS 

Características macroscópicas 

En el tallo de todas las especies se 
presenta un crecimiento secundario con 
variante cambial (cambia sucesivos), que 
da lugar a cilindros concéntricos com- 
puestos por bandas de xilema, de color 
blanco y de 1 mm a 2 mm de grosor, que se 
alternan con bandas de floema y tejido 
conjuntivo, de color castaño y de cerca de 
1 mm de grosor (Figura 3). No hay diferencia 
entre la albura y el duramen, la madera no 
presenta olor ni sabor característicos, el 
brillo es bajo, el veteado va de suave a 
pronunciado, la textura va de mediana 
a gruesa y el hilo es recto (figuras 4 y 5). 

Características microscópicas 

Xilema secundario. Forma un cilindro 
concéntrico que en la mayoría de las 


especies se interrumpe dando la apa- 
riencia de fascículos (Figura 6) o que en 
el caso de I. chilopsidis e /. rzedowskii es 
continuo (Figura 16). Los vasos general- 
mente se diferencian cercanos al anillo 
cambial, cuando el xilema ha alcanzado 
su máxima acumulación, son solitarios y 
la mayoría agrupados, redondos a elíp- 
ticos transversalmente (Figura 6), son 
medianos con un diámetro tangencial 
promedio de 114 pm a 142 pm y paredes 
de 4 pm a 5 pm de grosor. Los elementos 
de vaso son cortos, con una longitud 
promedio de 130 pm a 163 pm; las colas 
son frecuentes; la placa de perforación es 
simple; las punteaduras son areoladas, 
alternas, de forma poligonal, con la aber- 
tura tangencialmente alargada (Figura 7). 
Los análisis de varianza indicaron que no 
hay diferencias significativas en los carac- 
teres cuantitativos de este tipo celular. 

En todas las especies se observan 
fibrotraqueidas y son el elemento predo- 
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minante en el xilema secundario, son 
cortas en I. populina, I. chilopsidis e I. 
rzedowskii y medianas en el resto de las 
especies, con una longitud promedio de 
688 pm en /. populina a 1116 pm en I. 
arborescens, su diámetro medio va de 12 
pm a 25 pm, quedando en la categoría de 
finas, las paredes son delgadas con 4 pm 
a 5 pm de grosor (Figura 8). El análisis de 
varianza señaló diferencias significativas 
en la longitud de este elemento y la 
prueba de Tukey indicó que /. chilopsidis 
e /. populina son estadísticamente dife- 
rentes de todos los taxa con excepción de 
/. rzedowskii, la cual a su vez difiere de I. 
murucoídes, I. wolcottiana e I. arbores- 
cens y esta última difiere de todas las 
demás. Lo anterior se puede relacionar 
con el hecho de que /. chilopsidis, I. popu- 
lina e /. rzedowskii queden en la categoría 
de cortas, sin embargo, en el caso de L 
arborescens queda junto con las demás 
en la categoría de medianas, no refleján- 
dose la diferencia en la clasificación. 
También en el diámetro medio hay dife- 
rencias entre I. populina e I. chilopsidis 
con respecto a todas las demás especies 
y de /. arborescens e L murucoides con 
I. intrapilosa, e /. wolcottiana. 

El parénquima axial es de tipo vasi- 
céntrico, aunque ocasionalmente se 
presenta en bandas apotraqueales 
delgadas. Las células son elipsoides a 
oblongas, con paredes delgadas, poco o 
no lignificadas. Es frecuente que las 
células del parénquima adyacente 
puedan penetrar el lumen del vaso, 
dando lugar a tílides formadas por 
pequeñas células no lignificadas. 

Los radios se presentan uniseriados 
y multiseriados (figuras 10 y 11). Los 
uniseriados están formados por 3 a 8 
células erectas, con 6 a 8 radios por milí- 
metro lineal y con una altura promedio 
178 pm en /. populina a 302 pm en 
/. chilopsidis. Los análisis numéricos 
indican diferencias en el número de 
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células por radio entre /. chilopsidis y las 
demás especies y en la altura de I. popu- 
lina con respecto a I. arborescens, I. intra- 
pilosa, I. murucoides, I. wolcottiana e I. 
chilopsidis y esta última con /. rzedowskii. 

Los radios multiseriados (Figura 11) 
van de 6 a 10 series, con 156 pm de ancho 
en /. chilopsidis a 213 pm en /. wolcottiana, 
son muy bajos en I. chilopsidis e I. populina 
y bajos en el resto de las especies, con 
intervalos de altura que van de 746 pm en I. 
chilopsidis a 1226 pm en I. wolcottiana, son 
heterogéneos y están formados por células 
erectas, cuadradas y, en menor cantidad, 
procumbentes, mezcladas a lo largo y 
ancho de éstos (Figura 12), presentan 
abundantes granos de almidón y drusas 
(figuras 13 y 14). Los resultados indican 
diferencias estadísticas en el número de 
series de I. chilopsidis e I. populina con I. 
arborescens e I. wolcottiana mientras que 
en la altura /. chilopsidis es significativa- 
mente diferente de todas las demás espe- 
cies con excepción de I. populina, mientras 
que esta última difiere solo de /. wolcottiana 

Floema secundario. En un tipo de 
variante cambial como la que se presenta 
en este caso, es difícil dejar a un lado los 
elementos que constituyen al anillo en 
conjunto, por lo que a continuación se 
describen brevemente el floema y el 
tejido conjuntivo. 

El floema está compuesto por 
elementos de tubo criboso, células acom- 
pañantes y abundantes radios unise- 
riados. Los elementos de tubo criboso 
presentan áreas cribosas compuestas 
sobre las paredes laterales y placas 
cribosas simples o compuestas (Figura 
15), con una o más células acompa- 
ñantes que pueden alcanzar la misma 
longitud que estos. 

Laticíferos. Se presentan en el floema 
(Figura 9) y son del tipo articulado no 
anastomosante. 
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Anatomía de la madera de las especies arbóreas de Ipomoea (Convolvulaceae) 



Figuras 1-8. 1. Hábito de la planta {Ipomoea murucoides, J. Ceja s. n.); 

2. Flor (/. wolcottiana, J. Ceja 1565); 3. Tablilla en plano transversal (I. pauciflora, 

J. Ceja 1244); 4. Tablilla en plano tangencial (/. arborescens, J. Ceja 1242); 5. Tablilla 
en plano radial (/. pauciflora, J. Ceja 1244); 6. Distribución de los tejidos que conforman 
los anillos de xilema/floema y posición y forma de los vasos (/. pauciflora, J. Ceja 1244); 
7. Elementos de vaso con punteaduras poligonales (/. rzedowskii, J. Ceja 1242); 

8. Fibras (/. wolcottiana, J. Ceja 1565). XS= Xilema secundario, FS= Floema 

secundario, TC= Tejido conjuntivo. 
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Figuras 9-16. 9. Laticífero en corte radial [Ipomoea popuiina, J. Ceja 1400); 10. Radios 
uniseriados en corte tangencial (/. intrapilosa, JCR 1319); 11. Radios multiseriados, 
corte tangencial (/. pauciflora, JCR 1244); 12. Radios heterocelulares en corte radial 
(/. arborescens, J. Ceja 1232); 13. Granos de almidón en parénquima, (/. pauciflora, 

J. Ceja 1244); 14. Drusas en parénquima {/. murucoides, JCR 1237); 15. Floema en 
corte tangencial {/. popuiina, JCR 1400); 16. Anillos de xilema y floema en corte 
transversal (/. chilopsidis, JCR 1706). L= Laticífero. 




68 


Anatomía de la madera de las especies arbóreas de Ipomoea (Convolvulaceae) 


Tejido conjuntivo. Se le da el nombre de 
tejido conjuntivo a las bandas de parén- 
quima o esclerénquima que separan las 
distintas capas de tejido vascular. En este 
caso está constituido por células paren- 
quimáticas de forma isodiamétrica, no 
lignificadas, que presentan abundantes 
granos de almidón y drusas. 

DISCUSIÓN 

En la mayoría de las plantas con creci- 
miento secundario el cámbium vascular 
se origina entre el xilema y el floema 
primarios, manteniendo esta posición 
indefinidamente, produciendo xilema 
hacia el interior y floema hacia el exterior, 
sin embargo, hay patrones de crecimiento 
secundario que se desvían de este 
modelo, habiéndoseles considerado 
durante mucho tiempo formas de creci- 
miento anómalas o atípicas (Lujan, 1974). 
Actualmente se sabe que estas formas de 
crecimiento se presentan en alrededor de 
20 familias de dicotiledóneas (Metcalfe y 
Chalk, 1950), llegando a ser frecuentes 
en grupos como las Amaranthaceae 
(Rajput, 2001, 2002), Chenopodiaceae 
(Fahn y Zímmermann, 1982), Convolvula- 
ceae (Cariquisty Hanson, 1991), Nyctagi- 
naceae (Cariquist, 2004), Phytolaccaceae 
(Horak, 1981), entre otras. A este tipo de 
crecimiento se le asigna el nombre de 
variante cambial (Cariquist y Hanson, 
1991; Cariquist, 2001). 

El crecimiento secundario en las 
especies arbóreas de Ipomoea es resul- 
tado de la actividad de un tipo de variante 
cambial, en la que distintos cambia suce- 
sivos dan origen a xilema y floema secun- 
darios, además de tejido conjuntivo 
(Metcalfe y Chalk, 1950; Cariquist y 
Hanson, 1991; Cariquist, 2001; 2007). En 
cada uno de los cambia supernumerarios 
el xilema y el floema presentan caracte- 
rísticas y posición típicas, esto es, forman 
anillos concéntricos de bandas de xilema 


y floema secundarios, alternando con 
bandas concéntricas de parénquima. Sin 
embargo, el origen de las regiones meris- 
temáticas es distinto; de acuerdo con 
Luján (1974), para el caso de I. arbores- 
cens e I. pauciflora, cada una de éstas se 
forma a partir de las células parenquimá- 
ticas del floema del último anillo, produ- 
ciendo inicialmente derivados parenqui- 
máticos en ambos sentidos, diferencián- 
dose al final grupos de traqueidas y 
elementos de vaso en el xilema y laticí- 
feros en el floema, lo cual coincide con lo 
observado en los tejidos maduros de las 
especies estudiadas. 

Cariquist (1975) y Cariquist y 
Hanson (1991) han mencionado que las 
lianas con crecimiento secundario típico 
producen de manera característica abun- 
dante parénquima en sus tallos, lo que les 
proporciona una flexibilidad similar a la de 
un cable, facilitando el crecimiento de la 
planta. En las especies estudiadas en 
este trabajo, la variante cambial produce 
abundante parénquima (tejido conjun- 
tivo), entre las bandas de tejido vascular 
de cada anillo, funcionando de manera 
similar a las lianas, lo que apoyaría lo 
postulado por McDonald (1992) en el 
sentido de que en Ipomoea de hábito 
arbóreo, la presencia de cambia supernu- 
merarios es un carácter derivado que 
surge a partir de un ancestro con hábito 
trepador y crecimiento secundario típico, 
aspecto que puede evidenciarse por la 
forma en que se desarrollan las ramas 
jóvenes de algunas de estas especies 
arborescentes. Por su parte Manos et al. 
(2001) también apoyan lo anterior al 
señalar que la forma común en Ipomoea 
es la de una planta con tallo voluble, 
herbáceo o leñoso y que a partir de ella 
han surgido una diversidad de formas 
especializadas. 

Son diversas las implicaciones que 
puede tener el desarrollo de los cambia 
sucesivos en la biología de las plantas 
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que lo presentan, así por ejemplo, se ha 
relacionado la presencia de abundante 
parénquima conjuntivo y radial con la 
necesidad de almacenar agua suficiente 
para sobrevivir a una estación prolon- 
gada de sequía (Cariquist, 2001; OIson, 
2003), en este caso, cabe recordar que 
las especies del grupo en estudio 
habitan principalmente en selvas bajas 
caducifolias, donde se da la condición de 
sequía estacional. También se ha seña- 
lado que ante la escasez de parénquima 
axial, la presencia de grandes radios 
permite almacenar, además de agua, los 
granos de almidón necesarios para el 
sostenimiento de la planta (Rajput y 
Rao, 2003). Lo anterior coincide con lo 
observado en los taxa que integran la 
serie Arborescentes, en los que los 
granos de almidón son abundantes no 
sólo en el parénquima radial, sino 
también en el parénquima conjuntivo, 
constituyendo un reservorio energético 
importante para las distintas etapas 
fenológicas de la planta, como son la 
foliación, floración y fructificación. 

Si bien el estudio anatómico de la 
madera del grupo de Ipomoea arbores- 
centes realizado por McDonald (1992), 
apoyó en la delimitación de las especies y 
en la decisión de McPherson (1981), de 
transferir a Ipomoea carnea ssp. fistulosa 
(con hábito arbustivo y crecimiento 
secundario típico) a la serie Jalapae, 
nuestros resultados y su comparación 
con los obtenidos por Cariquist (2001) y el 
mismo McDonald (1992), muestran que la 
anatomía de las especies dentro del 
complejo estudiado es uniforme en el tipo 
y disposición de los tejidos. Sin embargo, 
cabe señalar que en I. chilopsidis los 
anillos de xilema y floema se disponen en 
cilindros continuos (Figura 16) y el tejido 
conjuntivo no es evidente, a diferencia del 
resto de los taxa, en los que los cilindros 
vasculares se organizan en fascículos 
interrumpidos radialmente por parén- 
quima y el tejido conjuntivo forma de 
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manera alterna bandas tangenciales 
continuas (Figura 6). 

En lo que se refiere a los distintos 
elementos mensurables (Tabla 2), son 
diversos los trabajos en los que se ha 
evaluado la relación que hay entre las 
características anatómico-cuantitativas y 
el tamaño de los individuos, encontrán- 
dose para algunos de ellos una correla- 
ción positiva (Noshiro et al., 1995; 
Terrazas y Loza-Cornejo, 2003; Cariquist, 
2005; Terrazas et al., 2008). En este 
estudio, se puede decir que los análisis 
estadísticos realizados permiten señalar 
que del conjunto de especies estudiadas 
/. chilopsidis e I. populina son las que 
presentan diferencias significativas en un 
mayor número de caracteres (longitud y 
diámetro de las fibras y número de series 
de los radios multiseriados en ambas 
especies; número de células de los radios 
uniseriados y altura de los radios multise- 
riadosen /. chilopsidis y altura de los 
radios uniseriadosen /. populina) con 
respecto a las demás especies, por lo que 
es interesante notar que I. populina es 
una liana e I. chilopsidis es un arbusto; sin 
embargo /. rzedowskii, que también tiene 
hábito arbustivo, si bien presenta simili- 
tudes en algunas características con los 
dos taxa señalados arriba (longitud de las 
fibras o altura de los radios uniseriados), 
en otras queda junto con las especies que 
exhiben un crecimiento arbóreo, por lo 
que consideramos que para establecer si 
existe alguna relación entre el hábito y las 
dimensiones de los caracteres, sería 
necesario ampliar el muestreo y ver si las 
características aquí señaladas perma- 
necen constantes. 

Por último, Cariquist (2001) señala en 
su trabajo que la diversidad anatómica de 
la madera registrada en los géneros que 
integran a la familia Convolvulaceae 
(distintos tipos de crecimiento y diferencias 
en algunos caracteres entre las especies 
arbóreas y lianoides), ofrece la posibilidad 
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de ser correlacionada con los sistemas 
taxonómicos existentes. Sin embargo, en 
este estudio, todas las especies presentan 
la misma variante cambial y a pesar de 
que algunas de las características cuanti- 
tativas consideradas tuvieron diferencias 
significativas, la variación que se presenta 
dentro de cada una de ellas limita su uso 
como carácter taxonómico. 

CONCLUSIÓN 

La serie Arborescentes de Ipomoea, 
presenta una serie de cambia supernumera- 
rios que producen anillos concéntricos de 
xilema, floema y tejido conjuntivo. En 
general las características anatómicas cuali- 
tativas y cuantitativas presentes en cada 
uno de los taxa, coinciden con lo señalado 
en otros trabajos sobre el grupo, sin 
embargo, el análisis estadístico de los carac- 
teres cuantitativos muestra que algunos de 
ellos (longitud y diámetro medio de las 
fibras, altura y número de células de los 
radios uniseriados, altura y número de 
series de los radios multiseriados) presentan 
diferencias significativas que pudieran estar 
relacionadas con el hábito. Estudios especí- 
ficos permitirán contrastar esta hipótesis. 
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lARTlCULO DE INVESTIGACION! 


Modelación espacial de área basal y volumen 
de madera en bosques manejados de 
Pinus patula y P. teocote 
en el ejido Atopixco, Hidalgo 

Spatial modeling of basal area and tree volume in managed 
Pinus patula and P. teocote forests in the ejido Atopixco, 

Hidalgo 

Isidro Ah med Cruz-Leyva'', José René Valdez-Lazalde2, 

Gregorio Ángeles-Pérez2 
y Héctor Manuel de los Santos-Posadas2 

RESUMEN 


Se modeló el área basal (AB m2/ha) y volumen de madera {VOL ms/ha) de rodales coetáneos 
de Pinus patula y P. teocote localizados en Hidalgo, México, utilizando variables de porcentaje de 
cobertura arbórea (COB), índice de área foliar en sitio (LAI5 m2/m2) y alométrico {LA/AL m2/m2), varia- 
bles cartográficas y de reflectancia obtenidas de una imagen multiespectral del satélite SPOT 5. Para 
la caracterización de las condiciones ambientales se digitalizó la carta de Efectos Climáticos (escala 
1; 250 000) de INEGI y el Modelo Digital de Elevación (MDE) correspondiente a la zona. Se integró 
un Sistema de Información Geográfica (SIG) constituido por las siguientes capas de información: 
altitud, pendiente, exposición, temperaturas, precipitaciones, bandas espectrales, y dos grupos de 
índices de vegetación basados en a) la relación espectral de bandas (NDVI - CTVI) y b) la distancia 
espectral entre vegetación y suelo desnudo (PVI - SAVI). Con datos de campo y del SIG se ajustaron 
modelos de regresión lineal múltiple. Para el AB la relación más significativa fue con LAIAL, Tempe- 
ratura Media Anual {TMA), ALTITUD y la Banda 4 de SPOT (P < 0,001, R2 = 0,9662). Para el VOL la 
relación más significativa fue con LAIAL, Temperatura Media Anual (TMA) y ALTITUD (P < 0,001, R2 
= 0,8856). 

PALABRAS CLAVE: 

Inventario forestal, SIG, variables cartográficas, variables espectrales. 


ABSTRACT 


Basal area (AB m2/ha) and tree volume (VOL m^/ha) were modeled for even-aged Pinus patula 
and P. teocote forest stands located in Hidalgo, México. Variables such as tree crown cover (COB), 
fieid recorded Leaf Area Index (LAI5 m2/m2), allometric Leaf Area Index (LAIAL m2/m2); GIS variables 
and reflectance data derived from a multispectral SPOT 5 image were used as independent variables. 
Environmental conditions were characterized by digitizing the INEGI’s chart of climatic effects (scale 
1: 250 000) and Processing the Digital Elevation Model (MDE) of the study area. Ageographic Infor- 
mation System (GIS) was set up with variables such as Elevation, Slope’s angle and aspect, Tempe- 
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rature, Precipitation, SPOT derived 
spectral bands, and two sets of vegeta- 
tion indices based either on a) the ratio 
of spectral bands (NDVI and CTVI) or on 
b) the spectral distance between vegeta- 
tion and bare soil {PVI and SAVI). 
Múltiple linear regression modeis were 
fitted based on fieid and GIS data. AB 
was significantly related to LAIAL, Mean 
Annual Temperatura {TMA), ALTITUD 
and SPOT Band 4 refiectance (P < 

0,001, R2 = 0,9662). VOL was signifi- 
cantly related to LAIAL, Mean Annual 
Temperatura (TMA) and ALTITUD (P < 

0,001, R2 = 0,8856). 

KEY WORDS: 

Forest inventory, GIS, cartographic 
variables, spectral variables. 

INTRODUCCIÓN 

Una de las razones principales que han 
impedido evaluar de manera confiable los 
recursos forestales es el alto costo 
asociado a los métodos tradicionales de 
mapeo, tanto en términos monetarios como 
de tiempo. Por ejemplo, el primer Inventario 
Nacional Forestal de México requirió 24 
años para completarse de 1961 a 1985 
(Sorani y Álvarez, 1996). Afortunadamente, 
la evaluación de los recursos naturales ha 
cambiado radicalmente en la última 
década, gracias a la tecnología de la 
geomática, en particular a los satélites de 
observación terrestre y los Sistemas de 
Información Geográfica (Camarero y Gutié- 
rrez, 1999; Wallerman, 2003). 

A pesar de que la geomática ha mejo- 
rado las metodologías de los inventarios 
forestales, sus resultados son válidos 
únicamente para escalas pequeñas 
(<1 :250 000) (Velázquez et al., 2002). Para 
escalas grandes (>1:25 000) a escala de 
rodal, la precisión alcanzada a la fecha no 
es suficiente para propósitos de manejo 
forestal (Mákelá y Pekkarinen, 2001). De 
aquí que los esfuerzos por modelar el 
comportamiento del volumen de madera 
(Lee et al., 1995), el área basal (González 
et al., 2004) o la geometría de claros a 


escalas grandes (Silbernagel y Moeur, 
2001), se encaminen a comprender los 
efectos de la absorción y refiectancia de las 
estructura arbórea (hojas y ramas) (McRo- 
berts et al., 2002), las limitaciones de la 
resolución espacial y espectral de los 
sensores remotos (Valdez et al., 2006) y a 
estudiar modelos matemáticos que mejor 
describan el comportamiento ecológico de 
las masas forestales (Franco-López et al., 
2001; García-García etal., 2006). En espe- 
cial, la importancia del análisis estadístico 
en el proceso de interpolación de datos 
espaciales y la dependencia del valor de un 
atributo en función de su localización 
geográfica (Miranda-Salas y Condal, 
2003). En este sentido, el objetivo de este 
trabajo es describir dos variables impor- 
tantes del bosque, el área basal (ab) y el 
volumen de madera (vol) del ejido 
Atopixco, por medio de técnicas de mode- 
lado espacial, con el fin de aportar herra- 
mientas que faciliten la administración 
racional de los recursos forestales. 

OBJETIVOS 

Modelar la distribución geográfica del 
área basal y volumen de madera de las 
especies Pinus patula y P. teocote. 

Obtener un inventario forestal del área 
basal y volumen de madera a través de la 
modelación espacial. 


METODOLOGÍA 
Área de estudio 

El ejido Atopixco tiene una superficie de 
1170,97 ha. Se ubica dentro del municipio 
de Zacualtipán, Hidalgo, entre las coorde- 
nadas 20° 37’ 26” N y 98° 37’ 48” W, y 20° 
35’ 20” N y 98° 35’ 23” W. Este municipio 
colinda al norte con Tianguistengo, al sur 
con Metzquititlán, al este con Veracruz, y al 
oeste con Metztitlán y Xochicoatlán (Figura 
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1). Está conformado por pendientes, 
mesetas y cañones; la altitud media es de 
2100 m. El municipio cuenta con tres ríos 
principales: Pánuco, Zoyatia y Miniahuaco, 
además de la laguna Chapultepec. Su 
sustrato geológico data del periodo Ceno- 
zoico (terciario superior) y está conformado 
por rocas ígneas extrusivas del tipo toba 
ácida (Ts - Ta, al norte) y basalto (Ts - B, al 
sur). El sustrato edáfico está representado 
por suelos acrisol órtico (Ao), luvisol 
crómico (Le) y feozem háplico (Hh). En la 
parte norte del ejido el suelo es rico en 
nutrientes; caso contrario es la parte sur 
donde el suelo es delgado y pedregoso 
(emm, 2000). 

En el área predomina el clima 
templado húmedo (C(m)) y templado 
subhúmedo (C(w2)) con temperatura 
media anual entre 12 °C y 18 °C, tempe- 
ratura del mes más frío entre -3 °C y 


18 °C y temperatura del mes más caliente 
bajo 22 °C. La precipitación media anual 
oscila entre 700 mm y 2050 mm, precipi- 
tación del mes más seco menor de 40 
mm; lluvias en verano; la lluvia invernal 
varía entre 5% y 10,2% del total anual 
(2050 mm) (emm, 2000). 

En los bosques manejados del ejido 
Atopixeo se encuentran especies como: 
P'mus patula, P. teocote y en menor abun- 
dancia P. greggii y P. montezumae, estas 
últimas especies son introducidas. En las 
áreas de vegetación natural existen espe- 
cies tales como: Quercus crassipes, Q. 
laurina, Q. excelsa, Q. castanea, Q. ru- 
gosa, Cornus disciflora, Viburnum sp., 
Cleyera theoides, AInus sp., Arbutus xala- 
pensis, Prunus serótina, Symplocos sp., 
Ternstroemia sp., Vaccinium leucanthum, 
Fagus grandifolia subsp. mexicana, entre 
otras (Ángeles-Pérez, 1995). 



Figura 1. Ubicación de ejido Atopixeo, municipio de Zacualtipán de Ángeles, Hidalgo. 
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Modelación espacial de área basal y volumen de madera en bosques manejados 


Adquisición y procesamiento en gabi- 
nete de datos de campo 

Para establecer los límites geográficos 
del área de estudio y de los diferentes 
rodales coetáneos que lo integran, se 
consultó “El plan general de manejo 
forestal del Ejido Atopixco”. Se utilizó 
además un modelo digital de elevación 
(resolución 30 m, inegi). 

Se definieron 21 anualidades que 
corresponden a 38 rodales coetáneos de 
edades que van de 0 a 26 años. Dentro 
de cada anualidad se ubicaron al azar 3 
sitios circulares de muestreo, cada uno 
de 400 m2 de superficie. Se marcó el 
árbol del norte y central, se midió el 
diámetro a la altura de pecho (dap) de 
todos los árboles y la altura total de al 
menos 5 árboles de dap variable. Con un 
GPS Trimble Geoexplorer III se georefe- 
renció cada uno de los 114 sitios muestre- 
ados (pao, 2004). 

Para estimar el volumen total de 
madera de Pinus spp. se empleó la ecua- 
ción de Carril lo-Anzurez et al. (2004), 
misma que fue ajustada para P. patula 
con datos del área de Huayacocotla, Ver., 
región muy cercana al área de estudio, 
donde P. patula es la especie dominante: 

VOL = exp(-9,7688)x^/(P- x ALtJ''''' 


donde; vol = Volumen total de madera 
(m3), DAP = Diámetro a la altura del pecho 
(cm) y ALT = Altura total (m). 

El volumen total, en especies 
hojosas, fue estimado con auxilio de la 
ecuación ajustada por sag (1976): 

VOL = exp(- 9.80m696)xDAP' ^™'’ x 


Inventario forestal 

Con la información obtenida en los 114 
sitios de muestreo, se estimó el inventario 
de los recursos maderables en el ejido 
Atopixco. Para esto se aplicaron tres 
técnicas de inventariado siguiendo la 
formulación expuesta por Schreuder et al. 
(2005): muestreo simple aleatorio (msa), 
muestreo estratificado al azar (mea) y 
muestreo de regresión (mrr). El propósito 
fue obtener una estimación base del área 
basal (ab) y el volumen de madera (vol) 
que sirviese de referencia para comparar 
las estimaciones realizadas con los 
modelos basados en las variables carto- 
gráficas, espectrales y biofísicas. 

Procesamiento de la información 
espectral (spots hrg) 

Se adquirió una imagen spots de la 
Secretaría de Marina a través del Colegio 
de Postgraduados. La escena fue tomada 
el 18 de abril de 2006, sus características 
se reportan en la tabla 1 . 

Corregida la geometría y la radiación 
de la escena (Lillesand y Kiefer, 2001) el 
siguiente paso fue calcular dos tipos de 
índices de vegetación {ivs). El primer tipo 
es conocido como “basados en la 
pendiente” (Slope-Based Vis), el cual 
resulta de la combinación de las bandas 
rojo e infrarrojo cercano. Bajo esta cate- 
goría se obtuvieron los siguientes índices: 
el índice Vegetación de Diferencia 
Normalizada (ndvi) y el índice de Vegeta- 
ción Transformado (corregido) {ctvi) 
(Eastman, 2001). 

El segundo tipo de ivs utilizados se 
denomina “basados en la distancia” 
(Distance-Based vis). Éstos se utilizan en 
ambientes donde la vegetación no es 
densa y es necesario separar la reflec- 
tancia del suelo de la reflectancia de la 
vegetación. 
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Tabla 1 . Parámetros de la escena spots 


Fecha 

2006-04-18 

17:05:10.5 I 

Insti'umento 

HRG2 

Número de bandas 

4 

Indicador de banda 
espectral 

HI1: Verde 

HI2; Rojo 

HI3: IR Cercano 

HI4: IR Medio 

Longitud de onda (pm) 

0,50-0,59 

0,61 - 0,68 

0,78-0,89 

1,58- 1,75 

Número de qanarcia 

7 

7 

5 

4 

Calibración absoluta de 
ganancias (W/m^/sr/pm) 

2,1 39452 

2,853960 

1,738550 

8,225268 

Angulo de orientación 

12,264436 grados 

Ángulo de incidencia 

20, 1 87256 grados 

Áng utos del sol 

Azimut: 11 1,469109 

Bevación: 66,197548 | 

Filas 

7302 

Columnas 

7784 


Para esto se hace una correlación entre 
las bandas roja e infrarroja, utilizando 
como factor de corrección en pendiente 
(,o) a las zonas de suelo desnudo (pvi) o, 
en su caso, transformando las bandas 
roja e infrarrojo en radianza absoluta 
(sAvi) (Eastman, 2001). Bajo esta cate- 
goría se obtuvieron dos índices: el índice 
Perpendicular de Vegetación {pvi) y el 
índice de Vegetación Ajustado a la 
Pendiente del Suelo {savi). 

Generación del sistema 
de información geográfica 

Para caracterizar las condiciones ambien- 
tales del área de estudio, se generó (en el 
paquete idrisi kilimanjaro) un Sistema de 
Información Geográfica (s/g) a partir del 
modelo digital de elevación (mde) de reso- 
lución espacial de 30 m y de las cartas 
Edafológica y Geológica (escala 1 :50000) 
y de Efectos climáticos (escala 1:250000) 
para los periodos de lluvias (mayo- 
octubre) y secas (noviembre-abril) de 
iNEGi. La información de las cartas se digi- 
talizó por medio de una tableta y del 
paquete arcgis. Las cartas digitales obte- 
nidas fueron transformadas a formato 
ráster con resolución semejante a la del 
MDE. A partir del mde se obtuvo la capa de 
elevación sobre el nivel medio del mar 
{altitud) y se construyeron las capas de 


pendiente {pendiente) y exposición {expo- 
sición) del terreno. Las cartas de isolí- 
neas climáticas (isotermas e isoyetas) se 
interpolaron mediante el módulo intercon 
del paquete idrisi kilimanjaro. De esta 
manera se obtuvieron las cartas de 
temperatura media anual {tma), tempera- 
tura máxima {tmax), temperatura mínima 
{tmin), precipitación media anual {pma), 
precipitación en lluvias {pmamo) y precipi- 
tación en secas {pmaná). Se consideraron 
además las variables biofísicas índice de 
área foliar en sitio (la/5), índice de área 
foliar alométrico {laial) y cobertura 
arbórea (coe) (Aguirre, 2007). 

Predicción de ab y vol mediante 
regresión lineal múltiple (Stepwise) 

Los datos de campo (inventario) corres- 
pondientes a área basal (ab) y volumen 
de madera {vol) se empataron con las 
nueve variables cartográficas del medio 
físico (información ambiental), con los 
datos espectrales obtenidos de la escena 
SPOT5 (Bandas 1-4), con los datos espec- 
trales transformados en índices de vege- 
tación -IV 's {ndvi, ctvi, pvi y saw) y con 
los índices de área foliar. Esta base de 
datos permitió ajustar, por separado, para 
las variables ab y vol, modelos de regre- 
sión lineal múltiple {rlm) de la forma 
general: 
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AB,VOL = i3o + iS,X, +¿ 32 X 2 + ... + + + ... + + ¡3-^ IV + f,. 

AB,VOL = )8o + p^COB + +p^X^+ ... + ¿3,^ X, + p^^B, + ...+ p^^B^ + i3¡ IV + e. 

AB, YOL = Po + P^LA15 + P,X^+ p,X^ + ...+ P^^X^ + P^^B^ +... + P^,B^ + P^ IV + £. 
AB,VOL = j8o + /3,I^/JL + j 82 X, + p^X^ + ... + P^^X^ +Pn^i + ... + P^^B ^ + p. IV + £,. 


Los modelos genéricos arriba descritos 
se ajustaron mediante el procedimiento 
STEPWiSE (Pindyck y Rubinfeid, 1981) 
para predecir los valores promedio de ab 
y voL que se esperaría alcance un sitio 
(píxel) donde las condiciones ambientales 
sean Xi, X 2 ,..., Xg; las propiedades 
espectrales sean Bi, B 4 ; y los IV's 
sean los indicados, al igual que las varia- 
bles biofísicas consideradas. Con el fin de 
identificar y corregir inconsistencias esta- 
dísticas en los modelos, además de las 
pruebas anova {P < 0,05) y del factor de 
inflación de varianza {VIF < 1 o 5), se 
buscaron observaciones aberrantes (O/A) 
con relación al conjunto de datos a 
predecir (Y,) y del conjunto de variables 
predictoras (X¡) (Belsley et al., 2005). 

Finalmente, los cuatro modelos que 
presentaron mejor ajuste, tanto para ab 
como para vol, se usaron para predecir la 
distribución espacial de tales variables en 
la totalidad de la superficie del ejido 
Atopixco. Para ello se resolvieron los 
modelos en cada uno de los pixeles que 
conforman el ejido. 


RESULTADOS Y DISCUSIÓN 

A través del inventario se obtuvo un valor 
promedio de área basal (ab) de 19,96 m2/ha 
bajo muestreo simple aleatorio (msa) y de 
22,54 m2/ha bajo muestreo aleatorio estrati- 
ficado (mea). El volumen de madera (vol) 
promedio obtenido fue de 118,89 m^/ha 
(msa) y de 140,69 m^/ha para las 21 anuali- 
dades o estratos (mea). Se estimó un inven- 
tario final de 42 019,11 m3, para el Ejido 
Atopixco, con una variación estimada en vol 
de ± 2 218,59 m3 (5,28%) (Tabla 2). 

Para mejorar la precisión porcentual (B/y) del 
inventario en vol, se procedió a aplicar los 
estimadores de muestreo de regresión (mrr). 
Si se utiliza solamente como valor de Px al 
promedio de los datos, el promedio estimado 
es igual al valor de la media, obtenida bajo 
MSA (118,89 m3/ha). Sin embargo, se tiene un 
estimador alterno de Px que es el valor de ab 
promedio obtenido bajo mea y que es mucho 
más preciso que el obtenido bajo msa. 
Usando este valor como Px para el estimador 
de regresión se tiene que la media bajo mrr 
es Yreg = 5,55 (138,8 m3/ha) (Tabla 3). 


Tabla 2. Estimadores por sitio de 400 m 2 en muestreo simple aleatorio (msa) 

y estratificado aleatorio (mea). 


ESTIMADORES 

ÁREA BASAL (m2) 

VOLUMEN (m3) 






Media 

0,7985 

0,90 

4,76 

5,63 

Varianza media estimada 

0,0021 

0,00033 

0,13 

0,02 

Límite superior de IC 

0,89 

0,94 

5,47 

5,92 

Límite inferior de IC 

0,71 

0,86 

4,04 

5,33 

Error permisible porcentual (B/y) 

11 , 53 % 

4 , 06 % 

14 , 94 % 

5 , 28 % 

Amplitud del IC a 95% (B) 

0,09 

0,04 

0,71 

0,30 

Inventario total (T) 

5961,76 

6730,69 

35507,37 

42019,11 

Media Alfa 

0,04 

0,05 

0,24 

0,28 

Tamaño de muestra, -(n) 

568,27 

65,07 

908,02 

100,78 

Alfa (a) 

0,05 

0,05 

0,05 

0,05 
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Tabla 3. Estimadores de volumen por 
sitio de 400 m2, en muestreo 
por regresión (mrr). 


ESTIMADORES DE REGRESIÓN 

VOLUMEN 

(m3) 

b (pendiente) 

7,38 

Correlación VOL vs. AB 

0,96 

Media 

5,55 

Covarianza de XY (Syx) 

1,29 

Varianza media regresión (SmedReg) 

0,01 

Límite superior de IC 

4,97 

Límite inferior de IC 

4,54 

Error permisible porcentual (B/y) 

3 , 80 % 

Amplitud del IC a 95% (B) 

0,21 

Inventario total (T) 

41454 


Este resultado sugiere que asumir msa 
tiende a subestimar considerablemente 
las existencias totales en volumen de 
madera y es mucho menos preciso, 
siendo los estimadores de mea y mrr más 
adecuados para estimar el inventario. 

Predicción de ab mediante regresión 
lineal múltiple 

La tabla 4 presenta las variables significa- 
tivas y los valores de los estimadores que 
describen los cuatro modelos de regre- 
sión lineal múltiple ajustados para 
predecir el ab mediante el procedimiento 
Stepwise. El coeficiente de determinación 
múltiple ( R 2 ) de los modelos propuestos 
es moderadamente alto, sus valores 
varían explicando de 73,1% a 96,62% de 
la variabilidad. 

Cabe resaltar que el modelo 3 es el 
que presenta la mayor capacidad predic- 
tiva ( R 2 = 0,96) mediante la inclusión de la 
variable biofísica laial y las variables 
cartográficas altitud y tma. Los modelos 
1 a 4 (Tabla 4) se usaron para predecir los 
valores de ab a escala de píxel. 


Modelado espacial del área basal [ab) 

El modelado espacial se obtuvo resolviendo 
los modelos para la totalidad de los pixeles 
que conforman el área de estudio, constru- 
yendo así la distribución espacial del ab en el 
ejido Atopixco. Por ejemplo, la figura 2 
muestra tal distribución con base en el 
modelo etiquetado como Modelo 1 , donde la 
relación entre ab y coe es positiva y alta 
(64,9% correlación semiparcial). La coe es 
un reflejo de la densidad. Márquez-Linares y 
Álvarez-Zagoya (1995) muestran que en 
Pinus cooperi de diferentes clases de copa, 
hay una relación negativa entre el diámetro 
y la densidad conforme la competencia 
aumenta. Es decir, la competencia tiene un 
efecto importante en la pendiente de la rela- 
ción tamaño (diámetro)-densidad, a partir 
del punto en el cual la competencia se hace 
más intensa (individuos intermedios). 

También la alta densidad, en rodales 
jóvenes, puede producir árboles con bajo 
porcentaje de ramas vivas a lo largo del 
tronco, coberturas del rodal muy altas y 
fustes cónicos (Black y Abrams, 2005). 
Esto puede ser captado por las fotos 
hemisféricas (Jonckheere et al., 2004) y 
explicaría porqué rodales jóvenes 
presentan altos valores de cobertura 
arbórea (Figura 2). 

En la figura 3 se muestra la distribu- 
ción espacial del ab pero ahora con base 
en el modelo 2 (Tabla 4), aquí la relación 
más alta del ab se da con el índice de 
área foliar (lais), misma que es positiva 
(68% correlación semiparcial). 

En los cuatro modelos ajustados 
para predecir ab, destaca la variable 
índice de área foliar alométrico (laial) 
(Figura 4), es ésta la que más aportó a la 
R 2 ajustada (Tabla 4). A diferencia de lais, 
el LAIAL está altamente relacionado con el 
área basal, pues su estimación depende 
del diámetro normal (dap) (Cano-Morales 
et al., 1996; Figueroa-Navarro, 2010). 
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Tabla 4. 

E 

R-cuadrada 

R-ajustada 

Intercepto 

COB 

LAI5 

LAIAL 

PENDIENTE 

ALTITUD 

ALTITUD2 

TTVIA 

TMAX 

TMAX2 

PMA 

PMAMO 

PMAM02 

PVI 

B4 


R-cuadrada 

R-ajustada 


Intercepto 

COB 

LAI5 

LAIAL 

SLOPE 

ALTITUD 

ALTITUD2 

TMA 

TMAX 

TMAX2 

PMA 

PMAMO 

PMAM02 

PVI 

B4 


Significancia {P< 0,001) del modelo de regresión lineal múltiple para estimar 
área basal (AB) y estimadores de los parámetros |3i. 


Modelo 1 


Modelo 2 



0,777 




0,768 



Estimador 

parámetro 

Pr>| t| 

Corr., tipo 1 
sem i parcial 
cuadrado 

FIV 

8,222 

0,0001 


0 

0,016 

0,0001 

0,64984 

1,67 

-0,018 

0,0029 

0,02113 

1,1 

-0,565 

0,0001 

0,09731 

1,68 

-0,054 

0,0468 

0,00912 

1,01 

Modelo 3 


0,9662 




0,9649 



Estimador 

parámetro 

i 

Pr>| t| 

Corr., tipo 1 
sem i parcial 
cuadrado 

FIV 

6,554 

0,0001 


0 

0,138 

0,0001 

0,947 

1,9 

-0,001 

0,009 

0,003 

1,86 

-0,27 

0,0001 

0,014 

2,77 

-0,002 

0,0306 

0,002 

1,11 



0,8058 




0,7938 



Estimador 

parámetro 

Pr> |t 1 

Corr., tipo 1 
semi parcial 
cuadrado 

FIV 

299,52 

0,004 


0 

0,39 

0,0001 

0,68 

1,82 

-0,01 

0,0296 

0,01 

1,24 

-0,29 

0,0051 

0,0001 

>5 

0,0001 

0,0052 

0,03 

>5 

■0,49 

0,0001 

0,08 

3,2 

-0,05 

0,0453 

0,01 

1,02 

Modelo 4 

0,7313 

0,7117 

Estimador 

parámetro 

Pr> |t 1 

Corr., tipo 1 
semi parcial 
cuadrado 

FIV 

1177,91 

0,0001 


0 

-0,576 

0,0001 

0,18 

>5 

0,0001 

0,0001 

0,07 

>5 

-58,827 

0,0001 

0,05 

>5 

1,656 

0,0001 

0,001 

>5 

0,183 

0,0001 

0,29 

>5 

-0,17 

0,0001 

0,001 

>5 

0,00002 

0,0001 

0,13 

>5 


COB = Cobertura arbórea, LAI5 = índice de área foliar en sitio, LAIAL = índice de área foliar alométrico, PENDIENTE = Pendiente del 
terreno, ALTITUD = Altitud sobre el nivel medio del mar, TMA = Temperatura media anual, TMAX = Temperatura máxima, PMA = Preci- 
pitación media anual, PMAMO = Precipitación media anual Mayo-Octubre, PVI = índice perpendicular de vegetación, B4 = Banda 4. 






sfnm 


Madera y Bosques 16(3), 2010:75-97 


83 



Area basaí. 

(m2/ha) 


[ooo- 
[ose- 
í 12<^- 
193 
25S- 
3 .. 1 ? - 

r~1 3 rs... 

^437. 
r ■; S04 .. 
|$7f. 
[ 633 - 

leas- 

753 . 
[ 610 - 
[ S6S- 
[925- 
I 992- 
i 10 79 


057 
125 
1 92 
2^ 
3ie 

3 74 

4 36 
S® 
570 
633 
6 94 
752 
Bm 
8 67 
924 
991 

m ?6 
- 12 te 


O 


9 ns frA 


«A 


*4 


Figura 2. Distribución del área basal (m^/ha) en el ejido Atopixco, Hidalgo 

con base en el Modelo 1 (Tabla 4). 
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Figura 3. Distribución del área basal (m2/ha) en el ejido Atopixco, Hidalgo 

con base en el Modelo 2 (Tabla 4). 
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La relación entre el ab y laial puede expli- 
carse como una extensión de la relación 
que existe entre el tejido vascular (xilema 
y floema) y la superficie foliar. Según 
Barrantes y Gracia (1989), el área foliar 
específica de una hoja (medida en 
cm2/mg) varía con la altura del árbol. Las 
hojas situadas a mayor altura (más 
expuestas a la radiación solar) presentan 
un área específica menor que las 
halladas en la parte baja de la copa. En 
otras palabras, el área foliar es propor- 
cional a la superficie transversal del 
tronco, por el cual pasa el agua necesaria 
para abastecerlo y que circula desde las 
raíces hacia las hojas. Este es un fenó- 
meno muy importante que tiene varias 
implicaciones; pues no solo relaciona la 
cantidad de área basal con el índice de 
área foliar, también muestra que la dispo- 
nibilidad de agua influye en la estructura 
arbórea. 


En todos los modelos, las variables 
cartográficas altitud, tma, pma y pmamo 
siempre se hacen presentes. Mención 
aparte merecen las variables de reflec- 
tancia en la banda B4 y el pvi. En el 
Modelo 3 (Tabla 4), el coeficiente de la 
variable B4 es negativo, esto significa que 
la cantidad de área basal aumenta 
conforme disminuye la reflectancia en el 
infrarrojo cercano (b4). En un estudio con 
plantaciones de Eucalyptus y Pinus, 
Lencinasa y Mohr-Bell (2007) determi- 
naron una correlación significativa (-0,58 
y -0,79) entre los valores de la b4 y la 
edad de dichas plantaciones. Estimaron 
la edad con un error medio de 2,2 años en 
eucaliptos y 3,5 años en pinos. Esa rela- 
ción inversa, entre valores espectrales y 
la edad para plantaciones de coniferas, 
también fue reportada por Sivanpillai et 
al. (2006). Tal fenómeno se explica por el 
agua de las hojas que absorbe la luz roja, 
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Figura 4. Distribución del área basal (m2/ha) en el ejido Atopixco, Hidalgo 

con base en el Modelo 3 (Tabla 4). 
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mientras que las demás estructuras celu- 
lares reflejan la mayor parte de la radia- 
ción infrarroja (s4). El incremento en área 
basal está acompañado por un aumento 
en el índice de área foliar alométrico 
(laial) y por tanto una mayor absorción 
del infrarrojo medio (b 4), razón que 
expone la relación negativa de AB-B4 
(Figura 4, Tabla 4). 

De acuerdo con Ante et al. (2005), 
en general los índices de vegetación que 
toman la razón entre la reflectancia del 
suelo y la vegetación (por ejemplo savi, 
pvi, entre otros) son mejores para discri- 
minar tipos de cobertura vegetal, en 
zonas con contrastes ambientales, que 
los índices de vegetación basados en la 
razón de bandas rojo-infrarrojo (ndvi, 
CTVi, entre otros). Es de llamar la atención 
que el coeficiente de pvi (Modelo 2) sea 
negativo, en tal caso habría que explorar 
su relación con la variable lai5 (Tabla 4). 

De lo anterior, Gilabert et al. (2002) 
escriben “el pvi se puede citar como un 
índice que presenta, dentro del rango de 
variación, una dependencia lineal con 
respecto al índice de área foliar (lai), y al 
NDVi, como aquél en el que se observa 
más claramente la saturación”. Para 
valores intermedios -en densidad de vege- 
tación- la dispersión de la radiación por 
parte de la masa vegetal puede disminuir 
la influencia del suelo, en este caso los 
índices de vegetación tradicionales (como 
el ndvi) presentan los valores más altos, 
indicando una eficacia menor para norma- 
lizar la perturbación introducida por la 
señal del suelo. Se manifiesta la limitación 
de estos índices en zonas que presentan 
vegetación dispersa (es el caso del ejido 
Atopixco). El PVI parece funcionar mejor, 
probablemente porque define con suficiente 
precisión la línea del suelo (Figura 3). 

Si bien las variables cartográficas 
ALTITUD, TMA, PMA y PMAMO S6 hacen 
presentes en los modelos 1, 2 y 3, es en 
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el Modelo 4 (Figura 5), donde cobran rele- 
vancia. Como su nombre lo indica, son 
variables obtenidas de información carto- 
gráfica, lo cual puede ser una virtud o una 
desventaja (Trotter y Dymond, 1997). En 
este sentido Miranda-Salas y Condal 
(2003) afirman que los atributos altitud y 
área basal (ab) satisfacen la condición de 
variables interpoladles (modelado espa- 
cial) y que el proceso de interpolación no 
puede ser exitoso si primero no se consi- 
deran los resultados obtenidos de un 
análisis estadístico preliminar de los 
datos a interpolar. Entonces se debe ser 
cuidadoso al considerar variables carto- 
gráficas en la modelación espacial. 

Como se aprecia en la figura 5 las 
variables cartográficas (Tabla 4) 
describen un rango de valores para ab 
que van de 12,31 m2/ha a 16,09 m2/ha. 
Es un rango estrecho en comparación 
con los modelos que incluyen a las varia- 
bles biofísicas COB e índices de área foliar 
{LAI5 y laial). Además se hace evidente el 
peso que tiene la precipitación y la tempe- 
ratura sobre la modelación espacial del ab 
(Figura 5). Esto indicaría que los modelos 
que sólo contemplan variables cartográ- 
ficas son menos eficaces para la modela- 
ción espacial de variables biofísicas (coe 
y lai) o dasométricas (ab) (Zurita-Milla, 
2000). Pero no es así, el mayor problema 
de utilizar tales variables es su construc- 
ción. Por ejemplo, las capas ráster de 
temperatura y precipitación fueron obte- 
nidas de cartas escala 1:250 000. Con 
una escala tan pequeña es complicado 
realizar operaciones de interpolación, aun 
cuando se utilice una base ráster óe 30 m 
X 30 m de resolución, el resultado es poca 
variabilidad espacial. Aún con ésta y otras 
desventajas se pudieron obtener coefi- 
cientes de correlación moderadamente 
altos (Tabla 4). De mejorar la calidad de la 
información contenida en mapas se daría 
un enorme avance en la modelación de 
variables dasométricas (Lillesand y 
Kiefer, 2001; Eastman, 2001). 
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Figura 5. Distribución del área basal (m2/ha) en el ejido Atopixco, Hidalgo 

con base en el Modelo 4 (Tabla 4). 


Sin olvidar lo mencionado arriba, queda 
de manifiesto la utilidad de emplear varia- 
bles cartográficas como altitud, precipita- 
ción y temperatura para describir la distri- 
bución del AB. Los coeficientes de altitud 
y temperatura guardan una relación posi- 
tiva con el incremento de ab. Esto 
concuerda con el trabajo de Camarero y 
Gutiérrez (1999) quienes describieron el 
patrón espacial de un ecotono de bosque 
dominado por Pinus uncinata, a través de 
variaciones en altitud, tipo de sustrato, y 
cobertura de herbáceas y arbustos. La 
descripción del ecotono permitió discernir 
la posible respuesta de las poblaciones 
de R uncinata ante cambios ambientales. 
Por su parte, Sánchez y López-Mata 
(2003) al realizar un análisis sobre la 
vegetación y los factores que influyen en 
su estructura y distribución a lo largo de la 
Sierra Nevada, notaron que la relación 
(correlación canónica) entre vegetación y 


altitud-precipitación fue positiva hasta un 
cierto rango, mientras que la temperatura 
presenta coeficientes canónicos nega- 
tivos. La temperatura en esa zona 
templada se asocia a sitios de cobertura 
vegetal seca y menor cantidad de precipi- 
tación. La misma tendencia fue obser- 
vada en el ejido Atopixco. Un hecho 
importante destacado por Sánchez y 
López-Mata (2003) fue que la distribución 
de las especies era continua a lo largo del 
gradiente. En consecuencia el patrón en 
la estructura y distribución de las comuni- 
dades arbóreas está definido por el 
gradiente altitudinal (temperatura y preci- 
pitación) y por las propiedades del suelo 
(pendiente del terreno, materia orgánica, 
cationes, nutrimentos y la profundidad del 
horizonte A y O). Elementos que podrían 
valorarse en futuras investigaciones para 
el ejido Atopixco. 
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Comparación de las estimaciones de 
área basal: inventario vs modelación 

De los cuatro modelos propuestos para 
estimar el área basal (ab), el Modelo 2 
(que incluye lais) se aproximó más a la 
estimación del inventario total. El valor de 
AB obtenido para el ejido por mea fue de 
22,54 m2/ha, mientras que la estimación 
mediante el modelo 2 fue de 20,47 m2/ha, 
es decir, con una diferencia en área basal 
de 2 m2/ha; un margen menor si se 
compara con la estimación mediante msa 
(2,6 m2/ha). Estos resultados concuerdan 
con los presentados por Aguirre (2007) 
quien observó que el ab fue el parámetro 
forestal que registró el error cuadrático 
medio (recm) más pequeño para pinos. 
Un dato muy importante, pues el área 
basal siempre fue la variable respuesta 
que mejor ajuste tuvo en relación con las 
variables independientes. Por otro lado, 
destaca el resultado obtenido con el 
Modelo 4, el cual está integrado única- 
mente por variables cartográficas. A pesar 
de tener la menor R2 ajustada (0,7117) 
fue el segundo modelo que arrojó estima- 
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dones más cercanas al inventario por 
MSA (tablas 4 y 5). 

Estimación de vol mediante regresión 
lineal múltiple 

La tabla 6 presenta las variables significa- 
tivas y los valores de los estimadores que 
describen los modelos de regresión lineal 
múltiple ajustados para predecir el vol. 
Estos modelos fueron la base para 
predecir la distribución espacial del vol 
en la superficie forestal del ejido Atopixco. 
Los cuatro modelos obtenidos (modelos 
5-8) ponen de manifiesto la importancia 
de utilizar variables biofísicas como coe, 
LAI5 y laial; y cartográficas como tma, 
PENDIENTE, ALTITUD, PMA y PMAMO en la 

estimación del vol. Merece la pena 
destacar también los valores obtenidos 
de los coeficientes de determinación 
múltiple (R2), mismos que oscilan entre 
0,72 y 0,88; este último se considera alto, 
representando una capacidad explicativa 
del 88%. 


Tabla 5. Comparación de las estimaciones de área basal (m2) con inventario tradicional 
y con inventario por siG. Intervalos de confianza (ic) a 95% de significancia estadística. 


Superficie (ha) 298,66 

Área Basal 

Límite Límite _ 

2 „ ■ í ■ Precision 

m /ha TOTAL inferior superior o/ /□/ \ 

delIC delIC 

Inventario 

Tradicional 

Inventario MSA 
Inventario MEA 

19,96 

22,54 

5961,76 

6730,69 

5274,66 

6457,56 

6648,85 

7003,83 

11,53% 

4,06% 

Inventario SIG 
SIG 

Modelo 1 (COB) 
Modelo 2 (LAI5) 
Modelo 3 (LAIAL) 

Modelo 4 
(Cartográficas) 

12,16 

20,47 

13,35 

16,09 

3631,76 

6113,67 

3987,18 

4805,52 



COB = Cobertura arbórea MSA = Muestreo Simple Aleatorio 

MEA = Muestreo Estratificado 
LAI5 = índice de área foliaren sitio Aleatorio 

LAIAL = índice de área foliar alométrico 
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Tabla 6. Significancia (P< 0,001) del modelo de regresión lineal múltiple para estimar 

voL y estimadores de los parámetros |3i. 


Modelo 5 


Modelo 6 


R-cuadrada 

R-ajustada 


0,7880 

0,7794 




Estimador 
parámetro [i 

Pr > 1 1 1 

Corr., tipo I 
semiparcial 
cuadrado 

FIV 

Intercepto 

1003,38912 

0,0001 


0 

COB 

IAI5 

lAIAL 

0,093© 

0,0001 

0,5639 

1,67 

PENDIENTE 

ALTITUD 

ALTITUD2 

-0,11661 

0,0112 

0,0143 

1,16 

IMA 

-133,53265 

0,0001 

0,16547 

>5 

1MA2 

IMAX 

1MAX2 

PMA 

PMAMO 

PMAM02 

4,4283 

0,0002 

0,04434 

>5 


Modelo 7 


R-cuadrada 

R-ajustada 


0,8856 

0,881 




Estimador 
parámetro |i 

Pr > 1 1 1 

Corr,, tipo I 
semiparcial 
cuadrado 

FIV 

Intercepto 

COB 

IAI5 

1048,44 

0,0001 


0 

lAIAL 

PENDIENTE 

0,77 

0,0001 

0,773 

1,86 

ALTITUD 

ALTITUD2 

-0,01 

0,0417 

0,005 

1,83 

IMA 

-137,06 

0,0001 

0,068 

>5 

1MA2 

IMAX 

1MAX2 

PMA 

PMAMO 

PMAM02 

4,59 

0,0001 

0,04 

>5 


0,8137 

0,8062 

Estimador 
parámetro |i 

Pr > 1 1 1 

Corr., tipo I 
semiparcial 
cuadrado 

FIV 

705,42 

0,003 


0 

2,56 

0,0001 

0,65 

1,82 

-0,13 

0,0029 

0,02 

1,17 

-92,36 

0,0049 

0,12 

>5 

3,02 

0,0076 

0,02 

>5 


Modelo 8 



0,7404 

0,7215 



Estimador 
parámetro |ii 

Pr > 1 1 1 

Corr,, tipo I 
semiparcial 
cuadrado 

FIV 

8093,86 

0,0001 


0 

-3,99 

0,0001 

0,18 

>5 

0,001 

0,0001 

0,06 

>5 

-393,47 

0,0004 

0,02 

>5 

11,11 

0,0004 

0,0001 

>5 

1,49 

0,0001 

0,33 

>5 

-1,39 

0,0001 

0,001 

>5 

0,0001 

0,0001 

0,15 

>5 


COB = Cobertura arbórea, LAI5 = índice de área foliar en sitio, LAIAL = índice de área foliar alométrico, PENDIENTE = Pendiente del 
terreno, ALTITUD = Altitud sobre el nivel medio del mar, TMA = Temperatura media anual, TMAX = Temperatura máxima, PMA = Preci- 
pitación media anual, PMAMO = Precipitación media anual Mayo-Octubre. 
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Modelado espacial del volumen de 
madera {vol) 

En la figura 6 se muestra la distribución 
del VOL en el ejido Atopixco con base en 
el Modelo 5 (Tabla 6). Es de destacar la 
importancia de las variables cob y tma en 
la distribución. 

La temperatura se correlaciona 
negativamente con el vol, es decir, a 
mayor temperatura menor cantidad de 
madera. Regionalmente esto puede 
asociarse a la disponibilidad de humedad 
(precipitación) y cómo ésta influye en la 
actividad del árbol y de algunos microor- 
ganismos. Como menciona Liétor- 
Gallego (2002), en bosques de Abies 
pinsapo, las condiciones extremas de 
temperatura y humedad del suelo 
retardan la actividad de los microorga- 
nismos relacionados con el ciclado del N 
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(mineralización y nitrificación). Al final del 
invierno y principios del otoño, coinci- 
diendo con los máximos de producción de 
hojarasca y con la recuperación de condi- 
ciones de temperatura y humedad favora- 
bles, se incrementan los niveles de N 
mineral en el suelo y los microorganismos 
son más efectivos para liberar nutrientes. 
Un fenómeno que podría explicar en 
parte la distribución arbórea en el ejido 
Atopixco. 

La COB es la variable que más aportó 
a explicar el comportamiento del vol 
(Tabla 6). Esta variable es sencilla de 
medir en campo, sin embargo, al 
momento de modelar se observa que la 
PENDIENTE influye sobre la eos y por ende 
sobre el vol, aún cuando su coeficiente 
de correlación semiparcial sea bajo. 
Como explica Valdez-Lazalde et al. 
(2006) una opción adecuada para estimar 
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Figura 6. Distribución del volumen de madera (m^/ha) en el ejido Atopixco, Hidalgo 

con base al Modelo 5 (Tabla 6). 
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la cobertura de copas, cuando las áreas 
tienen condiciones fisiográficas acciden- 
tadas, es contar las sombras registradas 
(en imágenes de satélite) como áreas de 
vegetación; contrariamente, en sitios de 
poca inclinación las sombras efectiva- 
mente son suelo sin cubierta vegetal. 
Esto de forma indirecta puede revelar 
porqué la pendiente influyó sobre el vol, 
las sombras formadas por las pendientes 
debieron contarse como coe (Figura 6). 

En el Modelo 6 (Figura 7) la 
pendiente también produce un efecto 
positivo sobre la cantidad de vol. Sin 
embargo, su contribución al modelo 
(Tabla 6) no se compara con la aportación 
del LAI5. Tanto la coe como el lai5 de- 
penden de la geometría de copas, de la 
arquitectura de hojas y ramas, así como 
de la cantidad de luz que éstas dejan 
pasar hasta la cámara (Frazer et al., 


1999), no obstante hay una mejor relación 
lineal del vol y lais que del vol y coe, 
razón por la cual el coeficiente de deter- 
minación múltiple del Modelo 6 es mayor 
que el coeficiente del Modelo 5 (Tabla 6). 

Al igual que en el Modelo 5, también 
en el Modelo 6 (Figura 7) se hace presente 
la TMA, misma que puede asociarse con la 
distribución de la humedad. De hecho, 
todos los modelos muestran una tendencia 
hacia el aumento creciente del vol, desde 
el suroeste hacia el noreste; lo cual coin- 
cide con el cambio del régimen climático de 
una zona árida (Valle de Metztitlán) hacia 
una zona templada (ejido La Mojonera). 
Esto viene apoyado por la presencia de la 
variable altitud que representa el gra- 
diente altitudinal (Figura 7). 

De los modelos ajustados para 
estimar vol, el Modelo 7 (Figura 8) es el 
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Figura 7. Distribución del volumen de madera (m^/ha) en el ejido Atopixco, Hidalgo 

con base al Modelo 6 (Tabla 6). 
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que presenta el mayor coeficiente de 
determinación (Tabla 6). En este caso, es 
evidente la alta relación existente entre la 
variable laial y el vol a través del ab. En 
este modelo, el gradiente climático 
(temperatura y elevación sobre el nivel 
del mar) ejerce un marcado efecto sobre 
la distribución del vol. Se esperaría que 
el modelo espacial que incluye al laial 
fuera mejor que los modelos que incluyen 
como variables dependientes a la cob y al 
LAI5, dado que su R2 fue la más alta de 
todos los modelos. No obstante, en 
términos generales, el Modelo 7 vuelve a 
mostrar la tendencia de aumento en vol 
desde el suroeste hacia el noreste (ejido 
La Mojonera). En este sentido, es posible 
perfeccionar las predicciones en vol si 
por un lado se aumenta la precisión de las 
variables cartográficas y por otro se halla 
un modelo matemático que emule la 
complejidad del paisaje. Esto es apoyado 


91 

por García-García et al. (2006) y McRo- 
berts et al. (2002) quienes comprobaron 
que los modelos espaciales son suscepti- 
bles de mejoras si se reduce la heteroge- 
neidad del bioma, lo cual hace que la rela- 
ción entre variables forestales y variables 
radiométricas no sea tan alta como las 
observadas en otros medios forestales 
más homogéneos. 

De los cuatro modelos ajustados 
para predecir el vol, el Modelo 8 (Tabla 6) 
generado a partir de variables cartográ- 
ficas presenta el menor coeficiente de 
correlación (R2 = 0,74). Esto se ve refle- 
jado en el intervalo del vol que va de 99,9 
m3/ha a 126,6 m^/ha (Figura 9), en 
comparación al Modelo 6 tiene un inter- 
valo de 0,6 m3/ha a 213 m3/ha. Al igual 
que en el modelo de área basal (Figura 
5), los coeficientes para altitud, tmax, 
PMA y PMAMO indican que la cantidad de 
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Figura 8. Distribución del volumen de madera (m3/ha) en el ejido Atopixco, Hidalgo 

con base en el Modelo 7 (Tabla 6). 
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VOL incrementa con la altitud y con el 
aumento en precipitación (pma, pmamo). 
Por el contrario, la cantidad de vol dismi- 
nuye al haber un aumento en la tempera- 
tura {TMAX), lo que concuerda con lo 
observado en los modelos 5, 6 y 7. 

Finalmente, el mayor reto al mo- 
delar el VOL fue la escala de las variables 
cartográficas (<1:250 000) de la zona de 
estudio. Esto produjo extensas áreas 
donde las variables {altitud, pma, pmamo, 
PENDIENTE, TMA y TMAx) no mostraban 
variabilidad espacial. Como describen 
Fazakas etal. (1999), Franco-López etal. 
(2001) y Mákelá y Pekkarinen (2001), la 
estimación del volumen de madera (con 
sensores remotos) puede aplicarse a 
zonas pequeñas, siempre y cuando se 


tenga presente que esa estimación 
depende dei cuadrado medio del error 
{rmse) de lo estimado y lo medido en 
campo. La pregunta relevante es ¿cómo 
reducir el error medio a la vez que se 
aumenta la varianza explicada? Al 
parecer las variables cartográficas son un 
tercio del camino, los otros dos tercios 
son representados por las variables biofí- 
sicas (coe, LAI5, laial). Al respecto Lee et 
al. (1995) mostraron que es posible ob- 
tener coeficientes de correlación múltiple 
elevados (R2 > 0,90) si se posee informa- 
ción geográfica de escala adecuada. En 
este caso, se tenían variables categóricas 
de especies, edades y cobertura para 
toda la zona de estudio. El secreto radicó 
en la calidad de los mapas. 
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Figura 9. Distribución del volumen de madera (m^/ha) en el ejido Atopixco, Hidalgo 

con base en el Modelo 8 (Tabla 6). 
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Comparación de las estimaciones de 
volumen de madera: inventario vs 
modelación 

Todos los modelos propuestos para 
predecir y describir la distribución espa- 
cial del volumen de madera (vol), 
presentaron sobreestimaciones con rela- 
ción al inventario obtenido mediante mrr 
(Tabla 7). Sin desatender este punto, fue 
el Modelo 8 (que considera variables 
cartográficas) el que mejor se ajustó a la 
estimación de inventario total. Si se toma 
como referencia el inventario obtenido 
por MRR (41453,16 m^), se observa que 
tal estimación es similar al obtenido por 
el Modelo 8 (37804,99 m^). La diferencia 
entre ambas estimaciones es de 
3648,16 m3. Por el contrario, si se toma 
como referencia al inventario obtenido 
mediante mea, la mayoría de los modelos 
(5 al 7) muestran sobreestimaciones de 
volumen de madera, en cambio el 
Modelo 8 está por debajo de ese inven- 
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tario. En este caso la diferencia entre el 
inventario por mea (42019,11 ms) y el ob- 
tenido a través del modelo 8 (37804,99 
m3) es de 4214 m3. 

CONCLUSIONES 

La variable dasométrica área basal (as) 
tuvo una correlación positiva con las 
variables cartográficas temperatura 
media anual (tma), altitud sobre el nivel 
medio del mar (altitud) y precipitación 
media anual en lluvias (pmamo). Las 
variables biofísicas de eos, lais y laial 
también mostraron alta correlación con 
AB. Destaca la correlación de ab con laial, 
esto se explica por la relación de AB con 
la biomasa seca del follaje, misma que 
fue parte importante para la estimación 
de LAIAL. La correlación de ab con la 
banda 4 (S4) siempre fue negativa. Las 
estructuras foliares de ios árboles influ- 
yeron en la relación de ab con b4. 


Tabla 7. Comparación de las estimaciones de volumen de madera (m3) obtenidas 
mediante un inventario tradicional y con inventario por sig. 

Intervalos de confianza (ic) a 95% de significancia estadística. 


fhaf'""'® 298,66 

Volumen de Madera Total 

Límite Límite _ 

3 ,. • í • Precisión 

m /ha TOTAL inferior superior o/ /di \ 

del IC del IC 

Inventario 

Tradicional 

Inventario MEA 
Inventario MRR 

140,69 

118,89 

42019,11 

35507,37 

39800,52 

33933,19 

44237,7 

37081,55 

5,28% 

4,43% 

Inventario 

SIG 

Modelo 5 (COB) 
Modelo 6 (LAI5) 
Modelo 7 (LAIAL) 

Modelo 8 
(Cartográficas) 

216,54 

212,9 

233,11 

126,58 

64672,88 

63585,74 

69621,76 

37804,99 



MEA = Muestreo 

COB = Cobertura arbórea estratificad oaleatorio 

MRR = Muestreo de razón y 
LAI5 = índice de área foliar en sitio regresión 
LAIAL = índice de área foliar 
alométrico 
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De los cuatro modelos propuestos para 
AB, el Modelo 2 (m/s) es el que mejor se 
aproximó a la estimación de inventario 
total. El inventario obtenido por mea fue 
de 22,54 m2/ha mientras que por Modelo 
2 fue de 20,47 m2/ha, la diferencia en 
área basal fue de 2 m2/ha. Destaca el 
resultado obtenido con el Modelo 4, el 
cual está integrado únicamente por varia- 
bles cartográficas. A pesar de tener la 
menor R2 ajustada (0,7117) fue el 
segundo modelo que más se aproximó al 
inventario final (16,09 m2/ha). 

Para la variable dasométrica volu- 
men de madera (vol) también se observó 
una correlación positiva con las variables 
TMA, ALTITUD y PMAMO. Las Variables biofí- 
sicas de coe, lai5 y laial también 
mostraron alta correlación con vol. 
Destaca la relación de vol con savi. 

De los cuatro modelos propuestos 
para vol, el Modelo 8 (que contiene varia- 
bles cartográficas) fue el que mejor se 
aproximó a la estimación obtenida 
mediante el inventario tradicional. Al 
tomar como referencia el inventario del 
predio obtenido por mrr (118,89 m3/ha) 
se observa que tal estimación es la más 
cercana al inventario obtenido por el 
Modelo 8 (126,5 m3/ha) y la única que no 
sobrestima volumen total. 

El mayor problema enfrentado en la 
modelación espacial del ab y el vol fue el 
empleo de información de escala muy 
pequeña (1: 250 000). Bajo estas circuns- 
tancias, por ejemplo, la interpolación de 
las isotermas e isoyetas fue difícil, apare- 
cieron amplias zonas con poca variación, 
razón que afectó a los modelos. En este 
sentido, puede afirmarse que de contar 
con información fisiográfica y climática de 
escala grande (1:25000) se habrían obte- 
nido resultados altamente significativos 
(Silbernagel y Moeur, 2001). De hecho, 
aun trabajando con información geográ- 
fica de poca calidad, se obtuvieron coefi- 


cientes de determinación múltiple mode- 
radamente altos. Razón que justifica el 
empleo de variables cartográficas en el 
inventariado de recursos forestales. 
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